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The  results  of  the  investigation  indicated  that,  for  simple 
span  bridges,  the  amplitudes  of  accelerations  which  psychologically 
disturbed  the  pedestrian  were  predominately  affected  by  the  span 
length  of  the  bridge,  the  weight  and  speed  of  the  vehicle,  and 
the  surface  roughness.   The  parameters  that  provided  the  minor 
effects  were  the  girder  flexibilities  and  the  transverse  position 
of  load.   For  two  and  three  span  continuous  bridges,  the  magnitudes 
of  acceleration  were  larger  than  the  recommended  limit  of  comfort 
only  when  the  surface  roughness  of  the  bridges  was  taken  into 
account. 

Deflection  limitations  and  maximum  girder  span-depth  ratios 
used  in  present  bridge  design  codes  do  not  assure  the  bridge 
users'  comfort.   Present  bridge  design  criteria  have  limited  the 
use  of  high  strength  steel  although  the  dynamic  characteristics 
of  the  bridge  may  be  satisfactory.   Some  reports  indicate  that 
bridge  acceleration  was  significant  in  producing  psychological 
effects.   The  results  of  this  study  indicate  high  strength  steel 
girders  could  be  used  for  highway  bridges  since  the  effect  on 
bridge  accelerations  of  flexibility  of  the  steel  was  relatively 
Insignif 1  cant. 
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CHAPTER  I 
INTRODUCTION 

The  purpose  of  this  investigation  is  to  develop  concepts 
and  simple  approximate  relations  that  may  be  used  to  estimate 
the  magnitude  of  the  dynamic  response  of  high^^;ay  bridges 
having  significance  in  producing  adverse  psychological  effects 
to  bridge  users.   The  investigation  is  aimed  at  a  better 
understanding  of  bridge  vibrations  and  ultimately  establish- 
ing a  basis  for  inclusion  in  bridge  specifications  of 
design  criteria  which  will  directly  regulate  dynamic  response 
characteristics. 

Dynamic  response  is  not  specifically  referenced  in 

(1)  * 
l^rcsont  bridge  design  codes     .   Deflection  limitations  and 

maximum  girder  span-depth  ratios  have  been  used  in  the  hope 
that  such  limitations  will  lead  to  structures  that  have 
satisfactory  vibration  characteristics.   Those  code  pro- 
visions do  not  attack  the  problem  directly  and  the  results 
obtained  by  their  observance  seem  to  reflect  this  fact. 
It  is  probable  that  other  phenomena,  for  example  acceleration, 
at  c  at  least  a.s  significant  as  deflection  amplitudes  in 


* 

Numbers  in  parentheses,  unless  otherv;ise  identified,  are 
reference  numbers. 


producing  adverse  psychological  effects.   Complaints  of 
disturbing  bridge  motions  have  come  from  pedestrians  and 
from  persons  in  halted  vehicles. 

The  economical  use  of  modern  high-strength  steels  in 
bridge  construction  has  been  hindered  somewhat  by  present 
code  limitations.   In  order  to  obtain  stiffnesses  which  meet 
AASHO  deflection  and/or  span-depth  ratios  an  increase  in 
the  quantity  of  steel  is  required.   The  increment  of  human 
comfort  thus  gained  is  of  questionable  significance  for 
bridges  not  located  in  urban  areas,  i.e.,  for  bridges  not 
serving  pedestrians  or  stopped  vehicles.   Significant 
savings  would  result  if  the  design  of  high  strength  steel 
girders  could  be  freed  from  arbitrary  restrictions  which 
do  not  guarantee  relevant  improvement  in  bridge  performance. 

Although  tlie  response  of  simple  span  and  multi-span 
continuous  highway  bridges  to  the  influence  of  moving 
vehicles  has  been  studied  at  some  length        ,  there 
has  been  comparatively  little  attention  given  to  the  control 
of  dynamic  response  of  bridge  structures  which  will  assure 
human  comfort.   The  design  of  bridge  structures  at  the 
present  time,  which  is  based  on  deflection  limitations  and 
limiting  span-depth  ratio,  still  does  not  guarantee  satis- 
factory vibration  characteristics  within  acceptable  limits, 
since  the  proper  parameters  influencing  the  vibration  have 
not  been  recognized  appropriately. 


The  dynamic  response  of  three  kinds  of  highway  bridges 
will  be  investigated  in  this  study:   the  simple  span, 
two  span  continuous  and  three  span  continuous  highway 
bridges.   It  is  necessary  that  the  method  of  analysis  which 
is  required  to  determine  the  response  of  these  bridges 
should  allow  consideration  of  the  major  parameters  which 

affect  the  response.   The  most  successful  analytic  studies 

(28) 
in  recent  years  have  been  presented  by  Oran     for  simple 

(4) (29) 
span  highway  bridges  and  by  Veletsos  and  Huang        for 

multi-span  continuous  highway  bridges.   Parameters  which  may 

be  considered  in  the  analysis  include  bridge  damping, 

roadway  unevenness  and  initial  vertical  oscillation  of 

the  vehicle.   The  theory  was  verified  by  comparisons  with 

the  results  of  laboratory  studies  on  simply  supported 

beams  and  with  results  of  the  AASHO  Road  Test  bridge 

(22) 
studies.       The  bridges  tested  in  the  AASHO  study  were 

simple-span  single  lane  structures  conforming  rather  closely 

to  the  assumptions  of  the  theory. 

There  have  been  a  few  reports  available 

discussing  criteria  for  human  response  to  vertical  vibra- 

(5) 
tion.   A  recent  paper  by  Wright  and  Walker     summarizes 

the  effects  of  bridge  flexibility  on  human  response  and 
presents  a  summary  and  discussion  of  the  previous  studies 
of  the  criteria.   It  is  suggested  in  the  paper  that  accelera- 
tion of  the  bridge  is  more  important  than  the  velocity  or 
jerk  in  producing  psychological  effects  since  the  natural 


frequencies  of  highway  bridges  range  between  1  to  10  cps . 
The  report  further  suggests  that  the  amplitude  of  the 

dynamic  component  of  acceleration  in  the  fundamental  mode 

2 

of  vibration  should  be  limited  to  100  in/sec   to  assure 

the  human  comfort.   Acceleration  criteria  for  human  response 
presented  by  the  paper  are  shown  in  Table  1.1. 

The  primary  purpose  of  this  investigation  is  to  study 
the  behavior  of  highway  bridges  under  a  moving  vehicle, 
as  affected  by  major  parameters  and  to  obtain  the  maximum 
accelerations.   By  comparison  with  the  available  acceleration 
criteria,  conclusions  about  the  significance  of  each  parameter 
can  be  made  and  finally,  the  accelerations  of  given  bridge 
structures  can  be  predicted  analytically.   All  parameters 
considered  in  this  study  are  kept  within  a  realistic  range 
of  expected  values. 

The  method  of  analysis  and  a  description  of  computer 
programs  to  compute  the  response  of  simple  span  and  multi- 
span  highv/ay  bridges  are  presented  in  Chapter  II.   All 
computer  programs  used  in  this  study  were  originally 
developed  at  the  University  of  Illinois  and  contributed 
by  Prof.  VI.    H.  Walker  of  the  Civil  Engineering  Department 
of  the  University  of  Illinois.   His  contribution  is  gratefully 
acknowledged.   The  programs  are  modified  to  handle  dimen- 
sional input  parameters  and  to  compute  the  accelerations 
at  any  specified  point  on  the  bridges.   In  Chapter  III  the 
effects  of  major  parameters  on  accelerations  of  multi-girder 


Table  1.1   Acceleration  Criteria  for  Human  Response  to 
Harmonic  Vertical  Vibration ^^^ 


Human  Response 


5- 

Amplitude  of  Acceleration  (in/sec  ] 

Transient         Sustained 


Imperceptible 
Perceptible  to  few 
Perceptible  to  some 
Perceptible 
Unpleasant  to  few 
Unpleasant  to  some 
Unpleasant 
Intolerable  to  few 
Intolerable  to  some 
Intolerable 
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simple  span  bridges  are  studied.   Chapter  IV  is  devoted 
to  a  study  of  the  accelerations  of  three  span  continuous 
highway  bridges  as  affected  by  major  parameters  and 
the  initial  conditions  of  the  vehicle.   The  effect  of  the 
same  parameters  are  studied  in  Chapter  V  for  two  span 
continuous  bridges.   A  summary  of  the  important  findings 
and  conclusions  are  given  in  Chapter  VI. 


CHAPTER  II 
METHODS  OF  ANALYSIS  AMD  COMPUTER  PROGRAMS 

2 . 1   General 

The  purpose  of  this  chapter  is  to  review  the  analyses 
of  simple  span  and  of  continuous  highway  bridges  under 
moving  vehicles.   General  theories  for  both  analyses  have 
been  summarized  and  the  methods  to  compute  the  accelera- 
tions at  any  specified  points  on  the  bridges  have  been 

added.   The  methods  of  analysis  for  simple  span  and  con- 

( 28 ) 

tinuous  highway  bridges  were  developed  by  Oran     and 

(29 ) 
Huang     ,  respectively.   The  descriptions  of  computer 

programs  for  both  analyses  are  also  presented  in  this 

(28) (29) 
chapter.   Both  original  computer  programs         have 

boon  modified  to  compute  the  accelerations  of  the  bridges 

in  addition  to  the  original  response  which  includes  the 

static  and  dynamic  deflections,  the  static  and  dynamic 

moments  and  the  interacting  forces  between  the  vehicle 

and  the  bridge.   The  input  and  output  form  of  both  programs 

have  also  been  revised  so  that  the  programs  can  handle 

either  the  dimensional  or  non-dimensional  data  input  and 

provide  printed  and  plotted  results.   Most  of  the  symbols 

used  in  the  analyses  are  defined  where  they  are  first 

introduced . 


2 . 2   Analysis  of  Static  and  Dynamic  Response  of  Simple- 
Span/  Multigirder  Highway  Bridges 

2.2.1   General 

In  this  analysis,  the  bridge  is  represented  as  a 
plate  continuous  over  flexible  beams,  and  the  bending  and 
torsional  stiffnesses  of  the  beams  are  taken  into  account. 
The  analysis  involves  two  major  steps: 

(a)  The  determination  of  the  instantaneous  values 
of  the  dynamic  forces  acting  on  the  bridge;  these  include 
the  interacting  forces  between  the  vehicle  and  the  bridge, 
and  the  inertia  forces  of  the  bridge  itself. 

(b)  The  evaluation  of  the  deflections  and  moments 
produced  in  the  bridge  by  these  forces. 

The  second  step,  which  is  strictly  a  problem  of 
statics,  is  solved  by  an  application  of  the  Rayleigh-Ritz 
energy  procedure.   The  deflection  of  the  structure  is 
expressed  as  a  series  combination  of  functions  that  is 
capable  of  approximating  any  deflection  configuration  in 
both  the  longitudinal  and  transverse  directions. 

The  method  used  to  evaluate  the  dynamic  forces  is 
essentially  an  extension  of  that  used  for  the  static  prob- 
lem.  In  contrast  to  the  expression  used  to  represent  the 
static  deflection,  however,  the  dynamic  deflection  of  the 
bridge  is  assumed  to  be  half-sine  wave  in  the  direction 
of  the  span. 


The  vehicle  is  represented  by  a  single-axle  loading 
consisting  of  a  sprung  mass  and  two  equal  unsprung  masses, 
or  wheels.   The  so-called  rolling  effect  of  the  vehicle 
is  thus  taken  into  account.   The  springs  are  assumed  to 
be  linearly  elastic.   No  damping  is  considered  for  either 
the  vehicle  or  the  bridge. 

2.2.2   Analysis  of  Static  Problem 

The  structure  considered  is  shown  in  Figure  2.1. 
It  consists  of  a  reinforced  concrete  slab  continuous  over 

parallel  steel  or  reinforced  concrete  beams  spanning 
in  the  direction  of  traffic  and  simply  supported  at  the 
ends.   The  beam  spacing  may  be  arbitrary.   The  dimensions 
of  the^  beams  may  vary  from  one  beam  to  the  next,  but  all 
beams  ,ire  assumed  to  be  prismatic.   The  slab  is  considered 
to  be  isotropic,  of  constant  thickness,  and  simply 
supported  at  the  abutments. 

The  assumptions  made  in  the  analysis  are  those 
embodied  in  the  ordinary  theory  of  medium-thick,  elastic 
plates  and  in  the  ordinary  theory  of  flexure  of  beams. 
In  addition,  it  is  assumed  that 

(1)  A  beam  and  the  slab  over  it  deflect  and  rotate 
alike . 

(2)  There  is  no  transfer  of  horizontal  shear  between 
the  beams  and  the  slab. 
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Figure  2.1    Simple  Span  Bridge  Model 
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The  span  length  of  the  bridge,  center  to  center  of 
supports,  is  denoted  by  L,  and  the  overall  width  of  the 
slab  by  b.   The  position  of  a  point  on  the  bridge  is 
specified  in  terms  of  dimensionless  cartesian  coordinates 
K    and  n,  defined  by  the  equations 


C  =  x/L 

(2-1) 
n  =  y/b 


Method  of  Analysis 

The  method  used  is  a  combination  of  the  Rayleigh-Ritz 
energy  method,  and  the  Levy  method  of  analysis  for  rectantu- 
lar  plates  simply  supported  along  two  opposite  edges.   Let 
the  vertical  load  on  the  structure,  p(f;,ri),  be  represented 
by  a  single  trigonometric  series  of  the  form 


PiK>T))    =      ^  p^  =  I     p^sinm/[^      •         (2-2) 
m=l      m=^l 


in  which  p   is  a  function  of  n  only.   The  deflection  of 

'^m  -^ 

the    structure,    w(5,ri),    can    then   be   expressed   as 


w(f;»n)    =      y     w       =      y     w     smmiif,  (2-3) 

'- ,     m  '" ,     m 

m=l  m=l 


where  v;   =  w  sinmnf;  is  the  deflection  component  correspond- 

m    m  -^ 

inq  to  the  load  component  p  ,  and  w   is  a  function  of  n 
"^       m       m 

only.   Let  the  deflection  functions  w   be  expressed  in  the 

form 


12 


w   =  y  ex   Y 
ra    ^  mn  n 


(2-4) 


where  Y   are  known  functions  of  the  n-coordinate,  and  a 

n  mn 

are  coefficients  which  will  be  evaluated  by  minimizing 

the  total  energy  of  the  system.   Let  the  functions  Y   be 

dimensionless ;  then  the  coefficients  a    have  the  dimension 

mn 

of  length. 

Let  the  Y   functions  be  taken  as  follows: 


n 


1      for  n  =  -1 
0  .5  -  n    for  n  =^  0 
sin  nn  n   for  n  2:  i 


(2-5) 


Note  that,  Y_,  ,  Y, ,  Y,,  ...  are  symmetric  and  Y-,  Y_,  Y .  ,  ... 

are  antisymmetric  with  respect  to  the  longitudinal  center- 
line  of  the  structure,  n  =  0.5. 

Energy  of  the  System 

The  energy  of  the  system  consists  of  the  strain 

energy  of  the  slab,  the  strain  energy  of  the  beams  and  the 

potential  energy  of  external  forces.   The  energy  expressions 

correspond  to  the  deflection  component  v;   =  w  sinm-nC 
*■  "^        m     m       ^ 

and  the  c'lssociated  load  component  p   =  p  sinmirC. 

'^      ^m    "^m 

Let  the  strain  energy  of  the  system,  v  ,  be  written 
in  the  form 


V  =  (V  )   +  (V  ). 
ro     m  s     m'b 


(2-6) 
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where  (V  )^  is  the  strain  energy  of  the  slab,  and  (V  ),  is 
m  s  ^-^  m  b 

the  strain  energy  of  the  beams.   (V  )  .,  and  (V  )^  are  qiven 


m    3 


by    the    equation 
DLb    rl   fl 


(V„ 


m    s 


2L       -"O    -0 


4    4 
m    It    D   bL 


m   b 


2-\2 


u«")  '~AhJ  "7LJ] 


y    y  a       a  Y 

^   ^     mn   ms  Jq      n 


Y    dn 
s 


dcdn       (2-7) 


-^~ rj    y    y  a       a  Y   'Y   '  dn 

2  ,       >  2    '■    '-     mn   ms   L      n     s 
IT     (iTic)       n   s  •'0 


+   -T-^^ T    y    y  a      rx  Y  "Y  "  dn 

4  ,       ,  4    '■    '-     mn   ms   L      n     s 
7T     (mc)       n   s  -"O 


(V^)h    -       I 


m  b 


i  =  0 


(E.  I.  )  .     rlfd^vj\^                (G^Jv  )  •    rl/a^w^  ^ 
i^   b   X    U        ^\     ^^     ,        b   b   1  [   ( m,      ^^^ 


3  \        2 

2L  ■'0    VH 


2Lb         -"O  \d  r,  3  n/  . 


(2-8) 


(2-9) 


4    4 
m    Ti    D   bL 


i=0 


A  .     y     y  a       a  (Y     )  .     (Y    )  . 

1    ''    '-■     mn   ms       n    1       si 


n   s 


+    -~, ^    y    y  a       a         (Y   '  )  .    (Y   •  )  . 

2,      .2^^mnmsni       si 
n     (mc)       n    s 


(2-10) 


v/here           c  =  b/L 

D  =  the  flexural  rigidity  per  unit  width  of  the  slab 

(E.I,).  =  the  flexural  rigidity  of  the  i^^  beam 

b  b  1  :3       J 

(G,  J,  )  .  =^  the  torsional  rigidity  of  the  i^"  beam 
b  b  1                    ■"           ■' 

and  a  prime  denotes  and  differentiation  with  respect  to  n. 


.lis  .,  <  •  !„i;,  c.i:. 
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The  strain  energy  of  the  system,  V  ,  can  be  rewritten 


in  the  form. 


4  4 
,,    m  TT  D  bL  —  ,       , 

V T —  — r-  I   (a   a   ) 

m      4    4   m   mn  ms 


(2-11) 


where  F„  (a   a   )  is  a   quadratic  form  of  the  unknown 
m   mn  ms 

coefficients  m       . 

mn 

The  potential  energy,  U  ,  of  the  load  component  p 

m  "^      '^m 

is  given  by  the  equation 


U     =   -bL    I  w  p   d^dn 


(2-12) 


=    -bL 


(w     sinmrr^)  (p     sinmirr)  d^dn 
m  m 


(2-13) 


bL 
2 


y  a         I      p    Y    dti 


(2-14) 


The  total  energy  of  the  system,  I  ,  is  the  sum  of 

m 

Eqs.  (2-11)  and  (2-14) 


4  4^  ,,  ,,        fl 

T     mTiDLb-=,  »    Lbr          „, 

4    4   m  mn  ms     2  '-  mn  L   m  n 

L  n     -^0 


m 


(2-15: 


The  condition  that  I   be  a  minimum  yields  a  system 
of  linear  algebraic  equations  of  the  form 


[F„] [A^]  -  [B^] 
mm      m 


(2-16) 
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in  which  [A  ]  is  a  column  matrix  of  the  unknown  coefficients 

m 

a     [B  ]  is  a  column  matrix  of  known  load  terms,  and 

mn,    m 

[F  ]  is  a  symmetric  matrix,  the  order  of  which  is  equal 

to  the  number  of  Y   functions  used  in  Eq .  (2-4).   The 
n  I   V    / 

element  f    ^  in  the  n^^^  row  and  s^^  column  of  the  [F  ] 
m  ,ns  m 


matrix  is  given  by  the  equation 


m,  ns 


Y    y     dn    +   —^ T-         Y  'Y  •  dn    +    —r-— j         Y  "Y  "dn 

„      n    s  2  ,       ,  2     L      n     s  4  ,       .  4    L      n     s 

0  IT     (mc)       •'0  IT     (mc         -"O 


P 

+    y 

i  =  0 


X  .  (Y    )  .  (Y    )  .     +   -TT 
1      nisi  2 


k. 

1 


71     (mc) 


^     (Y   •)  .    (Y^-)  . 
2         nisi 


(2-17) 


The  element  B    in  the  n    row  of  the  [B  ]  matrix  is 

mn  m 


2PL-^  sin  mTT?^ 


'mn     4^,,      4 
TT  Db    m 


^n^^l^ 


(2-18) 


The  solution  of  the  system  of  Equations  (2-16)  gives 

the  value  of  a   ,  which  are  then  used  to  determine  the 
mn 

deflection  component  w   =  w  sinmn^.   In  general,  Eqs . 

mm 

(2-16)  are  solved  for  as  many  values  of  m  as  may  be 
necessary  in  a  particular  application.   The  total  deflec- 
tion, w,  of  the  structure  is  then  determined  from  Eq .  (2-3) 
by  superposing  the  component  deflections.   The  latter 
equation  may  be  rewritten  in  one  of  the  following  forms 


w   =    )  [)  a      Y    1  sinm7i^=    V   [a    ]  •  [Y    (ti)]  sinmnf  (2-19! 

I  ■'     mn    n;  ^       m  n 

m  \n  /  m 
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where  Y  (n)  is  a  column  matrix  of  the  values  of  the  Y 
n  n 

functions  evaluated  at  the  point  under  consideration,  and 
a  dot  denotes  a  scalar  product. 

2.2.3   Analysis  of  Dynamic  Problem 

The  structure  analyzed  is  the  same  as  that  considered  in 
the  static  analysis.   In  addition  to  the  assumptions  made 
previously,  it  is  assumed  that  the  mass  of  the  slab  is  uni- 
formly distributed,  and  that  the  mass  per  unit  of  length  of 
the  beams,  although  it  may  vary  from  one  beam  to  the  next, 
is  constant  for  any  one  beam. 

The  vehicle  is  represented  by  a  single-axle,  two  wheel 
loading  consisting  of  a  sprung  mass  and  two  equal  unsprung 
masses,  as  shown  in  Figure  2.2.   The  center  of  gravity 
of  the  sprung  mass  is  assumed  to  be  located  halfway  between 
the  supporting  springs.   The  springs  are  considered  to  be 
linearly  elastic  and  to  have  identical  stiffnesses.   Damping 
for  both  the  vehicle  and  the  bridge  has  been  neglected. 

The  analysis  of  the  problem  involves 

(a)  The  determination  of  the  instantaneous  values  of  the 
interacting  forces  between  the  vehicle  and  the  structure,  and 
of  the  inertia  forces  due  to  the  mass  of  the  structure,  and 

(b)  The  computation  of  the  deflections  and  bending 
moments  produced  in  the  structure  by  these  forces. 

The  dynamic  deflection  configuration  of  the  structure 
in  the  longitudinal  direction  (^-direction)  is  assumed 
to  be  a  half-sine  wave  and  is  expressed  as 
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Static    equilibrium    position 
of   C.G.  of    sprung    mass 

Y 


^"'"l^i     '"^"l^J'l   m 


(a) 


(b) 


Figure  2.2    Vehicle  Model  for  Simple  Span  Bridge 


w  =  w  Sinn?  Yf  (t)Y  (n)  (2-20) 

o       '■  n     n 
n 

where    w  =  the  deflection  of  any  point  of  the  bridge  at 

any  time,  due  to  the  static  and  dynamic 

effects  of  the  vehicle 

w  =  a   quantity  with  the  dimension  of  deflection, 

chosen  arbitrarily  as  WL  /(E.I,). 

b  b  lo 

W  =  total  static  weight  of  the  vehicle 

(E,  I,  )  .   "  flexural  rigidity  of  a  reference  beam 
b  b  lo  -^  ^ 

f  (t)  -  dimensionless  coefficients  that  are  functions 
of  time;  these  are  the  generalized  coordinates 
for  the  bridge 

Y  (n)  =  dimensionless  functions  of  n,  as  previously 
discussed. 

The  coordinates  used  to  specify  the  configuration  of 
the  sprung  mass  are  the  vertical  displacement  of  the  center 
of  gravity  of  the  mass,  z,  and  the  rotation  of  the  mass 
about  an  axis  normal  to  the  transverse  vertical  plane, 
u,  (see  Figure  2.2).   The  vertical  positions  of  the  unsprung 
masses  are  determined  by  the  configuration  of  the  bridge. 
Thus  the  total  number  of  generalized  coordinates  of  the 
bridge-vehicle  system  is  equal  to  the  number  of  f  (t) 
functions  used  in  Eq .  (2-20)  plus  the  two  coordinates  z 
and  u  used  for  the  vehicle. 

The  vehicle  is  assumed  to  move  along  the  bridge 
v/ith  a  constant  velocity,  v,  and  the  wlieels  of  the  vehicle 


■;'  J  ?  '    i  .',    .I':o  .'i.i''.'-  _,.'w 


•■•*•      VT\i^     .iO      ^;i".J 
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are  assumed  to  be  attached  to  the  surface  of  the  bridge 
at  all  times.   The  system  under  consideration  possesses 
a  time-dependent  potential  energy  function,  and  it  is 
possible  to  formulate  the  equations  of  motion  by  applica- 
tion of  Lagrange's  equation 


JT_  ^  MU+V)  .  0  (2_2l) 


in  which     V  =  the  strain  energy  of  the  system 

U  =  the  potential  energy  of  the  gravity  forces 
T  =  the  kinetic  energy  of  the  system 

q   =  the  n^"  generalized  coordinate  of  the  system 

^n  ^  -^ 

■   .   __  ^ 

'■'n    dt  ■    .   . 

Energy  of  the  System 

The  datum  of  zero  energy  level  for  the  system  is 
defined  by  the  following  conditions:   the  structure  is  in 
an  unstressed  position,  and  the  springs  of  the  vehicle  are 
unde formed. 

Let  VI      represent  the  deflection  configuration 
of  the  bridge  when  loaded  with  its  own  weight:   this 
deflection  is  measured  from  the  unstressed  configuration 
of  the  bridge.   Then  the  total  deflection  of  the  bridge, 
measured  from  its  unstressed  position,  is  (w  +  w»).   The 
dead  load  deflection  configuration,  w-,  can  be  represented 
in  a  form 
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'      =   w      y    y  6       Y    1  sinimrC  (2-22) 

2  o  ^If:    mn    n/ 


-   w    sin  Tir    y  5    Y    (n)  (2-23) 

o  '      n    n 

n 

where  ^        and  6   are  constant  dimensionless  coefficients. 

mn      n 

Let  w^  be  a  deflection  function  representing  the 
deviation  of  the  deck  of  the  bridge,  when  loaded  with  its 
own  weight,  from  the  horizontal  plane  passing  through 
the  supports.   It  is  positive  when  downward.   This  quantity 
is  equal  to  the  sura  of  the  dead  load  deflection  configura- 
tion and  the  configuration  representing  any  possible 
unevenness  of  the  unstressed  bridge.   The  roadway  surface 
unevenness  function,  w  ,  appears  only  as  (w, ) . ,  i.e., 
with  its  ordinates  evaluated  at  the  transverse  location 
of  the  wheels,  can  be  represented  in  a  form 

(w,  )  .  =  w  e.  sinm.Tif  (2-24) 

1  ]     o  3      1  ^ 

v;here  e.  (j  =1  or  2)  are  dimensionless  quantities  and 
w  e.  denotes  the  amplitude  of  the  unevenness.   The  quanti- 
ties e,  and  e^  may  or  may  not  be  equal;  they  may  also 
be  positive  or  negative.   The  quantity  m,  is  a  positive 
integer. 

The  total  strain  energy  of  the  system,  V,  is  written 
in  the  form 


V  +  V,  +  V  (2-25) 

s     b     sp 
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where  V  ,  V,  ,  and  V    are  the  strain  energies  of  the  slab, 
s    b        sp  -=' 

of  the  beams,  and  of  the  springs,  respectively,  and  are 
given  by  the  equations 


'n 


d^dn 
(^-26) 


=  ! 


i-O 


(E,  I,  )  .   rl/9^(w+W^A  ^      (G.J.).   rl/a2(w+W,)V 


3      \     2 
2L     •'O  V   8f/ 


2Lb^   Jo  \   9^3n 


(2-27) 


sp 


=  il^  y   f2  +  z  -  (w+w,  )  .  +  (-l)^ub,  ] 
2   ■  -1       S        i  1  i 


(2-28) 


j  =  l 


The  complete  forms  of  the  equations,  after  substituting 

the  values  of  w,  Ww  w„  and  evaluating  the  integrals, 

are  given  in  Ref.  (28). 

The  potential  energy  of  the  gravity  forces  is 

v/ritten   in    the    form 

U   =    U+U,     +U        +U  (2-29) 

s  b  sp  u 

whore  U  ,  U,  ,  U   ,  and  U   are  the  potential  energies  of 

s    b    sp       u         ^  ■' 

the  slab,  of  the  beams,  of  the  sprung  mass,  and  of  the 
unsprung  masses,  respectively,  and  are  given  by  the 

equations 

•1  fl 

(2-30) 


U  =  -ugLb  I   I   (w  +  w-)  d^dn 
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"b  =  -^^ 


i=0 


(m.  )  .  (w  +  w_)  . 

D     X  ^  1 


dC 


U   ^  -Mg  (z  +  z  ) 
sp     -^      s 


(2-31) 


(2-32) 


U   =  -mg  I      (w  +  w,  ) 
j  =  l 


(2-33) 


The  complete  forms  of  the  equations,  after  substituting 
the  values  of  w,  v;,  ,  w_  and  evaluating  the  integrals, 
are  given  in  Ref.  (28). 

The  total  kinetic  energy  of  the  system  is  expressed 
in  the  form 


T  =  T   +T,  +T    +T 
s     b     sp     u 


(2-34) 


whore  T  ,  T,  ,  T   ,  and  T   arc  the  kinetic  energies  of 
s    b    sp       u  ^ 

the  slab,  of  the  beams,  of  the  sprung  mass,  and  of  the 
unsprung  masses,  respectively,  and  are  given  in  the 
forms 

T^  -  ^pLb  f   f  (^)      d5dn 
s    2    Jq  Jq   at 


(2-35) 


!   f\m.),  (|^).'  dC 


(2-36) 


„,     1  „  /dz>   ,  1  T  ^du> 
%=  2^  ^dt^  "■    2^   ^dt^ 


(2-37) 


•u  =  i"^  J,  ^^^"■'^'i^j^ 

j  =  l 


(2-38) 


"     i:<-- 
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The  complete  forms  of  the  equations,  after  substituting 
the  values  of  w,  w, ,  and  evaluating  the  integrals,  are 
given  in  Ref.  (28). 

The  Governing  Differential  Equations 

The  differential  equations  governing  the  motion  of 
the  bridge-vehicle  system  are  obtained  from  Eq .  (2-21) 
by  substituting  the  energy  expressions  derived  above.   The 
number  of  equations  obtained  is  equal  to  the  number  of 
generalized  coordinates  used  to  define  the  configuration 
of  the  system.   The  final  derived  equation  corresponding 
to  the  n    generalized  coordinate  of  the  bridge,  f  ,  is 


V7       }     f  " 
O  •'       s 


iiLb 
2 


Y  Y  d.,  +^     y    (m,  )  .    (Y  )  .   (Y   )  .  +mr,in^^     y    (Y  )  .  (Y 
L     n  s         2    .f_Q    b'l      n  1      si  L     -^t     ri  J 


s  J 


+  w  y  f ' 

o  '■    s 


^  <T'    .^i<*n':«s'j 


.       TlVt  llVt 

Sin  — ^ —  cos  -r — 


+  w^  y  f 

o  '-    s 
s 


ir  D  ab 


2  Tivt 


A      +    y     A  .  (Y   )  .  (Y  )  .  +     y     —^  (Y  ' )  .  (Y  ' )  . 
ns     >-,    1     n  1     s'l      .'-       2  2      n    i     si 

1=0  1=0    TT    c 


2         -1^2 

II V.  •^     .    2  iivt 


+  k  sm    — —    )    (Y  )  .  (  Y  )  .  -m  (tj— )    3in    -r—     )    (Y  )  .  (Y  )  . 
L      >,ni       SI  L  L      >,nisi 

j=l         ■"  ■'  j=l         -■  -' 


sinli^     y    (Y   ), 
^     j=l     "^ 


2^  +  k 


z-  (w^)  •  +  (-l)^ub-^ 
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mg  sm 


77 '/t 


,V>  ■^         .         TTVt 


n2 


I    iVj.+m   (f)     siji^     y    (y  )       -_t      :.o  (2-39) 


L     .-^  ■  n- D 


j=l 


"   T  \;,r2 


The  summations  on  s  should  extend  from  s  =  -1,  to  s  =  n, , 
the  maximum  value  of  n  used  in  Eq .  (2-20).  There  will  be 
a  totcil  of  n.  +  2  equations. 

The  equation  corresponding  to  the  7.   coordinate  of 
the  vehicle  is 


Mz"  +  k     I 
J-1 


ir  vt 


z  -  (w,  )  .  -  w  sin  —T —  y  f  (Y  ) 
1  1    o      L   '^  s   s 


=0   (2-40) 


aiid  the  one  corresponding  to  the  u  coordinate  is 


Ju"  H    kb,      y   (-1)^ 


irvt 


(-])^ub,-(w^).-w^sin-^If^(Y^). 


whore      A 


Y    Y    dn  +  -i^    I    Y   'Y  'dn  +  -i^  I    Y  "Y  "dn 
„ns  22    1-     ns  44l-,ns 

0  IT    c      •'0  Ti    c     •'0 


(2-41) 
(2-42) 


The  dimensionless  form  of  the  governing  equations  can 
Le  obtained  by  assuming  the  following  parameters 


2m  +  M 

P 
L[iib+   ):  ("Ik);] 


Total  weight  of  vehicle 
Total  weight  of  bridge 


i=0 


b'i- 


.'m 


2m  +  M 


(mj^). 


Unsprung  weight  of  vehicle 
Total  weight  of  vehicle 


(2-43) 


h  +     Y.      (m,  ) 
i=-0   D 


Weight  of  i^^  beam 
Total  weight  of  bridge 
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vT, 
b 

'217 


fv2     Tb2 
fb      Tv 


in  which 


f,  =  the  fundamental  natural  frequency  of  the  bridge 

evaluated  on  the  assumption  that  it  acts  as  a  one 
dimensional  beam. 

f  ~    the  natural  frequency  of  the  vehicle  for 
vertical  motion  on  its  springs. 
These  frequences  are  given  by  the  equations 

P 
^  Db  +  Z  (E.l.). 

f  2  ^   _1_  =   _L  i-0   ^  ^  ^ 

^   ~  T  2  ~  4  J   ^    P   ,   , 
b         ub  +  Z   (m,  )  . 

i=^0   ^  ^ 


(2-44) 


f  2  =  _i_  ^  _!_   2k 
•^   ~  T  2  ~  4.2  ^1 

V 


In  addition,  let 


vt 


^^(T)   =  f^(t) 


C(x)  =  ^^ 


(2-45) 


u(t)b, 
G(r)  =  — i 
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Now  by  multiplying  Eq .  (2-39)  by  -^ 


pbL  +  L  E  (m,  )  . 


i  =  0 


b'  i 


and  substituting  the  dimensionless  quantities  defined  by 
Eqs .  (2-43)  and  (2-44),  the  governing  differential 
equations  can  be  reduced  to  the  forms 


y    i|)"    (B    ^  +  C       sin     ttt)    +    y    ijj    '  D       sin  ttt  cos  ttt 
^    ^s        ns        ns  ^      s      ns 

s  s 


+    y    if)^(E    ^  +  F    ^sin     iit)+H     sin  71 
„      s       ns        ns  n 


+   L    (?,(), t)    +   R    siniTTsinm,   nr    =    0 
n  n  1 


(2-46) 


2a '"c"    +    ■"    K 


2c  -  s 


inm,TTT    I    e.-sinTTiJ^ij;        I  (Y    ). 

■^         1=1     ^  5^1  =  1^=" 


2a^p0"    +    Ti'^K 
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-  sin  m,  TTT      I      (-1) 


-^^:^, 


j=l 


-    sin  T,T    I    ^^   J^  (-l)J  (Y^). 


s  j  =  l 


=    0 


=  0 
(2-47) 


(2-48) 


where 


ns 


(l  -      y    Y-l   1      Y    Y    d     4      f    y  (Y    )  .  (Y    )  . 
V        iio    7J0      "    ^    ^       j  =  0^      nisi 


ns 


a^vo)       I      (Y^).(Y^). 
•     1        n    i       ST 


D         =    2a    TTvu       y      (Y    )  .  (Y    ) 


ns 


j=l 


n    D      s    J 


,  i  '    1  C  ^' '■    -»".',     •  '■*:-"• 
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ns 


1  +     z     \. 
i=0    ^ 


A         +       )      A  .     Y       .  (Y        . 
ns         .  '•  -     1      n    1      s    1 
1=0 


^  -T-2    .L'^i^V^i^V^ 


ns 


TT     c       i=0 


tt^v[k  (1  -  u))   -  a^u)]        ?     (Y    )   .  (Y     )   . 
j4l       "    3       s    J 


H       =    - 


lo 


1  +      E      X  . 
i  =  0     ^ 


j  =  l        '^    ^ 


L      =-TrKv(l-co)sinTTT       ?(Y).[f,+    (-l)^0] 


2  2         2         r 

R      =    n     v[K(l-a))     -    a    torn,]      )      e.(Y    ). 


j=l 


p    - 


Mb, 


Computation  of  Response 

The  procediire  used  to  evaluate  the  dynamic  response 
of  the  bridge-vehicle  system  may  be  summarized  briefly 
as  follows: 

(a)  The  governing  differential  equations  of  motion 
are  solved,  by  means  of  a  step-by-step  method  of  numerical 
integration,  to  determine  the  values  of  the  generalized 
coordinates  and  of  their  first  two  derivatives. 
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(b)  The  interacting  forces  between  the  vehicle  and 
the  bridge  and  the  inertia  forces  of  the  bridge  are 
evaluated. 

(c)  The  dynamic  deflections  and  bending  moments 
induced  in  the  bridge  are  determined  from  the  dynamic 
forces  acting  on  the  bridge,  instead  of  directly  from  the 
generalized  coordinates  computed  in  the  first  step,  and  the 
accelerations  of  the  bridge  are  evaluated. 

The  system  of  Eqs.  (2-46)  through  (2-48)  are  solved 

by  means  of  a  step-by-step  method  of  numerical  integration. 

The  time  required  for  the  vehicle  to  cross  the  span, 

0  <  T  <  1,  is  divided  into  a  number  of  small  intervals,  and 

tlie  governing  equations  are  "satisfied"  only  at  the  ends 

of  tliese  intervals. 

Let  q     represent  a  dimensionless  generalized 
'n,r    "^  ^ 

coordinate  -  it  may  refer  to  the  bridge  or  the  vehicle  -, 

e\nd  q     and  q     represent  its  first  and  second  derivatives 
^n,r      n,r    ^ 

v/ith  respect  to  t  .   Let  it  be  assumed  that  the  values  of 
r 

q    ,  q     and  q     are  known  for  each  generalized  coordinate 
^n,r   ^n,r      ^n,i:  ^ 

of  the  system,  and  that  it  is  desired  to  find  the  corres- 
ponding values  at  T  =  T  , ,  =  T   +  Ax  in  which  Ax  is  a 
^      ^  r+1     r 

short  interval. 

The  variation  of  q   within  the  time  interval  Ax  was 
■"n 

considered  to  be  linear;  with  this  assumption,  the  express- 
ions for  q    , ,  and  q    , ,  become 
■'nfr+l      ■'n,r+l 
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'^rirr+l    ^n,r    2      n,r   "n,r  +  . 


I    ,  -,  =  q     +  (At)  q      +  i(AT)^q     +  i(Ar)^q 


+  1 


(2-49) 


(2-50) 


The  iterative  procedure  may  be  suininarized  as  follows: 

1.  Assume  that  the  second  derivatives  of  the  general- 
ized coordinates  at  the  end  of  the  time  interval  are  the 
same  as  those  at  the  beginning  of  the  interval,  i.e.  take 

q    ,,  --  q     and  by  application  of  Eqs.  (2-49)  and  (2-50) 
""n/r+l     n,r      j      fi  i  ' 

evaluate  q    , ,  and  q    , ,  . 
^n,r4l      ^n,r-i-l 

2.  Substitute  the  values  of  q    ,,  and  q    ,,  thus 

^n,r+l      ^n,r+l 

obtained  into  the  governing  differential  equations,  and 

by  solving  the  resulting  system  of  algebraic  equations, 

obtain  improved  values  for  q    ,  ^  . 
'■  ^n,r+l 

3.  From  Eqs.  (2-49)  and  (2-50)  calculate  the  values 

of  q    , ,  and  q     ^  corresponding  to  the  values  of  q    , , 
'n,r  +  l      ^n,r+l       '^  ^  ^n,r+l 

just  determined. 

4.  Repeat  Step  2  by  using  the  latest  available 

values  of  q    , ,  and  q    , ,  . 
^n,r+l      'n,r+l 

5.  For  each  generalized  coordinate  compare  the  newly 
derived  values  of  q    ,,  with  the  previously  available  value. 
If  the  difference  between  the  two  values  for  each  coordinate 
exceeds  a  prescribed  tolerance,  repeat  Steps  3  through  5,  un- 
til all  differences  are  less  than  the  prescribed  tolerance. 
The  algebraic  equations  are  then  considered  to  be  solved,  and 

the  integration  for  the  time  interval  r  to  t  ,  completed. 

r     r+1    '^ 
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The  next  step,  to  obtain  the  dynamic  response  of  the 
system,  is  to  evaluate  the  interacting  forces  between  the 
veliicle  and  the  bridge.   The  static  value  of  the  interacting 
force  for  each  wheel  is  one-half  the  total  weight  of  the 
vehicle  or  W/2.   The  dynamic  increment  of  the  interacting 

force  for  the  j ^"  wheel  may  conveniently  be  stated  in  the 

W 
form  A.  -jr  where  A  ■  is  a  dimensionless  factor.   Now  let 
D  2^         1 


A.  =  (A^).  +  (A^)^ 


(2-51) 


where    (A 


I'j 


(A^). 


the  component  of  A .  due  to  the  dynamic 
increment  of  the  compression  in  the  spring 
the  component  of  A .  due  to  the  vertical 
acceleration  of  the  unsprung  mass  j . 


Those  quantities  may  be  determined  as  follows: 


W 
I'j  2 


(Aj^;^=  k[z  +  (-1)^  ub,  -  (^1+ w)  .] 


(A^)  .  2 


(2-52) 


(2-53) 


The  dynamic  increments  of  deflection  (AD) .  and  of 
moment  (AM) .  can  now  be  expressed  in  the  form 


(AD) .     =       y   D.   .  A .  +  y  D?    ^  " 
(AM)  .  =   y  M.   .  A  .  +  y  M°    V  " 


(2-54) 


(2-55) 
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where    D.  .  "  the  deflection  produced  at  a  specified  point 

1,3  '^  ^ 

of  the  i^"  beam  by  a  concentrated  force 

W/2  applied  at  the  position  of  the  j^"  wheel. 

M.  .  =:  the  bending  moment  corresponding  to  D.  .. 
ifD  ^  ^  x,3 

D?    =  the  deflection  produced  at  a  given  point 

1,8 

of  beam  i  by  a  static  load  which  is  distri- 
buted as  a  sine  wave  in  the  longitudinal 
direction. 
M.    =  the  bending  moment  corresponding  to  D . 

1  ,  S  1  ;  S 

The  expression  for  the  accelerations  of  the  bridge 
at  midspan  can  be  derived  directly  from  Eq .  (2-20)  and 
appears  in  the  form 


w  =  w^  y  f  (t)  Y  (r,)  (2-56) 

o  '-  n     n 


Substituting  the  second  derivative  dimensionless  function 
i|;  (n)  defined  by  Eq.  (2-45),  gives 


w  =  w^  (^)  I   i'^(T)  Yj^(n)  (2-57) 

n 

v  2 
v.'here  w  (-r-)   is  a  quantity  with  the  dimension  of  acceleration 

The  accelerations  of  the  beams  can  be  evaluated  at  the 

transverse  location  of  the  beams  in  tlie  n-direction. 


32 


2.2.4   Computer  Program 

The  method  described  in  the  preceding  sections  has  been 

(28) 
programmed      to  analyze  simple-span  bridges  having  from 

one  to  fifteen  uniformly  spaced  prismatic  beams.   The  beams 
are  assumed  to  be  arranged  symmetrically  with  respect  to 
the  longitudinal  centerline  of  the  structure.   All  interior 
beams  are  considered  to  be  identical.   The  exterior  beams 
on  either  side  of  the  bridge  are  assumed  to  be  located  along 
tlie  edge  of  the  slab  and,  while  identical  to  one  another, 
they  may  bo  different  from  the  interior  beams.   The  program 
is  limited  only  to  non-composite  sections  and  the  effect 
of  side  curbs  is  not  considered,   Tlie  vehicle  is  a  single- 
axle  load,  consisting  of  one  or  tv;o  v/heels.   The  roadv;ay 
surface  unevenness  is  represented  by  a  trigonometric 
function  in  the  longitudinal  direction. 

The  input  parameters  of  the  program  can  be  either 
dimensionloss  or  dimensional.   The  dimensionless  parameters 
have  been   defined  by  Eqs .  (2-1),  (2-24)  and  (2-43). 
The  dimensional  parameters  are  divided  into  three  major 
groups,  i.e.,  bridge  parameters,  vehicle  parameters  and 
solution  parameters.   After  the  dimensional  parameters 
are  entered,  the  frequencies  of  bridge  and  of  vehicle 
arc  computed.   Then  all  parameters  are  converted  to  a 
dimensionless  form  by  the  program  before  the  start  of  static 
and  dynamic  computations. 
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The  program  provides  results  for  the  complete  history 
of  the  response  of  the  system,  by  printing  and  plotting  out 
the  crawl  or  static  values  of  deflections  and  moments  in 
the  beam  at  midspan,  the  corresponding  dynamic  increments, 
the  dynamic  increment  of  the  wheel  reactions,  and  the 
accelerations  in  the  beams  at  midspan. 

Flow  Diagram 

A  general  flow  diagram  for  the  complete  program 
is  shown  in  Figure  2.3,  and  the  flow  diagrams  for  the 
subroutines  (QIMJ) ,  (DOMOIS) ,  (EQMON) ,  (RESLT) ,  (INTEG) 
are  shov/n  in  Figures  2.4,  2.5,  2.6,  2.7,  2.8  respectively. 

The  program  starts  with  the  reading  of  the  parameters 
specifying  the  characteristics  of  the  vehicle  and  bridge, 
and  necessary  parameters  to  carry  out  the  solutions. 
Next  the  subroutine  (QIMJ)  is  entered  to  evaluate  the 

quantities  Q   .  which  are  used  to  evaluate  the  deflections 

-m,  ] 

and  bending  moments  in  the  beams  at  midspan.   These  are 

defined  by  the  equations 

n 
o 


^  1  .        mil  r  I  ^7       /  X 

Q       .  =   —>r  sin-rr-       )        a         .Y(n.) 
m,i  2  2         '^,      mn,ini 

-"         m  n=-l  -^ 

^i    n  =     (^d^     (F^I    )    ■  (2-58) 

1,3  d    D     (E^Ij^)^ 


M.      .    =     (C    )  .  WL 
1/3  m    3 
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v.'here    D.  .  -  the  deflection  at  midspan  of  beam  i  due 
to  load  j  . 
M.  .  =  the  moment  at  midspan  corresponding  to  D .  .. 

m 

1  °  1 

(C.)  .    -   A-  y  A  Q       .  sinmTiT  (2-59) 

d    1  2         ,   ' ',  2m,-] 

■'■        V      m=l,  3.  .  .  m  -' 

A  .  O 

{C    )  .    ^   ~  J  Q       .  sinmTTT  (2-60) 

-'  o   ra=  1 ,  3  .  .  .  ■' 


Following  this,  the  subroutines  (DOMOIS)  and  (BNSCNS) 
are  entered  to  compute  and  store  the  values  of  D.    and 

1,5 

M.    ,  and  the  values  of  constant  coefficients  in  Eq.  (2-46), 
1 ,  s 

respectively . 

Before  starting  the  integration  of  the  differential 
equations,  the  second  derivatives  of  the  generalized 
coordinates  q   (i.e.,  i|)  " ,  c",  and  0"),  the  interacting 
forces,  and  the  deflections  and  moments  in  the  bridge 
arc  evaluated  for  t  =  0  by  entering  subroutines  (INTLS) 
and  (EQMON) .   The  next  part  of  the  program  is  to  integrate 
the  governing  differential  equations  numerically  and  to 
compute  the  interacting  forces,  the  deflections,  the 
moments,  and  the  accelerations  in  the  bridge  at  t  =  t  . 
The  procedure  is  done  at  each  time  interval,  t  ,  by  use 
of  subroutines  (EQMON),  (INTEG)  and  (RESLT) ,  until  t^  =  1, 
or  the  vehicle  leaves  the  bridge.   During  the  integration 
procedure,  the  maximum  results  of  the  response  are  chosen 
and  stored,  in  order  to  be  printed  out  afterward.   The 
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response  of  the  bridge  can  be  printed  out  at  the  desired 
step  of  integration  by  using  the  print-out  parameter  which 
is  read  in  at  the  beginning  of  the  program. 

The  interacting  forces,  the  deflections,  the  moments 
and  the  accelerations  of  the  bridge  at  each  step  of  time 
interval  can  be  plotted  by  using  the  subroutine  (OUT) . 

The  machine  time  required  to  obtain  a  solution 
depends  on  the  particular  problem  considered.   The  total 
time  consists  of  the  follov/ing  major  parts;  the  time 
required  to  read  in  parameters,  the  time  required  to 
calculate  the  necessary  constants  and  functions,  the  time 
required  for  the  integration  procedure  and  the  response 
computations,  and  the  time  required  for  printing  and 
plotting  out  the  results.   The  total  time  required  to  obtain 
the  response  of  60  ft.  span  bridge  with  400  integration 
steps  is  approximately  2  minutes  and  30  seconds. 

2 , 3   Analysis  of  Static  and  Dynamic  Response  of  Continuous 
Highway  Dridgos 

2.3.1  General 

In  this  analysis,  the  bridge  is  idealized  as  a  single 
continuous  beam  and  the  resulting  infinite  number  of  degrees 
of  freedom  is  replaced  by  a  discrete  system  having  a  finite 
number  of  degrees  of  freedom.   This  discretization  is 
effected  by  concentrating  the  distributed  mass  of  the  beam 
into  a  series  of  point  masses,  but  considering  the  flexibility 


^':  *.::     -'  i'U'     '"•■  '  '*  )   /I;      '''A     \\' 
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of  the  beam  to  be  distributed  in  the  actual  system.   A 

vehicle  of  the  tractor-trailer  type  is  represented  by  a 

three-axle  load  unit  consisting  of  two  interconnected  rigid 

masses.   Each  axle  is  represented  by  two  springs  in  series 

and  a  frictional  mechanism  which  simulates  the  effect  of 

friction  in  the  suspension  spring  of  the  vehicle. 

The  equations  governing  the  motion  of  the  bridge-vehicle 

system  are  formulated  in  general  terms.   They  can  be 

applied  to  continuous  bridges  of  any  number  of  spans 

as  v;ell  as  to  simple  span  bridges  or  cantilever  bridges. 

The  computer  program  has  been  developed  for  three- 

(29) 
span  bridges      having  a  uniform  cross  section  and  equal 

side  spans  and  for  a  load  unit  having  a  maximum  of  three 

axles.   This  program  has  been  modified  for  two-span 

bridges  having  ecjual  spans  for  the  same  load  unit.   The 

programs  can  handle  various  parameters  defining  the 

characteristics  of  the  bridge  and  the  vehicle,  the 

initial  conditions  of  the  vehicle,  and  the  roughness 

surface  of  the  bridge.   The  programs  provide  the  printed- 

out  and  plotted  results  including  the  interacting  forces, 

the  reactions,  the  deflections,  the  moments  and  the 

accelerations. 

2.3.2   Analysis  of  Three  Span  Continuous  Highway  Bridges 

Idealization  of  Bridge 

The  beam  used  to  represent  the  actual  bridge  is 


<  - 
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analyzed  as  a  system  having  a  finite  number  of  degrees 

of  freedom.   The  discretization  is  effected  by  replacing  the 

distributed  mass  by  a  series  of  concentrated  point 

masses  and  considering  the  beam  flexibility  to  be  distributed 

as  in  the  original  structure.   Bridge  damping  is  assumed 

to  bo  of  the  absolute  viscous  type,  and  is  approximated 

by  a  scries  of  dashpots  arranged  as  shown  in  Figure  2.9. 

The  magnitude  of  the  concentrated  mass  at  a  node  is 
taken  equal  to  the  total  distributed  mass  for  the  portion 
of  the  actual  beam  between  midpoints  of  the  panels  on  either 
side  of  that  node.   Similarly,  the  coefficient  of 
damping  for  the  associated  dashpot  is  taken  equal  to  the 
product  of  the  average  coefficient  of  damping  per  unit  of 
length  of  the  original  beam  multiplied  by  the  length  of 
the  tributary  section  for  the  node  under  consideration. 
Miisses  and  dashpots  attached  to  rigid  supports  have  no 
influence  on  the  re"'sponse  of  the  system  and  may  be  neglected. 
The  number  of  degrees  of  freedom  for  the  substitute  beam  is 
thus  equal  to  the  number  of  deflecting  masses  involved. 

The  analysis  is  based  on  the  ordinary  theory  of 
flexure  of  beams  which  assumes  the  material  is  linearly 
elastic  and  the  effects  of  shearing  deformation  are 
negligible.   Furtheirmore,  as  the  distributed  mass  of  beam 
is  replaced  by  concentrated  point  masses,  the  rotary  inertia 
effects  are  not  considered. 
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Figure  2.9   Three  Span  Continuous  Bridge  Model 


-Fiied  horizontal  planes 


Figure  2.10    Representation  of  a  Tractor-Trailer  Type 
Vehicle 
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Idealization  of  Vehicle 

Since  the  bridge  has  been  idealized  as  a  becun,  the 
width  of  the  vehicle  and  consequently,  the  rolling  effect 
cannot  be  considered  in  the  analysis.   Even  when  treated 
as  a  plane  system,  a  vehicle  is  a  very  complex  mechanical 
system.   Figure  2.10  shows  diagrammatically  the  detailed 
features  of  what  is  believed  to  be  a  complete  representation 
of  a  tractor-trailer  combination.   The  quantity  W, 

represents  the  weight  of  the  tractor  mounted  on  its  sus- 

(2) 
pension  system.   The  quantity  i,  is  the  dynamic  index 

of  the  tractor.   This  is  a  measure  of  the  rotary  moment  of 
inertia  of  the  weight  W, ,  and  it  is  defined  as  the  ratio 
of  the  radius  of  gyration  squared  to  the  product  of  the 
liorizontal  distances  between  the  two  supports  and  the 
center  of  gravity  of  the  weight.   The  dashpots  at  the  center 
of  gravity  of  W^  represent  damping  resistances  against 
vertical  motion  and  rotary  motion.   The  rigid  bar  repre- 
sents the  chassis  of  the  tractor  and  its  weight  is 
designated  as  v; .  .   The  point  masses,  with  weights  w,  and  w_, 
represent  the  mass  of  the  axles,  springs,  and  tires  for  the 
tv/o  axles.   The  quantities  W_ ,  i_,  and  W-.  refer  to  the 
trailer  and  have  the  same  meaning  as  that  of  the 
corresponding  quantities  for  the  tractor.   For  convenience 
in  presentation,  the  weights  W,  and  W„  are  referred  to  as 
"sprung"  weights  and  the  remaining  weights  are  referred  to 
as  "unsprung"  weights. 


,;'j   '  .'■  fj  J  ;' 
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The  dynamic  characteristics  of  the  tires  for  each 
axle  of  the  vehicle  are  represented  by  a  spring  and  a 
dashpot.   The  suspension  system  for  each  axle  is 
represented  by  a  massless  spring,  a  dashpot,  and  a 
frictional  device.   The  dashpot  accounts  for  the  effects 
of  shock  absorbers  or  air  suspension,  and  the  frictional 
device  accounts  for  any  frictional  force  that  may  develop 
in  the  suspension  system,  particularly  in  the  leaf  springs. 
Tlic  value  of  the  frictional  force  developed  at  any  time  is 
designated  by  F  and  the  limiting  or  maximum  possible 
value  is  designated  by  F'.   As  long  as  -F'  <F<F' 
for  a  particular  axle,  the  suspension  spring  for  that 
axle  is  inactive  (i.e.  only  the  tire  springs  deflects), 
and  the  effective  stiffness  of  that  axle  is  equal  to  the 
stiffness  of  the  tires.   On  the  other  hand,  if  F  =  +F', 
both  springs  are  active  and  the  effective  stiffness  is 
that  of  the  two  springs  acting  in  series. 

In  the  present  analysis  the  above  system  is 
further  simplified  by  (a)  neglecting  all  sources  of  viscous 
damping  and  (b)  replacing  the  "unsprung"  weights  by 
concentrated  "sprung"  weights  as  shown  in  Figure  2.11a. 
In  this  replacement  the  weight  of  the  chassis,  designated 
as  w.  in  Figure  2.10,  is  incorporated  into  the  weights  w^ 
and  w_.   This  replacement  is  justified  by  the  fact  that  the 
"unsprung"  weights  are  quite  small  in  comparison  to  the 
"sprung"  weights.   For  a  representative  tractor  the  ratio 
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of  the  total  "unsprung"  weights  to  the  "sprung"  weight 
is  about  ]./7,  and  for  a  trailer  it  is  for  all  practical 
purposes  negligible.   In  addition  to  the  three-axle  load 
unit,  in  Figure  2.11  are  shown  specialized  models  for  a  two- 
axle  and  a  single-axle  load  unit. 

With  its  velocity  specified,  the  three  axle  load  unit 
shown  in  Figure  2.11a  has  three  degrees  of  freedom.   The 
parameters  which  define  its  characteristics  are: 

(a)  the  weight  distribution  parameters  which  include 
the  weiglits  W,  ,  W„,  w,  ,  w_  and  w^,  and  the  dynamic  indices 
i^  and  i-,; 

(b)  the  geometrical  parameters  which  include  the 
axle  spacing?  1,  and  1_,  and  the  ratios  a^  through  a^. 
as  defined  in  Figure  2.11a; 

(c)  the  stiffness  parameters  for  the  tires  and  the 
suspension  springs;  for  the  i^h  axle  (i--  1,2,3),  the  stiff- 
ness of  the  tires  is  denoted  by  k^  .  and  the  stiffness 

^   t,i 

of  the  tires  and  the  suspension  springs  when  acting  in 

series  is  denoted  by  k^   . ; 
^   ts ,  1 

(d)  the  friction  parameters,  for  the  suspension  systems 
of  the  vehicle.   For  the  i^^  axle  this  is  the  limiting 
frictional  force,  F."  . 

The  Assumptions  for  the  Analysis 

The  analysis  is  based  on  the  ordinary  beam  theory, 
which  neglects  the  effects  of  shearing  deformation  and 
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Figure  2.11   Vehicle  Models 
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axial  forces.   In  addition,  since  the  distributed  mass  is 
treated  as  a  series  of  point  masses,  the  effect  of  rotary 
moment  of  inertia  does  not  enter  in  the  solution.   The 
vehicle  travels  with  constant  speed  and  is  assumed  to 
remain  in  contact  with  the  bridge  at  all  times,  and  its 
angular  displacements  are  considered  to  be  small.   It  is 
further  assumed  that  no  longitudinal  force  can  develop  at 
the  junction  of  the  tractor  and  trailer.   This  junction  is 
known  as  the  "fifth  wheel  pivot".   Finally,  all  springs 
of  the  vehicle  are  considered  to  be  elastic. 

Equations  of  Motion  for  Bridge  Model 

The  motion  of  the  vehicle-bridge  system  is  expressed 

in  terms  of  the  coordinates  z.  and  y   shown  in  Figure  2.12. 

1      -'  r  ^ 

The  coordinate  z.  denotes  the  vertical  displacement, 
measured  from  a  fixed  horizontal  plane,  of  the  point  of 
support  of  the  vehicle  mass  for  the  i^"  axle.   The 
coordinate  y   denotes  the  deflection  of  the  r^h  node  point 
of  the  beam.   This  deflection  is  measured  from  the  static 
equilibrium  position  when  the  bridge  is  subjected  to  its 
ov/n  weight  alone.   Both  coordinates  z  and  y  are  considered 
to  be  positive  when  downward. 

Figure  2.13  shows  a  three-span  continuous  bridge 
model  of  equal  side  spans  and  uniform  flexural  rigidity.  El 
The  length  of  the  center  span  is  denoted  by  L  and  the 
length  of  a  side  span  by  aL.   The  center  span  is  divided 
into  m  equal  panels  of  length  h,  and  each  side  span  is 


49 


apou  „.oj3z 


0) 

b» 

•H 
M 

n 

M-l 

o 
<\) 
>i 

Eh 

>i 
C 
< 


u 

•H 

x: 
> 

c 

QJ 
Cr> 

•H 
M 
CQ 

M-l 
O 

C 
O 
•H 

■p 
c 

•H 

e 
o 
u 


0 

u 
en 

•H 

U4 


50 


H 

apou  ^j(ui+ui) 

\ 
I 

\ 

PO 

CL 

1 
nt 

C 

ft 
\ 

3pOU^j(UJ+U) 

\ 
\ 

\ 

I 
T 

1 

1 
/ 
/ 

- 

E 

>f' 

.- 

/ 

/ 

t^ 

spou  ^,u 

IK 

O 

» 
/ 

14- 

u 

•H 

(1) 
> 

c 

(U 

'a 

m 
tn 

o 

13 
c 

•H 

■p 
c 
o 
u 

c 

a 

0) 
0) 

x; 


9pOU  n.0J9Z 


•rH    rH 

o  o 

u  s 


•H 


51 


divided  into  n  equal  panels  of  length  —ah.   The  nodes  are 
numbered  consecutively  starting  with  zero  at  the  left 
abutment  and  terminating  with  (2n+m)  at  the  right  abutment. 
The  panel  between  nodes  r  and  r  -  1  is  designated  as 
the  r^h  panel.   As  before,  the  mass  is  considered  to  be 
concentrated  at  the  node  points. 

Let  P .  be  the  interacting  force  between  the  bridge 
surface  and  the  i    axle  of  the  vehicle.   Then  the  equation 
of  motion  for  the  concentrated  mass  at  the  r    node  of  the 
beam,  m  ,  may  be  expressed  as  follows: 


2n+m 

(2-61) 

y 


^:n+m   .  -s  . 

ly^  +  cy   -   y   R^y.-  y  QP. 

r-^r     r-'r    .^,       r  -^  3  .  '•.  ^r  i 

1  =  1      ■'  1  =  1 


where  y   is  the  deflection  of  the  r^"  node,  as  previously 
r 

defined,  and  a  dot  superscript  denotes  one  differentiation 

with  respect  to  time;  c   is  the  damping  coefficient  for  the 

damper  at  the  r^"  node.   The  quantity  R -^  is  defined  as 

the  reaction-deflection  coefficient  and  represents  the  static 

reaction  at  the  r'^-^'^  node  point  induced  by  a  unit  deflection 

of  the  j th  node  point,  when  all  other  nodes  are  supported 

against  deflection.   A  reaction  is  considered  as  positive 

when  directed  upward.   In  an  analogous  manner,  Q    is  defined 

as  the  reaction-load  coefficient  and  represents  the  reaction 

at  the  rth  node  point  induced  by  a  concentrated  unit  load 

at  the  point  of  application  of  P.  when  all  nodes  are  supported 
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against  deflection.   Obviously,  when  the  unit  load  is  off 

the  bridge,  Q^=    0.      It  should  be  noted  that  the  interacting 

forces  P.  are  not  known  at  this  stage.   By  application  of 

Eq.  (2-61)  to  each  mass,  one  obtains  as  many  equations 

as  there  are  degrees  of  freedom  for  the  bridge  model. 

The  reaction-deflection  coefficients,  R -'  are 

r 

constants  for  a  given  bridge  model  and  can  be  evaluated 

by  a  modified  moment  distribution  procedure  introduced  by 

•  (37) 
T.  Y.  Lm     .   The  procedure  makes  use  of  the  concept 

of  the  effective  stiffness  and  effective  carry-over  factors 

which  are  defined  as  follows:   Consider  a  bar  ab  resting 

on  non-deflecting  supports  and  elastically  restrained  against 

rotation  at  end  b  by  a  coil  spring  having  a  stiffness  R. 

The  moment  at  end  a  necessary  to  produce  a  unit  rotation 

at  that  end  is  defined  as  the  effective  stiffness  of  that 

end  of  the  bar.   Denoted  by  K  '  ,  this  stiffness  is  given  by 

a 

the  equation, 

k   .  k,   K, 
K'=  [  1  -  _2.£-H_2-£j_£]  K  (2-62) 

^  K^  +  R      ^ 

b 

v;here  K   and  K,  are  the  Hardy-Cross  stiffnesses  of  the  bar 
a      b  -^ 

for  the  ends  a  and  b  respectively.   Similarly  k   ,  and 

k,    are  the  Hardy-Cross  carry-over  factors  from  ends  a 
D ,  a  — 

to  b  and  from  b  to  a,  respectively.   The  ratio  of  the 

moment  produced  at  end  b  to  the  applied  moment  at  a  is 

I 

defined  as  the  effective  carry-over  factor,  k   ,  ,  and 

is  given  by  the  equation: 
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k  '  ,     =   5- k^    ,  (2-63) 

^'^  (1-k      ,k,        )K.    +R      ^''° 

a,b   b,a      b 


For    a    prismatic    bar,    K      =    K,     =    K,    k      ^    ~   ^u        ~    ~  1/2  f 

a      D  SL  f  D  U  f  3. 


and  the  above  equations  become 


K'  =  [1  -i  — L--]  K  (2-64) 

K  +  R 


and 


^a  b  =  —T-^ ^  ^2-65) 

For  a  continuous  beam  the  coil  spring  symbolizes  the  continui- 
ty of  a  particular  span  with  the  adjacent  spans. 

In  the  course  of  calculating  the  coefficients  R  -^ 

r 

by  this  procedure,  one  calculates  also  the  moments  at  the 
nodes  due  to  a  unit  displacement  at  the  j t"  node.   These 
moments  arc  termed  as  moment-deflection  coefficients  and 
are  designated  by  J  -'  .   In  evaluating  the  coefficients 
R -^  and  J.  ,  the  following  quantities  are  used.   In  all 
cases,  it  is  assumed  that  the  bridge  model  is  supported 
against  deflection  at  the  node  points. 

(a)  Effective  Stiffness  Coefficients.   Consider 
the  portion  of  the  bridge  model  between  the  left  hand  abut- 
ment and  the  r^^  node  as  a  beam  continuous  over  non-deflec- 
tive supports  at  the  nodes.   Then  the  effective  stiffness 
of  the  beam  at  end  r  may  be  stated  as  the  product  of  a 


■..<JhJ.: 
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dimensionless  stif f ness~coef f icient  C   and  the  quantity 
4EI/h,  v;here  h  refers  to  the  length  of  a  panel  in  the  center 
span  of  the  bridge  model.   By  application  of  Eq .  (2-64) 
it  can  be  shown  that  the  coefficient  C  .  is  given  by  the 


following  recurrence  forrr^ula: 


1  - 


I  ~ 


+    C 


r-1 


(2-66) 


th 


wliere  h   is  the  length  of  the  r'-''  panel.   For  a  panel  on 

the  center  span,  h  =h;  and  for  a  panel  on  a  side  span, 

h  =  —  ah . 
r   n 

It  should  be  noted  that,  because  of  symmetry,  the 
dimensionless  coefficient  for  the  stiffness  at  node  r 
for  the  portion  of  the  beam  between  the  r  "  node  and  the 


right  hand  support  is  equal  to  C 


2n+m-r ' 


(b)  Effective  Distribution  Factors.   The  effective 
distribution  factor  for  the  right  hand  side  of  the  r 
node,  designated  as  d.,  is  given  by  the  expression, 


■2n+m-r 


r    C   +  C„  ^ 

r     2n+m-r 


(2-67) 


The  distribution  factor  for  the  left  hand  side  of  the  r 

node  is  1  -  d  .        '^ 
r 

(c)  Effective  Carry-Over  Factors.   The  effective 
carry-over  factor  from  node  r  to  node  r  -  1  is  designated 


th 
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ar  k     , .   By  npplication  of  Eq .  (2-65) ,  one  finds  that 
r ,  r-1     -ILL 

k  '    ,  -  .,  ,     (2-68) 

2   h       r-1 

r 

Gincn  the  beam  is  syirunetrical  about  the  center  line,  it 
follows  that 


k       =  k 
r,r+l     2n+m-r ,  2n+in-r-l 


For  the  sake  of  brevity,  in  the  following  discussion  the 

quantity  k  '    -.is    designated  as  k  '  . 
'       -*   r,r-l         ^  r 

To  determine  the  moment-deflection  coefficients  J  -' 

and  the  reaction-deflection  coefficients  R -^  ,  the  j  "  node 

r 

of  Lhf^  model  is  first  displaced  by  a  unit  amount,  and  by 

kccpinq  all  nodes  fixed  against  rotation  the  fixed-end 

moments  produced  at  the  nodes  (j-1),  j  and  (j+1)  are 

evaluated.   The  resulting  unbalanced  moments  (if  h.  =h.  ,, 

th.ere  is  no  unbalanced  moment  at  the  j^^  node)  are  then 

distributed  and  carried  over  by  use  of  the  quantities  given 

in  Eqs .  (2-67)  and  (2-68).   The  final  moments  at  the 

nodes  yield  the  coefficients  J  ^  .      The  reaction-deflection 
-'  r 

coefficients  R  -^  are  next  evaluated  from  the  equation 


.    J   ,3  -  J  ^      J    ^    -   J    J 

p3  -   r-1     r       r r+1  ,^  ^p» 

p.,         R-  ^,— (2-69) 

The  quantities  C   and  d   are  used  only  to  evaluate 

1  r      r             -^ 

the  coefficients  R  ■'  and  J  -'  ,  whereas  the  carry-over  factors 

r  r  '                 ^ 
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k'  and  the  quantities  R -'  and  J-'  are  used  repeatedly 
in  later  stages  of  the  solution. 

Another  quantity  needed  in  subsequent  computation  is 
the  total  angle  change  produced  at  the  r^h  node  when  the  beam 
is  cut  at  the  r^"  node  and  a  unit  bending  moment  is  applied 
on  the  two  sides  of  that  node.   As  before,  all  nodes  are 
assumed  to  be  held  against  deflection.   This  angle  change 
is  denoted  by  0_  and  is  given  by  the  expression, 

^  r    ^C     C„  ^    ■'  4EI 
r     2n+m-r 

By  use  of  Eq .  (2-68),  the  above  expression  may  be 
written  as 

It  should  be  emphasized  that  the  quantities  defined 
above  depend  only  on  the  characteristics  of  the  bridge 
model. 


Equations  of  Motion  for  a  Vehicle 

Let  P  ,  .  be  the  reaction  at  the  i'^^  axle  when  the 

St,  1 

vehicle  is  in  a  position  of  static  equilibrium.   With  P. 
denoting  the  dynamic  reaction  at  any  time  t,  the  disturbing 
force  for  the  i^^  axle  is  P .  -  P  .  .  and  the  equation  of 

1     St, 1  ^ 

motion  for  a  three-axle  vehicle  can  be  stated  in  the  form: 


57 


^11   ''^12   ^13 


^12   ^22   ^23 
^13   ^23   ^33 


1    st ,  1 

^2  ~  ^st,2 

3    s  t ,  3 


(2-71) 


v/here  g  is  the  gravitational  acceleration,  W  is  the  total 
weight  of  the  vehicle,  and  a^,  through  a^,  are  dimensionless 
coefficients  given  by  the  following  expressions: 


^11 


^  W  W    w 

^^1  ^^1^2^1^  -r+  ^5^^^3^-'^3^4i2^  if-'lf 


12 


W  W 


13 


^22 


a3a^a5(l-i2)  -^ 


W, 


,,  W    w 

(32  +a^a2i^)-^^+  (1  -  a^)  U^     +333^12)-^+  — 


(2-72) 


^23 


^3^4(1  -a^)  (1  -  12)  -^ 


'33 


I       2    ■    .  ,  ^2  ,  ^3 


The  details  of  derivation  are  presented  in  Ref.  (29). 

Premultiplication  of  Eq.  (2-71)  by  the  inverse  of 
matrix  A  yields. 
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'^U   ^12   ^13 
^12   ^22   ^23 


L^13 


'^23   ^33 


^1  "  ^st,l 
^2  ~  ^'st,2 
-3    s  t ,  3. 


(2-73) 


Let  u.  be  the  shortening  of  the  suspension-tire 
system  of  the  i^^  axle.  The  relationship  between  the 
interacting  force.  P.,  and  the  shortening,  u.,  of  the 
combined  suspension-tire  system  is  shown  in  Figure  2.14. 
Also  included  in  this  figure  is  a  diagram  showing  the 
relationship  between  u.  and  the  frictional  force,  F.. 

As  the  load,  P,  is  increased  above  its  initial  value, 
P  . ,  the  deformation  of  the  tire  spring  increases  linearly, 
but  the  suspension  spring  does  not  deform  at  first,  the 
initial  increase  in  the  load  of  the  suspension  system  being 
resisted  entirely  by  friction.   The  initial  paths,  oa ,  of  the 
P-u   diagram  and  the  F-u  diagram  in  Figure  2.14  are, 
therefore,  parallel,  and  the  slope  of  these  paths  is  equal 
to  the  stiffness  of  the  tire  spring,  k  .   This  relationship 
continues  until  the  frictional  force  attains  its  limiting 
value,  F'.   At  that  stage,  the  suspension  spring  engages, 
and  the  effective  stiffness  of  the  tire-suspension  system 
becomes  equal  to  the  stiffness,  k   ,  of  the  two  springs 
acting  in  series,  in  which 


1 
ts 


J_4-  JL 
•^t    ^S 


(2-74) 
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(a)      P-u    Diagram 


(b)      F-u    Diagram 


Figure  2.14    Force  Deformation  Relationships  for 
Tire  Suspension  System  for  Vehicle 
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and  k   is  the  stiffness  of  the  suspension  spring.   It 

follows  that,  along  the  segments  ab,  the  slope  of  the  P-u 

diagram  is  k  ,  the  slope  of  the  P-u  diagram  is  k   ,  and  the 
s  t  s 

frictional  force  has  its  limiting  value,  F'.   Now,  if 
at  the  load  corresponding  to  points  b  of  the  diagrams 
the  direction  of  loading  were  reversed,  the  tire  spring 
would  rebound  elastically,  but  the  deformation  of  the  sus- 
pension spring  would  not  change,  as  the  decrease  in  load 
would  be  resisted  entirely  by  friction.   The  force  in  the 
frictional  device  would  decrease  at  the  same  rate  as  the 
load,  and  the  unloading  paths,  be,  on  both  the  P-u  diagram 
and  the  F-u  diagram  would  be  parallel  to  each  other  and  to 
the  initial  loading  path,  oa .   This  condition  would 
continue  until  the  load  is  reduced  by  2F',  at  which  time 
the  frictional  force  would  attain  its  limiting  value  in  the 
negative  direction,  and  the  suspension  spring  would  once 
again  engage.   A  possible  path  on  these  diagrams  beyond  this 
point  is  represented  by  the  segments  ed  -  de  -  ef  -  fg. 

This  behavior  in  in  good  agreement  with  available  test 
data(25)(26)_ 

To  summarize,  the  parameters  that  define  the  behavior 
of  a  tire-suspension  system  are  the  stiffness  of  the  two 
springs,  k.  and  k^,  and  the  limiting  value  of  the  frictional 
force,  F'.   The  latter  force  may  conveniently  be  expressed 
in  the  form 
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liP 


St 


(2-75) 


in  which  p  =  a  dimensio riess  coefficient,  referred  to  as  the 
"coefficient  of  interleaf  friction". 


Evaluation  of  Interacting  Forces 

Equations  (2-61)  and  (2-71)  are  coupled  through  the 
interacting  forces  P.,  which  remain  to  be  evaluated. 
Let  time,  t,  be  measured  from  the  instant  the  first  axle 
enters  the  bridge.   Then  the  interacting  force  at  time 
t  is  given  by  the  equation 

rt    du. 


P.  =  P. 


t=0 


k.  -3 —  dx 


(2-76) 


where  P. 


.1    is  the  initial  value  of  P . ,  k-  is  the  instantane- 
ous  effective  stiffness  of  the  suspension-tire  system  for 


the  i^h  axle  at  any  time,  t,  and  u.  is  the  corresponding 
shortening  of  the  suspension-tire  system.   If  at  the 
instant  it  enters  the  bridge,  the  vehicle  is  at  the  posi- 
tion of  static  equilibrium,  the  initial  force  Pi]    =  P  .  • 
The  instantaneous  stiffness  k.  depends  on  the 
history  of  the  shortening  u.  as  shown  in  Figure  2.14. 
The  shortening  u.  can  be  expressed  in  the  form 


u .  =   z  .    +  d_, .  -  y„  ■ 
1    1    Pi   -'Pi 


(2-77) 


where  z.  is  the  coordinate  for  the  i^"  axle,  as  previously 


lJ  • 
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defined,  and  dp.  is  the  deviation  of  the  bridge  profile 

at  the  point  of  application  P..   This  deviation  is  measured 

from  a  horizontal  line  passing  through  the  first  abutment 

of  the  bridge,  as  shown  in  Figure  2.12.   The  quantity  yp . 

represents  the  deflection  of  the  bridge  at  the  point  where 

p.  acts.      . 
1 

The  term  "bridge  deviation"  includes  the  dead  load 
deflection,  initial  camber,  grade,  vertical  curve  and  road- 
way unevenness.   For  a  given  bridge,  these  quantities 
may  be  presumed  to  be  known.   This  deviation  is  considered 
to  be  positive  when  upward.   The  deflection  y  .  is  measured 
from  the  static  equilibrium  position  of  the  beam,  when 
acted  upon  by  its  own  weight,  and  it  is  positive  when 
downward. 


Computation  of  the  Response 

Application  of  Eqs.  (2-61)  and  (2-73)  to  each  con- 
centrated mass  of  the  bridge  model  and  to  each  axle  of 
the  vehicle  yields  a  set  of  simultaneous,  second  order 
differential  equations,  equal  in  number  to  the  number  of 
degrees  of  freedom  of  the  bridge-vehicle  system.   These 
e(]uations  can  be  solved  by  a  numerical  integration  method 
in  which  the  evaluation  of  the  interacting  forces  P- 
is  a  major  intermediate  step. 

As  the  integration  of  the  differential  equations  is 
carried  out,  the  values  of  all  the  coordinates,  of  all 
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accelerationG  at  nodes,  and  of  the  interacting  forces  are 
determined.   From  these  quantities  the  values  of  the 
corresponding  deflections,  moments  and  reactions  at  any 
desired  section  may  then  be  evaluated  by  statics. 

Integration  of  Equations  of  Motion 

The  differential  equations  of  motion  defined  by 
Eqs.   (2-61)  and  (2-73)  have  been  integrated  numerically 
using  an  iterative  procedure  within  each  integration 
step.   To  describe  the  method,  assume  at  some  time,  t  , 
the  values  of  all  displacements  and  of  their  first  two 
derivatives  are  known,  and  it  is  desired  to  determine  the 
correspoiiding  quantities  at  time: 

t  , ,  ^  t   +  At  (2-7  8) 

s  +  i     s 

in  which  At  is  a  small  time  interval.   In  the  method  used, 

an  assumption  must  first  be  made  concerning  the  manner  in 

\v'hich  the  accelerations  vary  within  the  interval,  At. 

This  variation  is  assumed  to  be  linear  in  the  present  study. 

The  velocities  and  displacements  at  t  , ,  can  then  be 

^  s  +  1 

determined  in  terms  of  the  known  quantities  at  t   and  the 

^  s 

unknown  acceleration  at  t  , ,  by  use  of  the  following  recur- 

s  +  1 

rence  formulas: 


>r,s  +  l  =  ^^r,s  '   i^^^r,s  +  i^^^r,s.l  ^2-79) 


^r,s  +  l  =  ^^r,s  -^  ^t^r,s  "^  I  ^^^^^  '  ^'r ,  s  ^  ^^^^^  '"^r  ,s-^l 

(2-80) 


r.  ,  •<  ■ 
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in  v;hich  q   is  either  the  displacement  of  a  node,  y.,  or  the 

displacement  of  a  point  of  support  of  the  vehicle,  z.;  a 

dot  subscript  denotes  differentiation  with  respect  to  time; 

and  the  subscripts  s  and  s+1  following  a  comma  identify 

quantities  corresponding  to  times  t   and  t   ,  respectively. 

Next,  the  unknown  accelerations  at  t  , ,  are  evaluated  by 

s  +  i 

satisfying  the  equations  of  motion  at  that  instant.   The 
details  of  the  iterative  procedure  used  are  as  follows: 

1.  Define  the  longitudinal  position  of  each  axle  of  the 

vehicle  at  time  t  , , . 
s  +  1 

2.  Assume  the  accelerations  y.  and  z.  at  t  ,,  are 

^2  1      s+1 

the  same  as  those  at  t  ,  and,  by  application  of  Eqs .  (2-79) 

and  (2-80)  determine  approximate  values  for  the  velocities, 

y.   ,,  and  i.        ,,  and  for  the  displacements,  y.   ,,  and 
-'3,  s  +  1       i,s+l  ^  ■^j,s+l 

1,S  +  1 

3.  Determine  improved  values  for  y.  ,t>Y-       .1 

3  I s  +  1    3 , s  +  ± 

and  y.   , i  by  proceeding  as  follows: 
-'  J  ,  s  )  1 

a)  By  application  of  Eq .  (2-61)  to  the  first 

node,  r=  1,  obtain  an  improved  estimate  for  y 

l,s  +  l 

The  major  operation  in  this  step  is  the  computation  of  the 

quantities  yp.-^y.  and  ]  Q^T?  .  .      The  values  of  P.  used  in  this 
■'  ^l-']      jli  1 

step  are  those  corresponding  to  the  beginning  of  the  time 

interval,  i.e.,  P.   ,  and  the  values  of  y.  and  y.  are  the 
1  f  s  3  3 

approximate  values  determined  in  Step  2. 

b)  From  Eqs.  (2-79)  and  (2-00)  calculate  the  values 

of  y,   ^  and  y,   , ,  corresponding  to  the  acceleration 
1 I s+1       1 , s+1 
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y,   .T  determined  in  Step  3a. 
^  1, s+1 

c)  Repeat  Steps  3a  and  3b  for  the  remaining  nodes 
(r = 2, 3 , . , 2n+m) ,  considering  successively  one  node  at  a 
time.   In  each  computation,  use  the  latest  available  values 
of  y . ,  z.  and  their  derivatives. 

4.  Determine  improved  values  for  z.   ,,,  z.   ,,  and 

"^  i,s  +  l    i,s  +  l 

z,   ,,  proceeding  as  follows: 
i,s  +  l  t-        ^ 

a)  Compute  a  first  approximation  for  P^    , ,  the 
value  of  the  first  interacting  force  at  the  end  of  the 
time  interval.   The  various  steps  involved  in  this  compu- 
tation are  described  in  detail  in  the  following  section. 

b)  Evaluate  z*     .  from  Eq.  (2-73),  using  the  latest 

available  values  of  P . .   The  value  of  P,  used  is  that 

1  1 

determined  in  Step  4a,  and  the  values  of  P„  and  P-, 
are  those  corresponding  to  the  beginning  of  the  time  inter- 
val. 

c)  From  Eqs .  (2-79)  and  (2-80)  determine  the  values 

of  z,    ,  and  z,   ,,  corresponding  to  the  acceleration 
l,s+l       l,s+l 

computed  in  Step  4b. 

d)  Ref)eat  Steps  4a  and  4c  for  the  remaining  axles, 

considering  one  axle  at  a  time  and  always  using  the  latest 

available  values  of  P.  and  z. . 

1       1 

5.  Compare  the  derived  accelerations  with  the  pre- 
viously available  values.   If  the  difference  between 
the  initial  and  derived  values  for  any  one  coordinate 
exceeds  a  prescribed  tolerance,  repeat  Steps  3  through  5, 
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always  using  the  latest  available  values  of  P . ,  z.,y.  and 

their  derivatives.   When  all  differences  are  less  than  the 

prescribed  tolerance,  the  integration  for  the  time  interval 

under  consideration  is  completed.   One  then  proceeds  to  the 

next  time  interval  and  repeats  the  process. 

In  the  computation  of  P.    ,  needed  in  Steps  4a  and  4d, 
'^  i,s+l  ^ 

it  is  assumed  that  the  effective  stiffness  of  the  suspension- 
tire  system  remains  constant  v;ithin  a  time  interval  of 
integration.   Under  this  assumption,  Eq.  (2-76)  may  be  re- 
written in  the  form 


P.       ^,    =   P.         +    (Au.  )k. 

1,S+1  1,S  11 


P.       .-,=?.         +    (u.       ^,  -  u.       )k.  (2-81) 

i,s+l  i,s  1,S+1  1,S       1 


The  value  of  u.   ,,  is  determined  by  Eq .  (2-77).   The 
i,s+l  -" 

value  of  y  .  in  this  equation  is  evaluated  by  superimposing 
the  following  tliree  components:   (i)  deflection  due  to  the 
moments  acting  at  the  two  ends  of  the  panel;  (ii)  deflection 
due  to  the  force  or  forces  P.  acting  on  the  panel;  and 
(iii)  deflection  due  to  a  rigid  body  displacement  of  the 
panel.   The  deflection  y  .  must  be  evaluated  for  each  cycle 
of  iteration  in  the  integration  process,  since  the  values 
of  y.  and  P.  vary  from  one  cycle  to  the  next. 

The  value  of  k.  in  Eq.  (2-81)  is  determined  by  making 
use  of  the  F-u  diagram  for  that  axle,  as  shown  in  Figure  2.14. 
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Let  the  frictional  force  corresponding  to  u .    be  denoted  by 

1  ,  s 

F.   .   In  the  F-u  diagram,  imagine  a  straight  line  which 
1 ,  s 

passes  through  the  point  (u.   ,  F.   )  and  is  parallel  to  the 

initial  line  oa.   Let  u.    be  the  abscissa  of  the  point 

i,s 

of  intersection  of  this  inclined  line  and  a  horizontal  line 

corresponding  to  the  positive  value  of  F'.   Similarly,  let 

Q 
u.  ^  represent  the  point  of  intersection  of  this  inclined 
1  ,  s 

line  with  a  horizontal  line  corresponding  to  the  negative 

value  of  r' .   Then  the  value  of  k.  is  determined  from  the 

1 

following  criteria: 


easel  u.  +Au.lu.  k.=k.  . 

^.Q     i,s       1    i,s  1    t,i 

case  2  u.  +  Au.  >  uV  k .  =  k . 

i,s      1    i,s      1    ts,i 

case  3  u .  +  Au .  >  u .  k .  =  k .  . 

^_  ^  Q     i,s       1    i,s       1    t,i 

case  4  u .  -f  Au.  lu.  k.  =k^ 

1,S      1    i,s      1    ts,i 


It  follov/s  that  the  selection  of  k.  depends  only  on  the  value 

of  Au,  u   and  u  .   The  value  of  F  need  not  be  computed. 

u      / 
For  cases  1  and  3  the  values  of  u   and  u   at  time  t  ., 

s+i 

are  the  same  as  those  at  time  t  ,  whereas  for  cases  2  and  4 
they  differ  by  the  amount  Au. 

The  time  interval  or  integration  employed  should  be 
sufficiently  small  so  that  successive  cycles  of  iteration 

converge  and  the  resulting  solution  is  stable  and  accurate. 

(35) 
For  a  linear  system,  it  has  been  shown   '   the  method  used 

is  convergent  if  the  integration  step 
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At  1  0.389  T  (2-82) 

in  which  T  is  the  shortest  natural  period  of  vibration  of  the 

system.   The  total  time  between  the  instant  the  front 

axle  enters  the  bridge  and  the  instant  the  last  axle  leaves 

the  bridge  is  (1  4  2a  + S,  + S^ ) L/v.   Let  N  be  the  number 

of  steps  used  for  a  complete  solution,  then 

(1  +  2a  +  S,  +  S^)  L 
^^     0.3^9  VT (2-83) 

2.3.3  Analysis  of  Two-Span  Continuous  Highway  Bridges 
The  method  of  analysis  is  the  same  as  that  used  to 

solve  the  three  span  continuous  bridges.   The  only  difference 
between  the  two  methods  is  the  bridge  model  used  in  the 
analysis.   The  bridge  model  in  this  case  is  a  two  span 
continuous  beam  of  equal  spans  and  uniform  flexural  rigidity, 
EI,  as  shown  in  Figure  2.15.   The  length  of  either  span  is 
denoted  by  L.   Both  spans  are  divided  into  n  equal  panels 
of  length  h  and  tfie  total  number  of  nodes  is  2n.   As  before, 
the  mass  is  considered  to  be  concentrated  at  the  node  points. 
The  total  time  for  the  vehicle  to  cross  the  bridge  is 
(2  +  Sj^  +  S2)L/v. 

2.3.4  Computer  Programs 

The  method  described  in  the  preceding  articles  has 

(29) 
been  programmed      for  digital  computation  of  the  dynamic 

response  of  uniform  continuous  bridges.   Two  different 
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programs  have  been  prepared.   The  first  program  was 
developed  to  compute  the  static  and  dynamic  histories  of 
uniform  three-span  continuous  bridges  with  equal  side  spans 
when  traversed  by  a  single  vehicle  load  having  either  one, 
two  or  three  axles,  while  the  other  has  been  modified  to 
compute  the  dynamic  response  of  uniform  two-span  continuous 
bridges  with  equal  span  lengths  when  traversed  by  the  same 
characteristics  of  a  vehicle  load. 

Both  programs  can  handle  either  the  dimensional  or  non- 
dimensional  input  and  provide  the  output  of  static  and  dynamic 
histories  including  the  reactions  at  the  supports,  the 
moments  over  the  interior  piers,  the  moments  and  deflections 
at  selected  points  on  the  spans,  the  interacting  forces  between 
the  axles  of  the  vehicle  and  the  bridge,  and  the  accelerations 
of  the  bridge  at  node  points.   The  programs  also  provide  the 
maximum  values  of  these  results  and  can  plot  out  the  static 
and  dynamic  histories  by  using  a  suitable  plotting  device. 

Description  of  Three-Span  Program 

A  general  flow  diagram  for  the  complete  program  is 
shown  in  Figure  2.16.   The  program  can  be  divided  into  three 
major  parts.   The  function  of  each  part  and  the  sequence 
of  operations  involved  are  described  in  the  following. 

(a)  Fart  I    The  program  starts  with  input  of  the 
parameters  specifying  the  characteristics  of  the  vehicle 
and  bridge,  excluding  those  parameters  which  specify  the 
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Yes 
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Figure    2.16         General    Flow   Chart    for   Complete    Program 
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initial  conditions  of  the  system.   Next  the  subroutine 

(G^ )  is  entered,  and  the  constants  including  the 

effective  carry-over  factors,  k  ',  the  moment-deflection 

coefficients,  J  ,  the  reaction-deflection  coefficients, 
r 

R -^ ,  and  the  angle  coefficients,  0  ,  are  computed  and  stored. 

The  coefficients  required  for  the  static  computation  are 

determined  by  taking  n  =  m=l.   Following  this,  subroutine 

(G_)  is  entered  to  compute  the  value  of  P  .  and  the 
^  st 

elements  of  the  matrix  B  in  Eq .  (2-73).   The  values  of  P  . 
are  determined  in  terms  of  the  parameters  specifying  the 
geometry  and  weight  distribution  of  the  vehicle.   The 
matrix  B  is  determined  by  first  forming  matrix  A  in  Eq . 

(2-71)  and  then  inverting  it. 

(b)  Part  II    This  part  performs  three  major  tasks. 
The  first  step  is  the  determination  of  the  maximum  static 
effects  using  subroutines  (SH)  and  (SS) .   Subroutine  (SH) 
determines  the  static  history  and  subroutine  (SS)  selects  the 
maximum  static  values.   Both  subroutines  make  use  of  a  number 
of  subroutines,  of  which  the  most  important  are  (a)  subroutine 

(DI'P),  which  defines  the  position  of  the  axles  on  the  bridge, 

(b)  subroutine  (SMD) ,  as  shown  in  Figure  2.17,  which  calcu- 
lates the  static  moments  and  deflections  at  any  specified 
point  on  tlie  bridge  by  use  of  subroutines  (SCM)  ,  (SDCP), 

(STMOME)  ,  (5TAREA)  ,  (SMRP)  and  (STREAC)"  . 

In  this  part  of  the  computation  the  characteristic 
coefficients  corresponding  to  n  =  m=l  are  used.   If  the 
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maximum  static  effects  are  knovm,  they  may  be  fed  into  the 
machine  at  the  beginning  of  the  problem  and  the  calculation 
of  the  maximum  static  values  is  bypassed  by  transferring 
control  directly  to  the  next  operation. 

The  second  function  of  this  part  of  the  program  is 
to  set  the  initial  conditions  of  the  pi'oblem  at  the  so-called 
"neutral  condition".   This  is  done  by  subroutine  (NIC) . 
If  the  initial  conditions  are  different  fiom  these,  the 
appropriate  parameters  are  read  in  at  this  stage. 


Th 


e  third  function  is  performed  by  subroutine  (G.) 
to  establish  for  each  axle  the  values  of  u   and  u  which 
are  consistent  with  the  initial  values  of  F. .   Three  values 
are  required  to  determine  the  value  of  the  effective  stiff- 
ness of  the  suspension-tire  system. 

(c)  Third  Part    The  principal  functions  of  the  third 
part  are  to  integrate  Eqs .  (2-61)  and  (2-73)  numerically 
and  to  compute  dynamic  deflections,  moments,  reactions  and 
accelerations.   The  major  operations  involved  are: 

(i)  To  determine  the  position  of  the  vehicle  at  the  end 
of  each  time  interval  by  use  of  subroutine  (DTP)  . 

(ii)  To  integrate  the  equations  of  motion  for  this  time 

u      / 
interval,  and  to  store  the  values  of  P,  u,  u   and  u   at  the 

end  of  this  time  interval.   This  operation  is  carried  out  by 

subroutine  (DINTE)  as  shov/n  in  Figure  2.18. 
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(iii)  To  evaluate  dynamic  deflections,  moinents  and  reac- 
tions by  use  of  subroutine  (DMD)  as  shown  in  Figure  2.19. 
Subroutine  (DMD)  makes  use  of  subroutines  (DMDM) ,  (DMDP) , 
(DDCM) ,  (DDCP),   (DYMOME) ,   (DYNREA) ,  (DYREAC),   (DMRP), 
(DMDIN)  and  (DDEFY) .   Since  these  quantities  are  computed 
at  the  end  of  each  time  interval  and  are  compared  with 
the  previous  values,  the  maximum  values  can  then  be  ob- 
tained.  The  foregoing  steps  are  repeated  until  the  last 
interval  is  reached. 

Description  of  the  Two-Span  Program 

The  three-span  program  was  modified  to  be  used  to  com- 
pute the  dynamic  response  of  two  span  continuous  bridges  by 
changing  some  subroutines  in  the  following  way: 

a)  The  stiffness  properties  of  the  two  span  bridge 
are  obtained  by  using  subroutine  (G, )  of  the  three-span 
program  using  as  input  three  equal  spans,  the  third 
of  which  has  a  moment  of  inertia  of  zero.   The  modified 
moment  distribution  is  used  and  the  coefficients  including 
the  effective  carry-over  factors,  k  ',  the  moment-deflection 
coefficients,  J  -^ ,  the  reaction-deflection  coefficients, 
R -^ ,  and  the  angle  coefficients,  0  ,  are  computed  and  stored 
for  use  in  the  rest  of  the  program.   Those  coefficients 
calculated  for  the  second  span  are  moved  into  the  third 
span  position,  therefore  the  second  span  is  skipped. 
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b)  All  subprograms  are  altered  so  the  truck  moves  off 
of  the  first  span  to  the  third  span  and  the  stiffness  proper- 
ties calculated  as  above  are  for  the  first  span  joined 
directly  to  the  third  span.   By  creating  the  two-span 
program  in  this  way,  the  bookkeeping  required  for  getting 
a  truck  on  or  off  the  bridge  does  not  need  to  be  altered. 


.   ■•  ■  5x.i^.    ".as 
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CHAPTER  III 
ACCELERATION  STUDIES  OF  SIMPLE  SPAN  HIGHWAY  BRIDGES 

3 . 1   General 

The  purpose  of  this  chapter  is  to  study  the  major 
parameters  which  affect  the  accelerations  of  the  simple- 
span  highway  bridges  under  a  moving  single-axle  vehicle. 
The  study  is  based  on  the  method  of  analysis  and  computer 
program  described  in  Art.  2.2.   The  sources  of  information 
used  in  this  chapter  include  the  Standard  Plans  of  Highway 

Bridge  Supnr;;tiuctures  of  the  U.S.  Bureau  of  Public  Roads 

(41 ) 
{h(.>reafter  abbreviated  as  B.P.R.)     .   The  study  is  re- 
stricted to  bridges  of  I-beam  type  composed  of  steel 
girders  and  a  reinforced  concrete  deck  with  non-composite 
action,  and  the  side  curbs  of  the  bridge  are  not  taken 
into  account.   All  frequencies  of  the  bridges  in  this 
chapter  are  natural  fundamental  frequencies  and  no 
damping  is  considered  for  either  the  vehicle  or  the 
bridge.   Throughout  this  chapter,  the  acceleration  under 
discussion  is  at  midspan. 

The  major  parameters  that  affect  the  acceleration 
of  the  bridge  can  be  classified  into  the  following 
four  groups: 
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a)  Solution  parameters  to  assure  the  convergence 
and  stability  of  the  solutions. 

b)  Bridge  parameters  including  the  span  length,  the 
width  and  the  flexibility  of  the  girder. 

c)  Vehicle  parameters  including  the  number  of  wheels, 
the  velocity  and  the  transverse  position  of  wheels 
on  the  bridge. 

d)  Parameters  expressing  the  initial  conditions  of 
the  bridge  including  the  surface  roughness. 

3  .  2   Solution  Parameters 

In  the  analysis  of  simple  span  bridges  described  in 
previous  chapter,  the  following  parameters  are  needed 
to  assure  the  convergence  £ind  stability  of  the  solutions: 

a)  The  number  of  integration  steps  (N) . 

b)  The  maximum  number  of  terms  used  in  the  longitudinal 
direction  in  the  computation  of  the  static  effects 

c)  The  maximum  number  of  Y   functions  used  in  static 

n 

computation  (n  ) . 
o 

d)  The  maximum  number  of  Y   functions  used  in  the 

n 

equations  of  motion  (n, ) . 

In  a  previous  study  by  Oran     ,  it  was  shown  that 

convergence  of  the  static  results  v;as  satisfied  by  using 

fifteen  terms  in  the  trigonometric  series  expression 

(m  =15)  and  eight  Y   functions  (n  =8).   Table  3.1 
o  ^     n  o 


,  (:  -1  ■•     :?../- 
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Table  3.1    Comparison  of  Accelerations  Obtained  by  Using 
Different  Numbers  of  Integration  Step 


Acceli 

erations  at 

Midspan 

x/L 

N 

Edge  Beam 

Beam  No,  2 

Center  Beam 

200 

-0.255277 

-0.297954 

-0.315632 

300 

-0.253380 

-0.298322 

-0.316935 

0.20 

400 

-0.252714 

-0.298451 

-0.317395 

500 

-0.252405 

-0.298511 

-0.317608 

600 

-0.252237 

-0.298543 

-0.317723 

200 

-0.159955 

-0.410796 

-0.514699 

300 

-0.159984 

-0.411008 

-0.514935 

0.50 

400 

-0.160029 

-0.411076 

-0.515064 

500 

-0.160055 

-0.411107 

-0.515096 

600 

-0.160072 

-0.411124 

-0.515113 

200 

-0.258546 

-0.331062 

-0.361099 

300 

-0.266470 

-0.329931 

-0.356219 

0.80 

400 

-0.269257 

-0.329532 

-0.354500 

500 

-0.270549 

-0.329347 

-0.353702 

600 

-0.271251 

-0.329246 

-0.353269 

200 

0.103452 

0.159115 

0.182171 

300 

0.112092 

0.156269 

0.174567 

1.00 

400 

0.115142 

0.155268 

0.171888 

500 

0.116558 

0.154804 

0.170645 

600 

0.117327 

0.154551 

0.169969 
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shows  the  solutions  for  the  accelerations  at  midspan 
considering  different  numbers  of  integration  steps.   The 
accelerations  are  obtained  at  the  different  vehicle  po- 
sitions on  the  bridge.   The  vehicle  is  a  single-axle  loading 
consisting  of  two  wheels  and  moves  in  the  longitudinal 
direction  of  the  bridge.   The  vehicle  enters  and  leaves 
the  bridge  when  the  values  of  x/L  are  equal  to  zero  and 
one,  respectively.   The  characteristics  of  the  bridge 
and  vehicle  are  defined  by  the  following  non-dimensional 
parameters:   c  =  0.4,  X=12.5,  k  =  0,  y^O.OS,  p  =  0.5, 
w  =  0,  (1^0.15,  p  =  0.2,  K  =  0.25,  n-,=0.35,  n2  =  0.65, 
m  =15,  n  =8,  n,  =4,  n- =  4  and  the  bridge  model  is  a  five 
girder  bridge.   Both  the  bridge  and  the  vehicle  are  assumed 
to  be  initially  in  their  static  positions  of  the  equilibrium, 
and  the  bridge  has  a  level  surface.   It  can  be  seen  from 
Table  3.1  that  the  accelerations  converge  with  the  increas- 
ing values  of  N.   Since  the  computer  cost  is  higher  when 
N  is  larger,  and  there  are  small  differences  of  the  accelera- 
tion values  between  400  and  600  steps  of  integration, 
400  integration  steps  were  used. 

Table  3.2  shows  the  results  for  the  accelerations 
at  midspan  which  are  computed  by  using  different  numbers 

of  Y   functions.   Since  the  Y   functions  have  been  used 
n  n 

to  specify  the  transverse  configuration  of  the  bridge, 
the  larger  number  of  n,  should  provide  better  results. 
The  parameters  defining  the  bridge  and  vehicle  system 
are  the  same  as  before.   It  can  be  seen  that  the 


^:J' ly*--"  .  :••-•  ifi-T 
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Table  3.2    Comparison  of  Accelerations  Obtained  by  Using 
Different  Numbers  of  nx  with  Wheels  Over 
Center  Beam 

Vehicle  at  center  beam  of  the  bridge,  N=  400 


'L 

Accelerations  at 

Midspan 

^/ 

"l 

Edge  Beam 

Beam  No.  2 

Center  Beam 

2 

-0.274466 

-0.274466 

-0.274466 

0. 

,20 

4 

-0.252714 

-0.298451 

-0.217295 

6 

-0.269274 

-0.294369 

-0.210257 

2 

-0.348272 

-0.348272 

-0.348272 

0. 

,50 

4 

•  -0.160029 

-0.411076 

-0.515064 

6 

-0.167445 

-0.406655 

-0.517048 

2 

-0.293332 

-0.293332 

-0.293332 

0. 

.80 

4  • 

-0.269257 

-0.329532 

-0.354500 

6 

-0.262305 

-0.338745 

-0.363971 

2 

0.090615 

0.090615 

0.090615 

1. 

,00 

4 

0.115142 

0.155268 

0.171888 

6 

0.118075 

0.161762 

0.176689 
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accelerations  converge  when  six  terms  of  the  Y   functions 

n 

are  used  (n^  =  6) .   Table  3.3  is  similar  to  Table  3.2  by 
showing  the  accelerations  at  midspan  converge  with  the 
increasing  values  of  n, ,  except  that,  at  this  time, 
the  vehicle  travels  over  the  edge  beam  of  the  bridge  and 
the  positions  of  wheels  are  defined  by  the  parameters: 
n-|^  =  0.10,  n^  =  0.40. 

For  the  numerical  results  presented  in  this  chapter, 
the  following  values  of  solution  parameters  are  used: 
N  =  400,  m^  =  15,  n^  =  8,  n  =  6,  and  n-,  =  6 . 

3.3   Bridge  Parameters 

3.3.1   Effect  of  Span  Length 

Several  previous  field  tests  have  indicated  that 
the  span  length  is  one  of  the  major  variables  affecting 
the  vibration  of  the  bridge.   In  this  study,  the  simple- 
span  bridge  models  with  different  span  lengths  are  used 
in  the  computation  to  evaluate  the  accelerations  of  each 
beam  at  midspan  while  keeping  the  vehicle  parameters 
constant. 

In  order  to  obtain  practical  values  for  bridge 
parameters,  the  Bureau  of  Public  Roads  bridges     with 
span  lengths  in  the  range  between  20  and  70  ft.  are 
considered.   These  bridges  are  of  the  I-beam  type  and  are 
designed  either  for  H  15-44   loading  or  for  HS  20-44 
loading.   The  bridges  designed  for  H 15-44  loading  have  a 
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Table  3.3    Comparison  of  Accelerations  Obtained  by  Using 
Different  Numbers  of  n^  with  Wheels  Over 
Edge  Beam 

Vehicle  at  left  edge  beam  of  the  bridge,  N  =  400 


x/L 

"l 

Acce 

ilerations  at  Midspan 

Left  Edge  Beam 

Right  Edge  Beam 

Center  Beam 

2 

-0.566199 

0.011037 

-0.277580 

0.20 

4 

-0.556376 

0.090094 

-0.303591 

6 

-0.558713 

0.087076 

-0.305147 

2 

-0.125344 

-0.499783 

-0.312563 

0.50 

4 

-0.118537 

-0.424207 

-0.340399 

6 

-0.112022  '■ 

-0.430594 

'■   -0.340199 

2 

0.011045 

-0.616277 

-0.252910 

0.80 

4 

0.208124 

-0.643312 

-0.272986 

6 

0.207292 

-0.634261 

-0.269409 

2 

0.481108 

-0.234304 

0.123401 

1.00 

4 

0.466486 

-0.164051 

0.192613 

6 

0.473863 

-0.169914 

0.195106 
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roadv/ay  width  of  2  8  ft.  and  a  concrete  slab  7  in.  thick. 
The  entire  deck  is  supported  on  four  wide-flange  steel 
girders.   For  the  bridges  designed  for  the  HS  20-44 
loading,  the  roadv.'ay  width  is  44  ft.,  and  the  slab 
thickness  is  7-1/2  in.,  and  a  total  of  six  wide-flange 
steel  girders  are  used  to  support  the  deck.   The  cross- 
section  of  both  bridge  types  are  shown  in  Figure  3.1. 
It  should  be  noted  that  the  effects  of  side  curbs  are  not 
considered  in  this  study. 

The  weight  per  unit  of  length  of  bridge  was  evaluated 
by  computing  the  weight  per  unit  of  length  of  an  interior 
girder  and  its  tributary  slab  and  multiplying  the  result 
by  the  number  of  girders.   The  fundamental  frequency  of 
vibration  of  each  bridge  was  computed  by  using  Eq.  (2-44) 
in  Article  2.2,  and  the  period  of  vibration  was  obtained 
by  inverting  the  corresponding  frequency.   These  results 
for  both  bridge  types  are  presented  in  Table  3.4.   It 
was  assumed  that  the  roadway  surface  was  horizontal  and 
smooth. 

The  vehicle  is  represented  by  a  single  axle  loading 
consisting  of  two  wheels  and  has  total  sprung  load  of  30 
kips  for  II  15-44  loading  and  of  72  kips  for  HS  20-44 

loading.   The  stiffness  of  each  tire  spring  for  both  vehicle 

(22) 
types  is  6  kips  per  inch     and  the  weight  of  unsprung 

loads  or  wheels  are  neglected.   From  this  information,  the 

frequency  of  vehicle  can  be  computed  by  using  Eq.  (2-44) 
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in  Article  2.2.   The  positions  of  both  v/heels  while  the 
vehicle  travels  on  the  bridge  are  shown  in  Figure  3.1, 
where  the  right  wheel  is  over  beam  no.  2  and  the  trans- 
verse spacing  of  the  v/heels  is  6  feet. 

The  maximum  acceleration  results  at  midspan  of  both 
B.P.R,  bridge  types  with  varying  span  lengths  are  evaluated 
and  shown  in  Table  3.5.   The  maximum  acceleration  of  each 
beam  is  shown  and  the  corresponding  value  of  x/L^ 
defines  the  longitudinal  position  of  the  wheels  at  that 
time.   Figure  3.2  shows  graphically  the  results  corres- 
ponding to  Table  3.5,  each  curve  representing  the  maximum 
acceleration  of  one  beam.   It  can  be  seen  that  the 
acceleration  increases  with  shortening  span  and  the 
rate  of  increase  for  spctn  lengths  between  20  ft.  and  40  ft. 
is  higher  than  that  for  span  lengths  between  40  ft.  and 
70  ft.   In  other  words,  the  acceleration  of  the  bridge 
decreases  if  the  span  length  of  the  bridge  increases. 
It  should  also  be  noted  that,  in  Figure  3.2,  the  maximum 
accelerations  of  the  B.P.R.  bridges  with  44  ft.  roadway 
width  are  greater  than  those  with  28  ft.  roadway  width. 
It  can  then  be  concluded  that  the  acceleration  increases 
for  heavier  vehicle  loads  on  the  bridge. 

Figure  3.3  shows  the  distributions  of  maximum  accelera- 
tion to  each  beam  corresponding  to  the  results  shown  in  Table 
3.5  and  Figure  3.2.   For  both  B.P.R.  bridge  types,  the  right 
wheel  is  over  the  second  beam  from  the  left  (beam  2) .   This 
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20.0     30.0     40.0     50.0 
SPAN  LENGTH  (FT J 


O  BEAM  1 
A  BERM  2 
+  BEAM  3 


60.0 


70.0 


28  FT  RQPDWflY   H15-44  LOADING  VELOCITY  60  MPH   LEVEL  SURFACE 
TWO  WHEELS   6  FT  WHEEL  SPACING   RIGHT  WHEEL  OVER  SECOND  BEAM 


Figure  3.2    Maximum  Accelerations  at  Midspan  of  B.P.R. 
Bridges 
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Figure  3.2,  cont, 
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Figure  3.3    Distributions  of  Maximum  Acceleration  to 
Each  Beam  of  28  ft  Roadway  B.P.R.  Bridges 
with  Different  Span  Lengths 
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causes  the  acceleralions  of  the  exterior  beiuns  to  be  larger 
than  tliose  of  the  interior  beams.   The  differences  in 
accelerations  between  these  beams  decrease  as  the  bridge 
span  increases. 

Figure  3.4  shov;s  the  acceleration  of  each  beam  in 
the  form  of  history  curves  for  a  bridge  having  a  60  ft. 
span  and  a  44  ft.  roadway  width  with  the  right  wheel 
of  the  vehicle  over  beam  2.   A  history  curve  is  a  plot  of 
acceleration  as  a  function  of  time,  or  in  tenns  of  longi- 
tudinal position  of  the  load  on  the  bridge.   It  can  be 
seen  that  accelerations  of  beam  1  and  beam  6  are  out  of 
phase  all  the  time.   It  is  of  some  interest  to  note  also 
that  the  maximum  acceleration  of  beam  1  occurs  when  the 
vehicle  enters  the  bridge  but  the  maximum  acceleration  of 
beam  6  occurs  when  the  vehicle  is  at  the  midspan  or  is 
leaving  the  bridge. 

Figure  3.5  shows  the  history  curves  of  acceleration 
of  beam  6  with  varying  span  lengths  corresponding  to  the 
results  presented  in  Table  3.5.   It  can  be  seen  that  the 
curves  are  similar  except  that  larger  ordinate  of  the 
curve  for  the  sliorter  span  lengths. 

3.3.2   Effect  of  Width 

In  order  to  investigate  the  effect  of  bridge  width 
on  acceleration,  three  bridge  models  as  siiovm  in  Figure  3.6(a) 
have  been  assumed.   All  three  bridges  have  the  same  parameters 
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except  the  width  and  the  number  of  beams.   Beam  1  and  beam 
4  are  designated  as  the  exterior  beams  and  beams  2,  3,  5  and 
6  as  the  interior  beams.   Vehicle  parameters  are  constant 
for  this  comparison.   The  study  is  divided  into  two  cases. 
In  case  A,  the  vehicle  travels  at  the  center  line  of  the 
bridges  as  shown  in  Figure  3.6(a),  while  the  left  wheel 
is  over  beam  1  for  case  B.   The  results  of  the  study  and 
descriptions  of  all  parameters  are  summarized  in  Table  3.6 
and  the  corresponding  history  curves  of  beam  1  are  shown 
in  Figure  3.7.   It  can  be  seen  from  these  results  that 
there  is  not  much  difference  in  beam  acceleration  with 
increasing  bridge  widths  although  the  accelerations  of 
narrower  bridges  tend  to  be  slightly  higher  than  those  of 
the  wider  bridges. 

3.3.3   Effect  of  Flexibility  of  Girder 

Figure  3.G(b)  shows  the  bridge  model  used  for  a  study 
to  obtain  the  accelerations  for  varying  girder  flexibilities 
v;hile  kcepinc;  other  parameters  constant.   The  transverse 
vehicle  position  is  arranged  so  that  its  left  wheel  is 
over  beam  1  at  all  times.   The  details  of  necessary  parameters 
used  in  this  study,  including  six  different  girders,  are 
presented  in  Table  3.7.   It  should  be  remembered  that  all 
five  girders  of  the  bridge  have  the  same  section.   The 
36  WF  230  v/ide  flange  girders  v/ould  be  the  normal  design  for 
this  bridge. 
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Table  3.6    Maximum  Accelerations  at  Midspan  of 

Simple  Span  Bridges  with  Different  Widths 

Girder  Bridge    60  ft. span    7.5  in.  Deck  Thickness 
36WF230    Girder  8  ft. Spacing    HS  20-44  Loading 
60  m.p.h.    Level  Surface 


Case  A     Load  at  Center  Line 


Maximum  Accel 

er 

ations 

at 

Midspan 

(in/sec^) 

Beam 

24 

ft  Width 

x/L 

32 

ft  Width 

x/L 

40  ft  Width 

x/L 

1 

-71.4  7 

0.50 

-69.76 

0.84 

-64.18 

0.25 

4 

-71.47 

0.50 

-69.76 

0.84 

-64.18 

0.25 

2 

-66.69 

0.81 

-49.79 

0.25 

-47.70 

0.23 

3 

-66.69 

0.81 

-49.79 

0.25 

-47.70 

0.23 

5 

70.65 

1.00 

-65.50 

0.54 

6 

-65.50 

0.54 

Case  B     Left  Wheel  Over  Beam  1 


Beam 

Max: 

imum  Acce 

ilerations 

at 

Midspan 

(in/sec^) 

24  ft  Width 

x/L 

32  ft  Width 

x/L 

40  ft  Width 

x/L 

1 

125.84 

0.98 

100.81 

0.05 

100.67 

0.05 

4 

120.15 

0.39 

119.14 

0.68 

107.20 

0.44 

2 

-80.53 

0.84 

-72.18 

0.20 

65.11 

0.07 

3 

-76.92 

0.51 

-69.97 

0.56 

62.89 

0.73 

5 

-56.91 

0.85 

-62.70 

0.24 

6 

59.82 

0.73 

100 
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Y       I^        5 


..1 8'. -.-I 
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II  I.-  f  13  f 
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:•      I'      r^      i^     13      f 
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(a)    Effect    of    Width     on     Acceleration 


I'  '"      I'  5  I'  I^ 


32' 


(b)      Effect    of    El     on     Acceleration 


Figure  3.6    Cross  Section  of  Bridges  Used  in  Study 
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Five  gj.rder  sections  smaller  than  the  36  WF  230  section 
are  used  and  the  smallest  section  (36  WF  135)  has  a  moment 
of  inertia  nearly  half  that  of  36  WF  230  section.   The 
results  of  the  study  are  tabulated  in  Table  3.7  and  the 
relationships  between  the  acceleration  of  each  beam  and 
the  moment  of  inertia  of  beam  are  plotted  in  Figure  3.8 
The  fundamental  frequency  and  the  total  weight  of  the  bridge 
for  each  type  of  girder  are  also  calculated  and  presented 
in  Table  3.7.   From  these  results,  as  might  be  expected, 
the  acceleration  is  lowest  when  the  largest  size  of  sec- 
tion has  been  used  and  it  is  higher  when  the  size  has 
been  reduced.   It  can  be  seen  that  the  accelerations  increase 
appro>:im<i toly  only  20  percent  v;hile  the  moment  of  inertia 
of  section  is  reduced  nearly  40  percent. 

3 . 4   Vehicle  Parameters 

3.4.1   Effect  of  Transverse  Position  of  Wheels 

The  }5ridge  model  used  for  this  comparison  is  the  same 
as  that  used  in  Article  3.3.3  and  shown  in  Figure  3.6(b). 
A  36  WF  230  girder  section  is  used.   Figure  3.9  shows  the 
transverse  position  of  the  wheels  on  the  bridge  in 
six  different  cases.   Case  A  has  the  leftmost  position  of 
wheels  on  the  bridge  while,  in  case  F,  the  wheels  are 
symmetrically  arranged  over  the  center  beam. 

Table  3.8  shows  the  acceleration  results  for  each 
case,  giving  the  maximum  acceleration  of  each  beam  and  the 
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32' 


1}       1=        I'        1^ 


Case    A      left   wheel    over 


edge    beam 


Case    B       Right    wheel     over 
beam  2 

Case    C       Symmetrical    wheels 
over    beam  2 

Case    D       Left     wheel    over 
beam  2 

Case   E       Right    wheel    over 
center    beam 

Case   F       Symmetrical    wheels 
over    center    beam 


Figure  3.9    Transverse  Position  of  Wheels  on  Bridge 
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longitudinal  position  of  the  vehicle  at  that  time.   Tlie 
details  of  all  necessary  parameters  are  also  provided 
in  Table  3.8.   The  relationships  between  the  acceleration 
of  beam  and  wheel  position  are  plotted  in  Figure  3.10. 
There  are  three  important  things  to  be  noted  from  these 
results.   First,  the  accelerations  of  edge  beams  are  the 
greatest  when  the  vehicle  travels  over  the  edge  beam  and 
tend  to  decrease  when  the  vehicle  travels  near  the  center 
line  of  the  bridge.   Second,  in  contrast,  the  acceleration 
of  the  center  beam  is  the  least  when  the  vehicle  travels 
over  the  edge  beam  and  the  greatest  when  the  vehicle  is  over 
the  center  line  of  the  bridge.   Center  beam  accelerations 
are  slightly  greater  than  the  accelerations  of  edge  beams. 
Third,  practically,  most  vehicles  travel  on  the  bridge 
with  a  wheel  position  corresponding  to  case  C  or  case  D 
and  the  accelerations  of, the  edge  beams  are  then  greater 
than  those  of  interior  beams . 

3.4.2   Effect  of  Number  of  Wheels 

Table  3.9  presents  the  results  of  acceleration  compari- 
son obtained  by  using  a  single  axle  having  one  or  two 
v/heels.   Both  vehicle  models  have  the  same  total  load  of 
72  kips.   As  shown  in  Tabic  3.9,  tliere  are  three  different 
transverse  positions  of  load  on  the  bridge:   first,  the  load 
is  over  the  edge  beam,  second,  the  load  is  over  beam  2  and 
third,  the  load  is  over  the  center  beam.   The  bridge 
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parameters  are  kept  constant  and  are  described  in  Table  3.9. 
It  can  be  seen  from  the  results  that  the  accelerations 
obtained  by  using  the  single  wheel  vehicle  are  greater 
than  those  obtained  by  using  the  two  wheel  vehicle.   It 
may  bo  noted  that  the  difference  in  accelerations  is  approxi- 
mately 10  percent   when  the  load  is  over  beam  2  and  over 
the  center  beam  but  the  acceleration  difference  ranges 
from  25  to  48  percent   when  the  load  is  over  edge  beam. 
This  can  be  explained  in  that  the  center  of  gravity  of 
two  wheel  vehicle  is  not  actually  over  the  edge  beam  as 
in  the  case  of  one  wheel  vehicle.   The  10  percent 
difference  may  be  expected  if  the  two  wheel  vehicle  has  its 
center  of  gravity  over  the  edge  beam. 

3.4.3   Effect  of  Transverse  Position  of  Single  Load  and 
Flexibility  of  Girder 
In  Article  3.3.3,  the  effect  of  girder  flexibility 
on  accelerations  was  studied  for  a  constant  position  of 
wheels  over  the  edge  beam  and  it  was  shown  that  reduced 
girder  section  had  relatively  small  effects  on 
acceleration.   In  this  study,  the  transverse  position  of 
single  load  on  the  bridge  is  varied  and  it  is  classified 
into  five  different  cases  as  shown  in  Figure  3.11.   In  each 
case,  the  comparisons  are  made  for  the  same  six  girder 
sections  as  in  Article  3.3.3  (properties  of  these  sections 
are  shown  in  Table  3.7).   In  Table  3.10,  the  remaining 
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parameters  are  summarized  and  the  results  of  the  study  in 
each  case  are  presented.   Figure  3.12  shov/s  the  relationship 
betv/ecn  maximum  acceleration  of  each  beam  and  the  variation 
of  inomc^nt  of  inertia  of  the  section  for  all  five  different 
load  positions  on  the  bridge.   It  can  be  seen  graphically 
that,  in  all  cases  except  case  A,  the  acceleration  is  slightly 
increased  when  the  section  of  beam  is  reduced.   To  be 
specific,  the  acceleration  is  not  increased  more  than  15 
percent  v;hile  the  moment  of  inertia  of  section  is  reduced 
over  30  percent.   In  case  A,  where  the  load  is  over  the  edge 
beam,  the  accelerations  of  the  edge  beams  increase 
significantly  v/hen  the  section  is  reduced.   However,  vehicles 
travelling  at  60  m.p.h.  tend  to  be  in  the  lanes  near  the  center 
line  of  bridge. 

The  curves  in  Figure  3.13  show  graphically  the  varia- 
tion of  the  maximum  accelerations  of  edge  beams  for  different 
transverse  load  positions  with  the  relative  flexibility 
of  beam.   For  a  position  of  60  m.p.h.  vehicle,  the  increased 
beam  flexibility  affects  the  maximum  acceleration  slightly 
if  the  moment  of  inertia  of  the  section  is  not  reduced  more 
than  3  5  percent  from  the  normal  size. 

3.4.4   Effect  of  Speed 

The  bridge-vehicle  system  used  in  this  article  is  the 
same  as  that  used  in  Article  3.3.3  and  has  been  shown  in 
Figure  3.6.   The  maximum  acceleration  of  each  beam  is 
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Figure  3.12,  cont. 


120 


CE 

to 
a 


cr 


en 


O  BEAM  1 
A  BEAM  2 
+  BEAM  3 
X  BEAM  4 
O  BEAM  5 


7000 


9000 


11000 


13000 


15000 


5  GIRDER  BRIDGE  60  FT  SPAN  72  KIPS  SINGLE  LOAD  LEVEL  SURFACE 
CASE  DJ   SINGLE  LQAD  HALFWAY  BETWEEN  BEAM  2  AND  BEAM  3 


Figure  3.12,  cont. 


121 


o. 


UJ  o  — 


tx 
a. 

o 

CO 

a 

o 

in 

n 

h- 

CL 

■21 

tj 

o 

»— i 

t— 

Q 

CT 

o 

Qi 

-^ 

UJ 

_J 

UJ 

O 

O 

cr 

O  BEAM  1 
A  BEAM  2 
+  BEAM  3 
X  BEAM  H 
<^  BEflH  5 


1 — 

13000 


7000 


9000 


11000 


15000 


5  GIRDER  BRIDGE  60  FT  SPAN  72  KIPS  SINGLE  LOAD 
CASE  E)   SINGLE  LOAD  OVER  BEAM  3  (CENTER  BEAM) 


LEVEL  SURFACE 


Figure  3.12,  cont, 


122 


o_ 


CM 

o 


rr 

n. 

<_3 

to 
a 

O 
n 

»— 1 

r\j 

i- 

1— 

(X 

t: 

o 

o 

t-t 

h- 

o 

cr 

to 

a: 

^^ 

LU 

_j 

UJ 

CJ 

o 

rr 

a 

■V 

• 

o 
o 

»— • 

X 

CL 

BEAM  5 


-+- 


-H 


-H 


O  36WFE30 
A  36WF194 
+  36WF102 
X  36WF170 
O  36WF160 
36WF135 


-h 


-+- 


BEAM  1  BEAM  2 

TRANSVERSE  PQSITIQN  QF  WHEELS 


BEAM  3 


5  GIRDER  BRIDGE 
LEVEL  SURFACE 


60  FT  SPAN  SINGLE  LOAD  72  KIPS  60  MPH 


Figure  3.13    Effect  of  Flexibility  of  Girder  and  Transverse 
Position  of  Load  on  Acceleration 


123 


en 

Q- 

a 


o  _ 


I— 

o 

en 

lO 

(H 

•— • 

iij 

_) 

LU 

LJ 

tJ 

(T 

o 

V 

' 

— ( 

o 

i: 

•-4 

H- * 

X 

(X 

BERM  1 


O  36WF230 

A  36WF194 

+  36WF182 

X  3GWF17Q 

O  36WF160 

4^  36WF135 


BEAM  1  BEAM  2 

TRANSVERSE  POSITION  QF  WHEELS 


BEAM  3 


5  GIRDER  BRIDGE 
LEVEL  SURFACE 


60  FT  SPAN  SINGLE  LOAD  72  KIPS  60  MPH 


Figure  3.13,  cont, 


124 


evaluated  for  vehicle  speeds  ranging  from  20  mph  to  70  mph. 
The   values  of  the  remaining  parameters  are  the  same  as 
previously  and  are  presented  in  Table  3.11.   The  numerical 
acceleration  results  with  varied  speeds  of  vehicle  are 
also  presented  in  Table  3.11  and  the  corresponding  curves 
are  plotted  in  Figure  3.14.   As  might  be  expected,  the 
values  of  the  maximum  accelerations  increase  for 
vehicles  with  higher  speeds.   It  should  be  noted  also  that 
the  rate  of  change  of  acceleration  in  each  beam  is 
approximately  the  same  as  that  of  the  vehicle  speed. 

3 . 5   Surface  Roughness  Parameter 

Several  previous  test  reports  ^"^  ^  ^^^^  ^"^^  ^  have 
indicated  that  the  surface  roughness  is  one  of  the  most 
important  parameters  which  affect  the  vibration  of  highway 
bridges.   These  reports  have  recommended  that  the  bridge 
surface  should  be  as  smooth  as  possible.   In  this  study, 
the  effect  of  the  surface  roughness  on  accelerations  is 
investigated  and  the  corresponding  numerical  results 
are  presented. 

The  bridge  models  used  in  the  previous  articles  in 
this  chapter.  Article  3.2  through  Article  3.4,  were  assumed 
to  have  a  level  surface.   In  this  article,  the  same 
bridge  and  vehicle  models  as  used  in  Article  3.3  and 
shown  in  Figure  3.6  are  used  but,  this  time,  the  initial 
bridge  surface  is  rough.   The  parameters  representing 
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the    surface   roughness   have   been   designated   as   w,  ,    v/    ,    in 
^  ^  loo 

and  e,  which  are  defined  in  Article  2,2.3.   The  term  w, 
represents  the  deviation  of  the  deck  of  the  bridge  and  has 
been  defined  by  Eq .  (2-24),  the  term  w^e .  denotes  the 
amplitude  of  unevenness  and  the  term  m   indicates  the 
number  of  half  sine  waves.   It  should  be  noted  that  the 
shape  of  bridge  surface  is  ideolized  as  a  number  of  half 
sine  waves  with  amplitudes  at  both  supports  of  the  bridge 
equal  to  zero.   The  amplitude  of  roughness  can  be  upward  or 
downward  depending  on  the  sign  of  the  we.  quantities. 

Throughout  the  study,  it  is  assumed  that  the  shape  of 
the  bridge  roughness  surface  and  its  amplitude  under  both 
v/heels  of  vehicle  are  the  same.   The  study  of  the  effect 
of  the  roadway  unevenness  parameters  on  accelerations  can 
be  classified  into  two  cases.   First,  varying  numbers  of 
half  sine  waves  were  studied  while  keeping  the  amplitudes 
constant,  and  second,  roughness  amplitudes  were  varied  while 
keeping  the  numbers  of  half  sine  waves  constant.   The  remain- 
ing parameters  of  the  bridge  and  vehicle  are  considered 
to  be  the  same  as  before. 

Table  3.12  shows  the  numerical  results  for  accaleration: 
obtained  by  using  different  numbers  of  half  sine  waves  to 
represent  tl^e  roadway  unevenness  of  the  bridge.   A  number 
of  half  sine  waves  equal  to  zero  indicates  that  the  bridge 
has  a  level  surface  and  the  larger  number  of  half  sine 
waves  indicates  that  the  bridge  surface  is  rougher.   The 


128 


0 

4J 

<U 

0 

a 

• 

X 

(U 

• 

en 

Ch 

•n 

• 

•H 

B 

M 

CQ 

0 

VX3 

c 

(« 

a. 

CO 

Cn 

K   10 

CJ 

•H    Ifi 

r-l 

T5   0) 

a 

tl   C 

R 

0  x: 

•H 

ij  tp 

w 

3 

en 

rr   0 

M-l    0) 

^  Ci 

o  > 

1 

nj 

0  i+-( 

C  S 

o)   0 

fO 

Q4  0) 

CO    dJ 

U)    C 

K  Tl 

TI  -H 

3 

■  H  CO 

4J 

■•-t 

*"  iW 

0    rH 

4J  .H 

rn   O. 

CM   B 

r:  ^ 

1X1  x: 

0  0 

m   u 

■H 

C 

-P    W 

-H 

(0   u 

u  <u 

C  Ln 

0)  J3 

m    • 

•H    G 

0,0 

OJ    3 

CO 

U  ^ 

. 

0 

+j 

-<  -u 

y-i  r^ 

c 

e  0) 

0  e 

3  j^ 

VD   nj 

e  <u 

0) 

■H  M-l 

m 

X  U-i 

(fl    •r^ 

0)   M 

s:  Q 

en  (1) 

T3    > 

■H  0 

V^ 

rsi 

pa  .H 

r-i 

0) 

• 

^^    0) 

ro 

0)  x: 

•d  £: 

0) 

>-l 

tH 

■H  -P 

XI 

0  ^ 

(0 

0) 

tH 

m  i_q 

u 

■H 


Q) 

CQ 

m  (u 

0  c 

•H 

M  CO    w 

0)          0) 

XI  •+-!    > 

6  -H    (0 

3  t«  S: 

z  a 

H 

oooocTiOCTirocnQO'xiO'xir-^ 
ix>v£)U3r^|---criix)cricrir-m»x) 


000000000000 


•>3'a>>JDOOrHLn'^'Lnco'*or) 
oor^aDcri^ojOiHinoootn 

I     I  I 


v£)V£)rOrHn(N'*no>x>'^Ln 
LnLnoor~~criCT\oooor^Ln'«:}'m 


000000000000 


r~-vDrHr^r-~iHtMiH.HO-<^o 

OrO'^Ta'crirMCNiHf^'NiHO 

oorHLn<rioor--tNr~orNi'^ 
r^r-~r-~rH(NOooooooocr\ 

I        I        I      iH   rg    rH      I      rH    rH    ^  I 

I  I  I 


oocNcnoavDoooocoo^incoLn 


000000000000 


oo>i)crinooLr)<NOOcHr-H 
(TimfMU3in<x>rsironrHv£)cr< 

usvocNcn^mcHCNnrHCTi';!' 
LnLni^ocricrirjcrirnrHv£i>^ 

I        I        I      rH  I      (N    (N    fN    rH      I        I 

I  rill 


oooocTiroovo'^fMocri 
(N<NCM(NVD<Xiin'<TrocMrvJ.H 


000000000000 


nr^voous-^rroOfHroinr^ 
rHvoTLnoocornrHor-crvr^ 

(NCNrHVXl^rMr^OM'rHrrio 

r^r-~cocriooMD(N>Hcricrioo 

I         I         I         1      rH    rH    M    rH  I      rH    rH 

I         I         I         I  II 


LnLnLninrorH'^OiH<X)'a'M3 
ooomCT\cx5a>crir^ir)^a'oo 


000000000000 


nooOLn(NrH"^Ln(X)CMcx)r- 
t^ir)r^oorHr^rocNr--(Ncocr> 

otNinoDr^vDnmcrirHMDoo 
ooO'^cr\0'*vx>rHin>HLn 

iH«HtHrHCNCNrO'^-^OJ(NtH 
I       I       I       I      I  I 


OiHroinr^cTirHt^jnuni^iTi 

rH   .H  .H  .-I   .H  rH 


129 


amplitude  of  roughness  measured  at  the  top  of  the  waves  is 
assumed  to  be  constant  and  has  the  value  of  0.5  inches. 
The  number  of  half  sine  waves  used  ranged  from  0  to  19 
and  the  corresponding  maximum  accelerations  together  with 
the  longitudinal  position  of  vehicle  when  it  happens  are 
listed.   These  results  are  also  plotted  in  Figure  3.15. 

It  can  be  seen  from  these  results  that  the  accelera- 
tions were  not  influenced  by  the  effect  of  the  surface 
roughness  for  less  than  three  half  sine  waves.   This  effect 
increases  markedly  when  the  number  of  half  sine  waves 
is  greater  than  five  and  reaches  its  peak  when  the  bridge 
surface  consists  of  twelve  half  sine  waves.   The  maximum 
results  affected  by  the  number  of  half  sine  waves  are  as 
much  as  five  times  the  values  of  bridge  with  level  surface. 
It  should  also  be  noted  that  the  accelerations  of  edge 
beams,  designated  as  beam  1  and  beam  5,  are  greatest  for 
every  number  of  half  sine  waves. 

Figure  3.16  shows  the  comparison  of  the  history 
curves  of  interacting  force  between  the  bridge  and  the 
vehicle  with  three  different  number  of  half  sine  waves, 
the  level  surface,  seven  and  twelve.   The  results  of  the 
interacting  forces  are  in  the  terms  of  static  v/heel  force 
and  correspond  to  the  results  in  Table  3.12.   It  can  be 
seen  that  the  interacting  forces  on  the  level  surface 
bridge  are  very  much  different  from  the  rough  surface 
bridge.   The  periods  of  the  interacting  forces  with  seven 
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and  twelve  half  sine  waves  are  equal  to  the  first  and 
second  natural  periods  of  bridge  vibration  respectively. 

Figure  3.17  shows  the  comparison  of  the  history  curves 
representing  the  acceleration  of  beam  5  with  three 
different  numbers  of  half  sine  waves  corresponding  to  the 
results  shown  in  Figure  3.16  and  Table  3.12.   All  curves 
are  very  similar  when  the  longitudinal  position  of  the 
vehicle  is  in  the  first  quarter  of  the  bridge  length. 
It  is  then  seen  that,  for  the  rest  of  the  bridge  length, 
the  curves  representing  the  rougher  bridge  will  have  more 
cycles  and  greater  amplitudes  of  acceleration.   This  is 
caused  by  the  resonance  between  the  bridge  and  the  vehicle. 

The  remaining  parameter  which  also  represents  the 
roadway  unevonness  is  the  amplitude  of  roughness.   The 
practical  range  of  this  parameter  is  between  zero  to  one 
inch.   The  study  of  the  acceleration  affected  by  the 
varied  amplitudes  of  roughness  is  classified  into  two 
different  cases.   In  case  A,  a  bridge  profile  consisting 
of  three  half  sine  waves  which  had  little  effect  on  the 
acceleration   is  used,  while  in  case  B,  twelve  half  sine 
waves  having  the  most  effect  are  used.   In  each  case, 
five  different  amplitudes  of  roughness  have  been  used  to 
evaluate  the  results.   These  five  amplitudes  have  the 
values  of  0.,  0.25,  0.50,  0.75  and  1.0  inches.   The 
remaining  parameters  are  considered  to  be  the  same  as 
before . 
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The  numerical  results  in  both  cases  are  presented  in 
Table  3.13  and  the  corresponding  curves  are  plotted  in 
Figure  3.18.   From  these  results,  it  can  be  seen  clearly 
that  the  effect  of  the  amplitude  depends  on  the  number 
of  half  sine  waves.   For  the  increased  values  of  roughness 
amplitude,  the  acceleration  of  each  beam  does  not  change 
at  all  in  case  A  but  it  increases  as  much  as  ten  times  from 
the  original  value  with  the  level  surface  in  case  B. 
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CHAPTER  IV 
ACCELERATION  STUDIES  OF  THREE  SPAN  CONTINUOUS  HIGHWAY  BRIDGES 

4 . 1   General 

The  purpose  of  this  chapter  is  to  study  the  major 
parameters  which  affect  the  accelerations  of  three-span 
continuous  highway  bridges  subjected  to  moving  vehicles. 
The  method  of  analysis  and  the  description  of  the  computer 
program  used  in  this  study  have  been  presented  in  Articles 
2.3.2  and  2.3.4  respectively.   In  order  to  be  able  to 
compare  the  investigated  accelerations  of  the  bridge 
model  with  the  available  acceleration  criterion  for 
human  response,  all  input  data  and  output  results  are  in 
dimensional  form  and  in  a  practical  range  of  bridge 
parameters.   The  parameters  can  be  classified  into  the 
following  four  groups: 

a)  Solution  parameters  to  satisfy  the  convergence  and 
stability  of  the  solution. 

b)  Bridge  parameters  including  the  span  length, 
the  ratio  of  side  span  to  center  span  and  the 
flexibility  of  girder. 

c)  Vehicle  parameters  including  the  vehicle  load, 
the  axle  spacing  and  the  speed  of  vehicle. 
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d)  Initial  conditions  of  vehicle  and  bridge  including 
the  initial  oscillation  of  vehicle  and  the  surface 
roughness  of  the  bridge  deck. 

Throughout  this  chapter,  the  bridge  is  assumed  to  be 
an  I-beam  type  composed  of  steel  girders  and  a  reinforced 
concrete  deck  v/ith  non-composite  action.   It  has  a  pris- 
matic cross-section  and  equal  side  spans.   Damping  is 
considered  in  both  bridge  and  vehicle.   The  bridge  and 
vehicle  models  are  referred  to  Figures  2.9  and  2.11 
respectively.   Most  of  the  vehicle  parameters  such  as 
the  coefficients  of  friction  and  frequencies  of  the  tire- 
suspension  system,  weight  ratios  of  tire  to  the  vehicle 

load  were  obtained  from  the  information  presented  in  Table 

(29) 
4.1.^      Table  4.1  shows  the  characteristics  of  "typical" 

vchicJ.es  including  three  axle  vehicles,  two  axle  vehicles 

and  single  axle  loading.   These  characteristics  are  average 

values  and  were  obtained  from  information  given  in  Reference 

2  and  from  manufacturers'  data.   The  vehicle  is  assumed 

to  travel  from  left  to  right  with  constant  velocity.   The 

calculated  acceleration  of  the  bridge  occur  at  specified 

node  points  of  the  prismatic  beam  as  shown  in  Figure  2.9. 

The  maximum  values  of  acceleration  together  with  the 

corresponding  position  of  axle  at  that  time  are  considered. 

The  natural  frequencies  of  the  continuous  beam  used 

to  approximate  the  actual  bridge  can  be  evaluated  by 

obtaining  the  value  of  the  dimensionless  parameter  X  . 
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Table  4.1    Data  for  "Typical  Vehicles 


(29) 


Quantity 

Unit 

Three-axle 

Two-axle 

Sinig 

[le-axle 

Vehicle 

Vehicle 

Vehicle 

W 

kips 

72 

64 

64 

h 

.  ^ft..  ,  ■ 

12 

14-35 

- 

^2 

ft. 

14-35 

- 

- 

W^/W 

0.08 

0.90 

- 

w^/w 

0.80 

- 

- 

w^/W 

0.03 

0.05 

- 

w^/W 

0.05 

0.05 

- 

W3/W 

0.04 

- 

- 

il 

0.5-1.0 

0.9-1.7 

- 

^2 

0.9-1.7 

- 

- 

^1 

0.602 

0.5 

- 

^3 

0.494 

- 

- 

^5 

0.083 

- 

- 

Ml- 

0.05-0.10 

0.12-0.28 

0. 

,12-0.28 

•^2 

0.12-0.28 

0.12-0.28 

- 

^3 

0.12-0.28 

- 

- 

't,l 

ops 

3.13-3.72 

3.13-3.72 

3. 

,13-3.72 

^,2 

cps 

3.13-3.72 

3.13-3.72 

- 

't,3 

cps 

3.13-3.72 

- 

- 

^ts,l 

cps 

1.57-1.65 

1.99-2.14. 

1. 

,74-2.36 

'ts,2 

cps 

1.99-2.14 

1.74-2.36 

- 

^ts,3 

cps 

1.74-2.36 

- 

- 

^t,l 

kips 

8 

32 

64 

^st,2 

kips 

32 

32 

- 

^st,3 

kips 

32 
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which  is  plotted  as  a  function  of  span  ratio  "a",  in 
Figure  4.1.   The  subscript  j  is  an  integer  denoting  the 
order  of  the  natural  frequency  under  consideration.   These 
data  were  evaluated  by  application  of  the  method  described 
briefly  in  Reference  36.   The  natural  frequency  in  cycles 
per  second  can  then  be  calculated  by  the  equation 

(f,  )  .  =  -3-  ^  (4_]) 

where  EI  is  the  flexural  rigidity  of  the  beam  cross  section, 
m.  is  its  mass  per  unit  of  length  and  L  is  the  length  of 
the  center  span. 

4 . 2   Solution  Parameters 

Solution  parameters  which  are  necessary  to  assure  that 
successive  cycles  of  iteration  converge  and  the  solution  is 
stable  consist  of  the  number  of  integration  steps  for  a 
complete  solution  and  the  number  of  mass  concentrations 
which  affects  the  accuracy  of  the  results.   Table  4.2 
shows  the  convergence  of  the  acceleration  at  mid-center 
span  due  to  the  variations  in  the  number  of  integration 
steps.   These  accelerations  are  evaluated  at  the  different 
positions  of  the  front  axle  of  the  vehicle  on  the  bridge. 
The  bridge  model  used  is  a  64  ft.-  80  ft.-  64  ft.  span,  I- 
beam  girder  type  v;ith  a  bridge  damping  factor  of  0.02. 
The  damping  factor  is  the  ratio  of  the  viscous  damping  of 
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Table  4.2   Comparison  of  AccGlerations  at  Mid-Center 

Span  Obtained  by  Using  Different  Numbers  of 
Integration  Steps 

64-80-64  ft.  Continuous  Span    36WF170  Girders 
n  =  4    m=4    c/Cj- =  .  02    HS  20-44  Loading 
Smooth  Vehicle    60  m.p.h.    Level  Surface 

2 
Accelerations  at  Mid-Center  Span  (in/sec  ) 

N  =  600    N  =  1000    N  =  1600    N  =  2000    N  =  2200 

-1.483  -1.479  -1.476 

-0.813  -0.830  -0.829 

11.347  11.343  11.338 

-9.237  -9.230  -9.205 

-1.009  -0.980  -1.036 

-2.986  -3.017  -2.993 

-8.083  -8.122  -8.041 

5.517  5.548  5.473 

-0.227  -0.305  -0.297 

-2.552  -2.490  -2.457 

-0.017  -0.044  -0.021 

-1.499  -1.730  -1.661 


0.1 

-1. 

.616 

-1. 

,510 

0.2 

-0. 

,591 

-0. 

,716 

0.3 

11. 

,543 

11. 

,375 

0.4 

-8. 

.525 

-9, 

.300 

0.5 

-1. 

,581 

-0. 

.831 

0.6 

-3. 

,499 

-3. 

.187 

0.7 

-7. 

,970 

-8. 

,251 

0.8 

5, 

.922 

5, 

,945 

0.9 

0. 

.073 

-0, 

.459 

1.0 

-2. 

.740 

-2, 

.760 

1.1 

-0, 

.169 

0, 

.184 

1.2 

-1. 

.289 

-1, 

.636 
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bridge  to  the  critical  damping  value.   The  value  of  0.02 
is  a  practical  value  for  most  three  span  continuous  highway- 
bridges,  according  to  References  14,  16,  19,  22  and  25. 
The  vehicle  model  is  a  three  axle  vehicle  with  HS  20-44 
loading  type  and  has  a  velocity  of  60  m.p.h.   The  initial 
conditions  are  that  the  vehicle  has  no  initial  vertical 
oscillation  before  entering  the  bridge  (hereafter  referred 
to  as  smooth  vehicle)  and  the  bridge  has  a  smooth  and 
level  surface. 

Observation  of  the  calculated  accelerations  obtained 
by  using  from  600  to  2200  integration  steps  indicates  that 
2000  integration  steps  are  sufficient  for  the  stability 
of  the  solutions.   Since  the  cost  of  the  computer  solution 
increases  as  N  increases,  the  "minimum"  number  should  be 
chosen  for  economy. 

It  was  shown  in  Reference  29  that  the  results  obtained 
from  the  bridge  model  with  n  =  4  and  m  =  4  (where  n  and  m 
represent  the  number  of  divisions  of  end  spans  and  center 
span,  respectively)  were  sufficiently  accurate  and  the 
natural  periods  of  this  model  were  close  to  those  of  the 
continuous  beam.   The  cost  of  the  computer  solution  is 
increased  for  a  larger  number  of  mass  concentrations. 
Since  a  large  number  of  solutions  were  to  be  made  for 
the  complete  parameter  study,  the  bridge  model  with  n=4 
and  m  =  4  was  considered  to  be  satisfactory  and  was  used 
throughout  this  study. 
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4  .  3   Bridge  Parameters 

4.3.1   Accelerations  of  B.P.R.  Bridges 

The  characteristics  of  the  Bureau  of  Public  Roads 
bridges  presented  in  Reference  41  are  for  I-beam  type 
bridges  which  are  designed  either  for  H  15-44  loading 
or  for  HS  20-44  loading.   The  lengths  of  the  individual 
spans  are  in  the  ratio  of  4:5:4,  and  the  overall  length 
ranges  from  130  ft.  to  260  ft.   Thus  the  shortest  bridge 
has  spans  of  40  ft.  -  50  ft.  -  40  ft.  and  the  longest 
bridge  has  spans  of  80  ft.  -  100  ft.  -  80  ft.   The  bridges 
designed  for  the  H  15-44  and  HS  20-44  loadings  have  roadway 
widths  of  28  ft.  and  44  ft.,  respectively.   Since  it  was 
shown  in  Article  3.3.1  that  the  accelerations  of  simple 
span  B.P.R.  bridges  with  roadway  width  of  44  ft.  were 
greater  than  those  with  roadway  width  of  28  ft.  because 
of  the  heavier  moving  vehicle  load,  the  three  span  continu- 
ous B.P.R.  bridges  with  roadway  widths  of  44  ft.  are  the 
only  type  to  be  studied  in  this  article. 

The  cross  section  of  the  B.P.R.  bridge  is  shown  in 
Figure  4.2  and  the  single  beam  with  the  mass  concentrations 
used  to  represent  the  bridge  in  the  analysis  is  idealized 
as  shown  in  Figure  4.3.   Also  shown  in  this  figure  is  tlie 
HS  20-44  vehicle  load  which  is  assumed  to  be  smooth 
vehicle.   It  should  be  remembered  that  the  bridge  and 
vehicle  models  shown  in  Figure  4.3  are  not  the  complete 
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models  used  in  the  analysis.   Figures  2.9  and  2.11  show 
the  complete  models. 

Table  4.3  shows  the  characteristics  of  the  three 
span  continuous  B.P.R.  bridges  including  the  v/eight  per 
unit  length,  the  fundamental  natural  frequency  and  the 
corresponding  period  of  vibration  of  the  bridges.   The 
quantities  shown  are  necessary  as  input  data  to  evaluate 
the  corresponding  accelerations  at  nodes  1  to  9 .   The 
parameters  of  three-axle  vehicle  are  considered  to  be  the 
same  as  in  Article  4.2  and  the  roadway  surface  is  hori- 
zontal and  smooth.   The  bridge  is  assumed  to  have  a  damping 
factor  of  0.02  and  this  value  is  used  for  all  bridge 
models  throughout  this  chapter. 

Table  4.4  shows  the  maximum  accelerations  at  each 
node  obtained  by  using  B.P.R.  bridge  models  with  different 
span  lengths.   The  corresponding  position  of  front  axle 
as  the  maximum  acceleration  occurs  is  also  shown.   These 
results  are  plotted  graphically  and  shown  in  Figure  4.4. 
It  can  be  seen  that  these  accelerations  increase  with 
decreasing  values  of  span  length.   In  other  words,  the 
accelerations  are  larger  for  shorter  spans.   This  is  in 
agreement  with  the  results  from  simple  span  bridges 
presented  in  Article  3.3.1  except  that  the  accelerations 
of  simple  span  bridges  are  four  times  greater  than  those 
for  the  three  span  bridges.   It  should  be  also  noted  that, 
for  most  of  the  span  lengths,  the  acceleration  of  node  2 
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Table  4.4    Maximum  Accelerations  of  Three  Span 
Continuous  B.P.R.  Bridges 

HS  20-44  Loading    60  m.p.h.    Level  Surface 
n=4    m=4    a=0.8    c/c^=0.02 


Maximum 

5 

Acceleration  (in/sec  ) 

«,   J 

40-50-40 

x/Lj^ 

48-60-48 

x/L^ 

56-70-56 

x/L^ 

1 

19.00 

0.34 

-16.59 

0.33 

-15.79 

0.46 

2 

25.12 

0.34 

-17.57 

0.38 

14.82 

0.51 

3 

16.28 

0.10 

-15.27 

0.32 

-13.83 

0.45 

4 

-12.56 

1.13 

11.11 

0.48 

-10.38 

0.63 

5 

18.90 

0.40 

-15.42 

1.12 

-10.46 

0.62 

6 

.  14.53 

0.40 

13.12 

0.36 

09.30 

0.62 

7 

-14.63 

0.41 

12.38 

1.09 

-14.38 

0.33 

8 

18.03 

0.46 

14.92 

0.55 

-14.30 

0.32 

9 

16.86 

0.46 

-17.08 

0.51 

-14.12 

0.34 

Maximum 

Acceleration  (in/: 

sec^) 

»1   J 

64-80-64 

x/L^_ 

72-90-72 

x/L^ 

80-100-80 

x/L^ 

1 

13.30 

0.23 

11.10 

0.21 

-9.55 

0.44 

2 

15.61 

0.22 

9.98 

0.07 

-10.38 

0.64 

3 

-13.00 

0.32 

10.52 

0.06 

-9.48 

0.45 

4 

10.65 

0.29 

9.77 

0.29 

8.00 

0.29 

5 

11.49 

0,30 

10.36 

0.28 

12.25 

0.29 

6 

10.59 

0.30 

7.69 

0.42 

7.28 

0.31 

7 

9.96 

1.07 

-10.10 

0.31 

-9.91 

0.31 

8 

13.24 

0.59 

-8.55 

1.05 

-9.70 

0.95 

9 

-11.53 

0.33 

-9.23 

0.30 

-10.07 

0.29 
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is  greatest  and  the  accelerations  at  the  center  of  either 
end  span  or  center  span  are  greater  than  those  at  other 
points  on  that  span. 

Figure  4.5  shows  the  history  curves  of  the  accelerations 
of  three  niidspr'in  nodes.   The  accelerations  of  node  2 
increase  rapidly  when  the  vehicle  enters  the  bridge,  while 
the  accelerations  of  nodes  5  and  8  are  not  significant  at 
this  step.   The  same  is  true  for  the  other  two  nodes,  that 
is,  the  accelerations  of  node  5  and  node  8  have  their 
largest  values  when  the  vehicle  is  on  the  center  span  and 
third   span  respectively. 

4.3.2   Effect  of  Span  Ratio 

The  bridge  models  in  the  previous  articles  had  a  span 
ratio  of  0.8  wliicli  is  commonly  used  for  three  span  contii^u- 
ous  highway  bridges.   In  this  article,  the  span  ratios 
range  from  0.6  to  1.0  and  the  resulting  effects  on  accelera- 
tion are  studied.   The  bridge  model  used  as  shown  in 
Figure  4.6(a)  has  the  center  span  L  constant  at  80  ft. 
while  both  end  spans  range  from  48  ft.  to  80  ft.   The 
remaining  parameters  of  the  vehicle  and  bridge  models 
are  considered  to  be  the  same  as  before.   The  fundamental 
frequency  and  period  of  vibration  of  each  bridge  for  the 
corresponding  value  of  span  ratio  was  evaluated  from  the 
plot  shown  in  Figure  4.1.   These  frequencies  and  the  maxi- 
mum acceleration  at  each  node  as  affected  by  the  span  ratio 
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parameter  are  presented  in  Table  -1 . 5 .   Thet;e  results  are 
also  shown  graphically  in  Figure  4.7.   It  can  be  seen 
that  the  span  ratio  has  no  effect  on  the  magnitudes  of 
acceleration.   It  follows  that  the  accelerations  of  the 
end  spans  are  greater  than  that  of  the  center  span  when 
the  span  ratio  is  less  than  0.8  and  these  accelerations 
tend  to  be  equal  when  the  span  ratio  approaches  the  value 
of  1.0.   This  can  be  explained  by  the  fact  that  the  acceler- 
ation is  larger  for  the  shorter  span. 

4.3.3   Effect  of  Girder  Flexibility 

Figure  4.6(b)  shows  the  bridge  model  used  for  the  study 
of  the  accelerations  as  affected  by  using  different  values 
of  girder  stiffness.   The  bridge  has  64  ft.  -  80  ft.  -  64  ft. 
span  lengths  witli  a  damping  factor  of  0.02  and  a  level  sur- 
face.  The  remaining  parameters  are  considered  to  be  the 
same  as  before.   Five  different  girder  sections  are  used 
to  evaluate  the  accelerations  at  each  node.   Tlio  36VJF170 
is  supposed  to  be  the  proper  section  for  the  design  of 
this  bridge.   The  section  is  reduced  and  four  smaller 
sections  are  used  to  study  the  effect  of  girder  flexibility. 
The  smallest  section  used  is  the  33WF118  v/hich  has  a 
moment  of  inertia  nearly  half  that  of  the  36WF170.   The 
properties  of  each  section  v/ere  obtained  from  the  informa- 
tion given  in  Reference  4  2  and  shown  in  Table  4.6. 
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Maximum  accelerations  at  nodes  are  obtained  by  using 
these  five  different  sections  while  keeping  all  other 
parameters  constant.   These  results  are  shown  in  Table 
4.6  and  are  plotted  graphically  as  a  function  of  girder 
moment  of  inertia  in  Figure  4.8.   It  can  be  seen  from 
these  results  that  the  accelerations  of  the  bridges  analyzed 
are  not  materially  affected  by  the  increased  girder 
flexibilities . 

4.4   Vehicle  Parameters 

4.4.1   Effect  of  Vehicle  Loads 

Three  different  types  of  commercial  trucks  have  been 
considered  as  the  vehicle  model  used  on  the  64  ft.  -  80  ft.  - 
64  ft.  B.P.R.  bridge  at  a  speed  of  60  m.p.h.   These  trucks 
are  US  20-44,  HS  15-44  and  H  20-44.   The  total  load  on  each 
axle  of  each  type  was  obtained  from  the  information  given 
in  Reference  43.   Since  the  bridge  was  idealized  as  a 
single  beam  in  the  analysis,  these  vehicle  loads  are  also 
idealized  as  a  single  line  and  the  resulting  load  on  each 
axle  is  shown  in  Figure  4.9(a).   The  total  vehicle  loads 
for  US  20-44,  HS  15-44  and  H  20-44  are  72,  54  and  40  kips 
respectively . 

Each  truck  was  simulated  to  traverse  the  bridge  model 
and  the  m.aximum  accelerations  at  nodes  were  obtained  as 
shown  in  Table  4.7.   As  might  be  expected,  the  accelerations 
of  the  bridge  are  greatest  for  the  heaviest  vehicle  load. 
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Table  4,7    Maximum  Accelerations  of  Three  Span  Con- 
tinuous Bridges  Obtained  by  Using  Different 
Types  of  Vehicle 

64-80-64  B.P.R.  Bridge    n  ==  4    m=4    c/c  =,02 
Level  Surface    Smooth  Vehicle    60  m  .  p  .  li . 


Maximum 

Accelerat 

ion  (in/sec^] 

1 

Node 

HS  20-44 

x/L^^ 

HS  15-44 

x/Lb 

H  20-44 

x/L^ 

1 

13.30 

0.23 

9.75 

0.23 

10.60 

0,17 

2 

15.61 

0.22 

11,37 

0.22 

12.83 

0.16 

3 

-13.00 

0,32 

-9.11 

0.32 

11.12 

0.16 

4 

10.65 

0.29 

8,3  6 

0.29 

7.67 

0,2  3 

5 

11.49 

0.30 

9.27 

0.30 

9,50 

0,24 

6 

10.59 

0.30 

8,39 

0.30 

8.13 

0.2  4 

7 

9.96 

1.07 

7.29 

1.07 

7.62 

0.69 

8 

13.24 

0.59 

-9.13 

0.92 

-8.11 

0,26 

9 

-11.53 

0.33 

-8.52 

0.33 

-9.62 

0.27 
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It  can  be  seen  that  the  accelerations  for  IIS  15-44  and 
H  20-44  are  close  although  HS  15-44  has  the  larger  total 
load.   It  can  be  explained  that  HS  15-44  has  a  larger 
number  of  axles  than  H  20-44  and  the  number  of  axles  has 
an  effect  on  acceleration.   This  will  be  discussed  further 
in  Article  4.4.2.   Study  results  reported  in  this  article 
indicate  that  the  heavier  vehicle  load  will  cause  greater 
accelerations  of  the  bridge. 

4.4.2   Effect  of  Number  of  A.xles  and  Axle  Spacings 

Figure  4.9(b)  shows  the  typical  vehicles  with  different 
number  of  axles.   The  parameters  of  typical  vehicles  were 
given  in  Table  4.1.   The  bridge  model  and  the  remaining 
parameters  are  the  same  as  before.   Although  the  values 
of  the  axle  load  for  tv/o-axle  and  single  axle  vehicle 
are  not  practical,  the  results  of  the  accelerations  of 
the  bridge  obtained  will  be  useful  for  the  further  studies 
of  design  criteria.   These  results  can  also  be  compared 
with  those  for  simple  span  bridges,  for  which  the  vehicle 
model  was  a  single  axle  load. 

The  maximum  accelerations  at  various  nodes  due  to 
different  moving  typical  vehicles  are  compared  in  Table 
4.8  and  shown  graphically  in  Figure  4.10.   It  can  be  seen 
that  there  is  not  much  difference  between  the  acceleration 
magnitudes  for  three  axle  and  two  axle  vehicles.   The  front 
axle  load  of  8  kips  does  not  have  much  influence  on 
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Node 


Table  4.8    Maximum  Accelerations  of  Three  Span 

Continuous  Bridges  Sugjected  to  Vehicles 
with  Different  Numbers  of  Axles 

64-80-64    B.P.R.  Bridge    n=4    m=4    c/c  =.02 
Level  Surface    Smooth  Vehicle    60  m.p.h. 


Maximum  Acceleration  (in/sec^) 


3-axle   ^/Li        2-axle  ^/^u  1-axle   ^/Lu 


11.56   0.45        -12.18   0.14 


2 

-11.56 

0.27 

13.14 

0.02 

3 

10.92 

0.07 

12.12 

0.02 

4 

10,50 

0.30 

10.93 

0.24 

5 

14.33 

0.31 

13.89 

0.25 

6 

9.01 

0.30 

8.69 

0.24 

7 

-12.21 

0.27 

-12.93 

0.21 

8 

-13.26 

0.42 

-13.89 

0.21 

9 

13.16 

0.22 

14.14 

0.16 

21. 

.12 

0. 

.03 

26, 

.16 

0. 

,02 

24, 

.24 

0, 

.02 

14, 

.03 

0, 

.24 

-17, 

.99 

0, 

.43 

-14, 

.66 

0, 

.04 

18, 

.34 

0, 

.54 

21, 

.15 

0, 

,54 

-21, 

.51 

0, 

.50 
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acceleration.   As  might  be  expected,  the  magnitudes  of 
accelerations  under  single  axle  loads  are  almost  double 
the  values  of  two  and  three  axle  vehicles. 

Table  4.9  and  Figure  4.11  show  the  effect  of  vehicle 
axle  spacing  on  bridge  accelerations  .  The  vehicle 
model  is  the  HS  20-44  loading  with  trailer  axle  spacing 
varying  from  20  ft.  to  35  ft.  the  tractor  axle  spacing 
is  kept  constant  at  12  ft.  and  the  vehicle  speed  is  60 
m.p.h.   All  remaining  vehicle  parameters  are  the  same  as 
in  Article  4.3.3.   The  bridge  model  is  the  64  ft.  -  80  ft.  - 
64  ft.  B.P.R.  bridge  with  a  damping  factor  of  0.02  and 
a  level  surface. 

From  the  results,  the  magnitudes  of  the  maximum 
accelerations  are  not  affected  by  varied  axle  spacing 
although  these  values  are  slightly  higher  when  the 
axle  spacing  is  between  30  ft.  and  35  ft.   It  should 
be  noted  that  these  results  will  be  changed  if  different 
span  lengths  of  the  bridge  model  are  used.   In  order 
to  consider  the  results  in  a  non-dimensional  form,  the 
ratio  of  the  axle  spacing  of  the  vehicle  to  the  length 
of  the  center  span  is  introduced.   It  can  be  seen  that 
for  the  three  span  continuous  bridge,  the  accelerations 
slightly  increase  when  the  axle  spacing  ratio  ranges  between 
0.37  to  0.43. 

4.4.3   Effect  of  Vehicle  Speed 

Table  4.10  and  Figure  4.12  show  the  maximum  accelerations 
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of  the  three  span  continuous  bridge  affected  by  the  vehicle 
speeds  between  30  m.p.h.  to  7  0  m.p.h.  The  bridge  model 
is  the  64  ft.  -  80  ft.  -  64  ft.  B.P.R.  bridge  and  has  the 
remaining  parameters  the  same  as  in  Article  4.4.2.  The 
vehicle  model  has  HS  20-44  loading.  As  might  be  expected, 
the  accelerations  of  the  bridge  increase  as  the  vehicle 
speed  increases. 

4.4.4   Effect  of  Frequency  Ratio 

Table  4.11  and  Figure  4.13  show  the  maximum  bridge 
accelerations  obtained  by  using  the  values  of  frequency 
ratio  f  /f,  in  the  range  between  0.5  and  2.0.   The  bridge 
model  used  is  64  ft.  -  80  ft.  -  64  ft.  B.P.R.  bridge  as 
used  in  Article  4.4.2  and  has  a  fundamental  natural  fre- 
quency of  3 . 0  c.p.s.   The  vehicle  has  HS  20-44  loading 
and  is  assumed  to  oscillate  on  its  tires  only.   The 
frequency  of  vehicle  f   used  is  the  tire  frequency  of  the 
rear  axles  and  the  corresponding  coefficients  of  friction 
of  the  suspension  spring  have  an  infinite  value. 

It  can  be  seen  from  the  results  shown  that  the 

frequency  ratio  has  no  appreciable  effect  on  the  bridge 

accelerations  and  the  magnitudes  of  the  accelerations  in 

these  curves  arc  about  the  same  for  all  values  of  frequency 

ratio.   In  particular,  the  accelerations  for  f  /f,  =^  1.0 
'^  V   b 

are  no  larger  than  those  for  the  other  frequency  ratios. 
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Figure  4.13    Effect  of  Frequency  Ratio  on  Acceleration 
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'1 .  5   Initial  Conditions  of  VGhicle  and  Bridge 

4.5.1   Effect  of  Initially  Oscillating  Vehicle 

It  is  seldom  that  a  vehicle  enters  the  bridge  with  its 
suspended  mass  in  equilibrium  for  vertical  motion.   It  is 
therefore  necessary  to  investigate  the  effect  of  initially 
oscillating  vehicles  on  bridge  motion.   The  initial  oscilla- 
tion of  the  vehicle  may  arise  either  from  the  ever  present 
irregularities  of  the  approach  pavement  or  from  a  sharp 
discontinuity  between  the  approach  pavement  and  the  entrance 
of  the  bridge. 

The  method  of  analysis  assumes  that,  while  the  vehicle 
is  on  the  approach  pavement,  the  vertical  oscillation  of 
each  axle  is  of  the  simple  harmonic  type.   The  initial 


th 
i'  --—--'■  -"-  -^ 

and    tlio   pavement    is   expressed   as 


value  of  the  interacting  force.  P.,  between  the  i^"  axle 


P.  =  (1  +C  .  cos  0.)P^.  .  (4-2) 

1  1      1   st , 1 


where  C.  is  the  amplitude  of  the  initial  force  variation, 

e.  is  the  phase  angle  betv;een  the  time  at  which  the  force 

attains  its  maximum  value  and  tlie  time  at  which  the  front 

axle  enters  the  bridge,  and  P  .  .  is  the  static  reaction 

-'         s  t ,  1 

on  the  axle.   For  a  two-axle  loading,  the  phase  angles 

0^  and  0-  will  in  general  be  different.   The  quantity 

AO  =  0,  -  0„  will  be  referred  to  as  the  phase  angle  difference 

of  the  vehicle. 
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The  values  of  C.  and  0.  depend  upon  such  variables 
as  the  dimensions  and  location  of  the  irregularities  and 
discontinuity,  and  upon  the  speed  of  the  vehicle.   The 
initial  value  of  the  frictional  force  depends,  in  addition, 
on  the  past  history  of  the  deformation  of  the  suspension 
spring,  and  may  have  any  value  bctv/een  F'  and  -F '  . 

It  was  shown  in  Reference  30  that  the  height  of 
disco  itinuity  required  to  produce  a  variation  of  interacting 
force  of  0.15  P  ^,    or  C.  =  0.15,  is  about  0.12  in.   The 
value  of  C.  =  0.15  is  considered  to  represent  the  vehicle 
with  an  initial  oscillation  due  to  the  roughness  of  the 
approach  pavement.   In  addition,  the  value  of  C.  =0.50 
is  considered  to  define  the  effect  of  a  large  discontinuity 
at  the  abutment,  or  of  a  large  irregularity  on  the  approach 
pavement  located  close  to  the  entrance  of  the  bridge. 
I'he  value  of  phase  angle,  0.,  was  considered  to  vary  between 
0°  and  360°  and  the  initial  value  of  interleaf  friction 
was  assigned  the  values  of  F ' ,  zero  and  -F ' , 

4.5.1.1   Effect  of  Amplitude  of  Initial  Oscillation 

Maximum  accelerations  of  the  bridge  as  affected  by 
the  different  amplitudes  of  initial  oscillation  of  vehicle 
are  studied.   The  bridge  model  is  considered  to  be  the 
same  as  in  Article  4.4.3.   The  vehicle  model  has  an 
HS  20-44  loading  with  the  initial  oscillation  at  the 
approach  pavement  and  travels  at  the  constant  speed  of 


17^ 


60  m.p.h.   The  amplitudes  of  initial  oscillation  of 
vehicle,  C.,  are  assigned  to  have  the  values  of  0.0,  0.15, 
0.30  and  0.50.   These  four  values  represent  the  smoothly 
moving,  the  oscillation  due  to  the  surface  roughness,  the 
oscillation  due  to  a  larger  amplitude  of  roughness  and 
tlie  oscillation  due  to  a  large  discontinuity,  all  at 
the  approach  pavement,  respectively.   The  values  of  C- 
are  considered  to  be  the  same  for  each  axle.   The  initial 
values  of  the  frictional  force  and  phase  angle  for  each 
axle  are  all  taken  equal  to  zero. 

The  maximum  accelerations  at  nodes  evaluated  by 
using  the  different  values  of  C.  are  compared  in  Table  4.12 
and  also  are  shown  graphically  in  Figure  4.14.   It  can 
be  seen  that  the  accelerations  increase  for  the  larger 
amplitude  of  initial  oscillation.   The  magnitude  of 
maximum  acceleration  of  all  nodes  increases  as  much  as 
1.5  times  the  original  value  for  a  smooth  vehicle.   It 
can  be  noted  that  the  acceleration  of  the  mid-center 
span  is  the  largest  value  when  the  initial  oscillation  of 
vehicle  is  taken  into  account. 

4.5.1.2   Effect  of  Initial  Phase  Angle  Difference 

It  is  seldom  that  the  initial  oscillating  three  axle 
vehicle  which  passes  over  the  irregularities  at  the  approach 
will  have  the  maximum  initial  forces  for  all  axles.   In 
other  words,  all  axles  would  not  have  the  same  value  of 
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initial  phase  angle  of  zero  to  obtain  the  maximum  force 
referred  to  Equation  (4-2) .   It  can  be  explained  that  the 
vehicle  has  two  axle  spacings,  the  spacing  between  the 
first  and  second  axles,  and  the  spacing  between  the  second 
and  third  axles.   It  will  take  quite  a  time  between  first 
and  the  following  axles  to  pass  the  same  irregularities 
at  the  approach.   Assuming  that  the  force  of  the  first 
axle  attains  the  maximum  value  at  the  approach  or  0,  =  0° , 
the  forces  of  second  and  third  axles  will  or  will  not 
attain  their  maximum  values  at  the  approach  or  0_  and  0^ 
can  vary  from  0°  to  3  60°. 

The  effect  of  initial  phase  angle  difference  of  the 
axles  on  accelerations  was  investigated  and  the  results 
are  shov/n  in  Table  4.13  and  Figure  4.15.   The  bridge 
model  is  considered  to  be  the  same  as  in  Article  4.5.1.1. 
The  vehicle  model  has  HS  20-44  loading  with  15%  initial 
oscillation  for  all  axles.   The  value  of  initial  coefficient 
of  friction  is  taken  equal  to  zero  for  each  axle.   The 
initial  phase  angles  for  first  and  second  axles  are  kept 
constant  at  0°  while  the  initial  phase  angle  of  third 
axle  is  varied  from  0°  to  180°.   The  reason  for  keeping 
the  phase  angle  of  the  first  two  axles  at  zero  while 
studying  the  accelerations  due  to  the  difference  between 
the  second  and  third  is  that  the  first  spacing  is  somewhat 
shorter  than  the  second  spacing,  so  the  time  required  for 
the  force  to  change  its  value  is  less. 
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Table  4.13    Maximuin  Accelerations  of  Three  Span 
Continuous  Bridges  Subjected  to  15o 
Initially  Oscillating  Vehicle  with 
V^lried  Phase  Angle  Differences 

64-80-64  B.P.R.  Bridge    n=4    m  =-  4    c/c  =.02 
Level  Surface    HS  20-44  Loading    60  m.p.ii. 


Oscillating    Vehicle         C.  =--  0 . 

15 

"i  " 

0 

^1="2= 

0°          AO    =    O^-O^ 

Maximum  Acceleration 

(in/ 

sec^ ) 

Node 

AO  =-  0° 

x/L^ 

A0  =  30°       x/Lj^ 

AG  = 

60° 

x/Lj^ 

AG  = 

90° 

X/Lj^ 

1 

17.08 

0.23 

17.3  6         0.2  3 

17, 

.33 

0.2  3 

17, 

.84 

0.23 

2 

18.99 

0.22 

19.18         0.22 

19. 

.84 

0.22 

19, 

.18 

0.2  2 

3 

-13.20 

0.3  2 

-14.73          0.32 

-12. 

.05 

0.32 

-13, 

.18 

0.17 

4 

-11.12 
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-12.62          0.40 

-10, 

.51 

0.4  0 

11. 

,15 
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-17. 

,96 
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.17 
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-15.15         0.39 
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,92 
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7 
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-14, 
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-16, 
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8 
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-19, 
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-21, 

.93 
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9 
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-18, 

.09 
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-19, 

.19 

0.33 

AG  =  120° 

x/Lj^ 

AO  = 

150° 

x/Lj^ 

A  0  ^ 

180° 

x/L^ 

1 

17.18 

0.23 

17. 

,80 

0.23 

18, 

.39 
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2 

20.12 

0.22 

22. 
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0.2  2 
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3 
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0.  32 
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,06 
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4 
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15, 

.36 
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From  the  results  shown  in  Table  4.13  and  Figure  4.15, 
the  accelerations  are  not  much  affected  by  the  initial 
phase  angle  difference  ranging  from  0°  to  180°  .   The 
accelerations  of  midspan  nodes  are  still  the  dominant 
values  for  whole  range.   It  should  be  noted  that  the  effect 
of  the  phase  angle  difference  ranging  from  180°  to  360° 
is  not  shown  because  the  values  of  initial  forces  were 
similar  to  those  in  0°  to  180°  range. 

The  effects  of  initial  phase  angle  difference  were 
also  investigated  for  a  vehicle  model  having  50%  initial 
oscillation  due  to  a  discontinuity  at  the  approach.   All 
remaining  parameters  are  the  same  as  before.   The 
accelerations  obtained  are  shown  in  Table  4.14  and  Figure 
4.15.   Although  these  accelerations  were  markedly 
increased  for  the  angle  differences  equal  to  90° 
and  180° ,  it  can  be  seen  that  the  initial  phase  angle 
difference  has  no  appreciable  effect  on  acceleration. 

4.5.1.3   Effect  of  Coefficient  of  Interleaf  Friction 

In  the  previous  studies  of  this  chapter,  the  coefficient 
of  friction  of  the  suspension  springs  for  the  smoothly 
moving  vehicle  was  assumed  to  have  the  value  of  0.15  or 
)i  =  0.15,  except  when  the  effect  of  frequency  ratio  was 
studied,  for  which  \i    was  taken  as  infinite.   In  this  study, 
the  vehicle  has  an  initial  oscillation  and  initial  coef- 
ficient of  friction  of  suspension  spring  (y.).   The 
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Table  4.14   Maximum  Accelerations  of  Three  Span 
Continuous  Bridges  Subjected  to  50% 
Initially  Oscillating  Vehicle  with 
Varied  Phase  Angle  Differences 

64-80-64  B.P.R.  Bridge    n=4    m=4    c/c  =.02 
Level  Surface    HS  20-44  Loading    60  m.p.fi. 
Oscillating  Vehicle    C.=0.50    y-=0 
0j^  =  ©2  ==  0°    A0  =  02"®3"^ 


Node 

Maximum 

Acceleration  (in/sec^) 

1 

A0  = 

0° 

x/L^ 

AG  =  30° 

x/L^ 
b 

A0  =  60° 

x/Lb 

A0  = 

90° 

x/Lj^ 

1 

18. 

,08 

0.28 

15.97 

0.28 

17.73 

0.23 

19. 

,12 

0.28 

2 

19. 

,04 

0.28 

-20.29 

0.38 

19.86 

0.28 

21, 

.63 

0.28 

3 

17. 

.05 

0.29 

-16.01 

0.38 

17.09 

0.29 

-16. 

,20 

0.32 

4 

17. 

.87 

0.36 

-18.36 

0.45 

17.74 

0.36 

18. 

.94 

0.36 

5 

-20. 

.82 

0.46 

-24.41 

0.46 

-20.99 

0.46 

26, 

,25 

0.35 

6 

-17. 

,74 

0.46 

-18.77 

0.46 

-15.35 

0.46 

16. 

,58 

0.35 

7 

-15. 

,50 

0.39 

-12.21 

0.38 

-16.47 

0.38 

19, 

.34 

0.42 

8 

-20. 

.34 

0.39 

-17.07 

0.39 

-20.43 

0.39 

28, 

,57 

0.42 

9 

-17. 

.87 

0.40 

-15.07 

0.40 

-18.89 

0.40 

-22, 

.04 

0.47 



AO  =  120° 

x/Lb 

A9  =  150° 

x/Lb 

A0  =  180° 

x/L^ 

- 

1 

-14. 

.43 

0.18 

-14.43 

0.18 

-15.65 

1.05 

2 

-15, 

.18 

0.25 

-16.98 

0.39 

-19.31 

1.06 

3 

12, 

.72 

0.29 

-13.03 

0.40 

-14.48 

0.32 

4 

13, 

,53 

0.36 

14.88 

0.88 

13.66 

0.35 

5 

-18, 

.00 

0.47 

16.57 

0.35 

21.27 

0.35 

6 

-14, 

.31 

0.46 

13.07 

1.01 

16.03 

0.35 

7 

-9, 

.98 

0.38 

15.06 

0.50 

-16.18 

0.46 

8 

-13, 

.54 

0.39 

-23.58 

0.55 

23.04 

0.71 

9 

-13, 

,08 

0.40 

-19.63 

0.55 

19.92 

0.51 
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numerical  solutions  of  the  accelerations  of  bridge  are 
obtained  for  three  different  values  of  the  coefficient 
of  interleaf  friction:  y  =  ~,  p  =  0  and  y  =  0.15.  It 
is  assumed  that  the  initial  coefficients  of  friction  of 
suspension  springs  are  infinite  for  p  -  <«  and  equal  to 
zero  for  p  =  0  and  p  =  0.15. 

As  previously  explained,  for  p  ~  ■»  the  suspension 
spring  remains  idle  and  the  vehicle  oscillates  on  its 
tires  only.   for  p  =  0,  the  suspension  spring  acts  in 
series  with  the  tire  spring.   Finally,  for  p  =  0.15  the 
frictional  force  exists  and  the  energy  dissipates  in 
tlie  suspension  spring.   The  vehicle  model  is  assumed  to 
hcivc    15  percent  of  initial  interacting  force  and  initial 
phdse  angle  equal  to  zero  for  each  axle.   The  vehicle  has 
US    20-44  loading  and  travels  at  60  m.p.h.   The  bridge 
model  and  the  remaining  parameters  are  the  same  as  in 
Article  4.5.1.2. 

The  results  of  investigations  using  three  different 
values  of  ii  are  obtained  and  shown  in  Table  4.15.   It 
can  be  seen  that  the  magnitudes  of  accelerations  are  not 
much  dif f eient  for  the  three  different  values  of  p .   The 
accelerations  for  p  =  0.15  are  slightly  less  than  those 
for  p  "  ""  because  of  the  energy  dissipated  in  the  suspension 
spring  and  the  amplitude  of  interacting  force  for  p  -    0.15 
was  reduced,  as  described  in  Reference  30.   The  accelerations 
for  p  =  0  are  supposed  to  be  higher  than  those  for  p  =  0.15 
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Table  4.15    Maximum  Accelerations  of  Three  Span 
Continuous  Bridges  Subjected  to 
Initially  Oscillating  Vehicle 
with  Different  Coefficients  of 
Friction  of  Suspension  Spring 

64-80-64  B.P.R.  Bridge    n=4    m=4    c/c  =.02 
Level  Surface    HS  20-44  Loading    60  m.p.n. 
Initial  Oscillating  Vehicle    C-  =0.15 
0j_  =  0  Vj^=     °"    for  u  =  °°  y^  =  0  for  p  =  0  and  ij  =  0.15 


Maximum   Ace 

deration 

(in/sec^; 

) 

Node 

\>    - 

-    °° 

x/Lj^ 

p  =  0.15 

x/Lj^ 

V  =  0 

x/Lb 

1 

16. 

,07 

0.89 

18.08 

0.23 

-15.06 

0.56 

2 

-20. 

.94 

0.94 

19.04 

0.22 

-16.09 

0.25 

3 

13. 

,51 

0.90 

17.05 

0.32 

13.34 

0.51 

■    4 

15. 

,62 

0.29 

17.87 

0.40 

-14.34 

0.50 

5 

23. 

.19 

0.30 

-20.82 

0.40 

-18.19 

0.63 

6 

19. 

,78 

0.30 

-17.74 

0.39 

-13.46 

0.63 

7 

14. 

.24 

0.67 

-15.50 

1.07 

-13.43 

0.41 

8 

19, 

,52 

0.66 

-20. 34 

0.92 

-14.06 

0.41 

9 

17, 

.06 

0.75 

-17.87 

0.33 

-13.16 

0.40 
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but  they  are  the  least.   It  can  be  explained  that  the 
vehicle  frequency  is  different  for  each  case. 

4.5.2   Effect  of  Bridge  Surface  Roughness 

In  the  preceding  study  of  three  span  continuous 
bridge,  the  surface  of  the  bridge  model  is  assumed  to 
be  smooth  and  level.   Generally,  this  condition  is  seldom 
the  case  for  most  highway  bridges.   The  surface  roughness 
ol  the  bridge  specified  in  this  study  includes  the 
initial  camber,  grade,  vertical  curve  and  roadway 
unovonness.   It  is  designated  as  the  term  d  .  described 
i)i  Article  2.3.2.   For  a  given  bridge,  this  quantity 
may  be  presumed  to  be  known  and  its  amplitude  practically 
ranges  from  0.0  to  1.0  inch. 

Tlio  surface  roughness  can  be  of  any  shape  due  to 
factors  described  above.   For  this  study,  the  profile 
of  the  surface  roughness  is  idealized  as  a  series  of 
continuous  half  sine  waves  along  the  span  length.   The 
amplitude  of  the  roughness  is  always  equal  to  zero  at  the 
support  and  is  maximum,  at  the  middle  of  the  wave.   For 
the  three  span  continuous  bridge,  each  span  has  an  equal 
number  of  half  sine  waves.   Figure  4.16  illustrates  the 
examples  of  the  profiles  of  bridge  models  having  one, 
two  and  four  half  sine  waves.   The  first  wave  of  each 
span  has  a  positive  or  upward  amplitude.   For  a  given 
bridge  model,  all  peak  amplitudes  of  half  sine  waves  are 
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assumed  to  be  the  same.   In  this  study,  two  parameters 
representing  the  surface  roughness  arc  investigated: 
the  number  of  half  sine  waves  and  the  amplitude  of  rough- 
ness. 

The  effect  of  different  numbers  of  half  sine  waves 
on  bridge  accelerations  are  investigated  while  keeping 
the  peak  amplitude  of  half  sine  waves  constant.   The 
bridge  model  is  a  64  ft.  -  80  ft.  -  64  ft.  damped  B.P.R. 
bridge.   The  vehicle  model  has  HS  20-44  loading,  smoothly 
moving  type  v/ith  the  speed  of  6  0  m.p.h.   The  peak  amplitudes 
of  roughness  used  are  kept  constant  at  0.5  in.   All 
remaining  parameters  are  considered  to  be  the  same  as 
before. 

The  maximum  accelerations  at  the  nodes  are  evaluated 
by  using  from  zero  to  12  half  sine  waves.   The  results 
are  listed  in  Table  4,16  and  shown  graphically  in  Figure 
4.17.   It  can  be  seen  that  the  accelerations  for  the  number 
of  half  sine  v/aves  ranging  from  0  to  2  are  nearly  constant. 
These  accelerations  increase  sharply  and  attain  the 
maximum  for  four  half  sine  waves  and  gradually  decrease 
for  more  than  four  half  sine  waves.   The  magnitudes  of 
accelerations  for  four  half  sine  waves  are  very  much 
greater  than  those  for  the  level  surface  condition.   This 
can  be  explained  in  Figures  4.18-1  and  4.18-2  which  shows 
the  histories  of  interacting  forces  of  the  second  axle 
and  accelerations  at  mid-center  span  respectively.   In 
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O  NODE 


NODE 
+  NODE 
X  NODE 
O  NODE 
^  NODE 
«  NODE 
)K  NODE  8 
X   NODE  9 


2       4       6       8       10 

NUMBER  OF  HALF  SINE  WAVES  ON  EACH  SPAN 


64-80-64  BPR  BRIDGE   C/CR=0.02   N=H=4   ROUGH     SURFACE 
HS  20-44    SMOOTH  VEHICLE   60  MPH 


Figure  4.17   Effect  of  Surface  Roughness  on  Acceleration 
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each  figure,  the  results  obtained  by  using  the  level 
surface  condition  and  four  half  sine  waves  are  compared. 

It  can  be  seen  that  the  period  of  the  oscillation  of 
interacting  force  for  the  value  of  four  half  sine  waves 
is  coincident  with  the  second  natural  period  of  the  bridge 
and  resonance  between  the  bridge  and  the  vehicle  occurs. 
The  results  show  that  the  acceleration  of  the  bridge 
increases  indefinitely  although  the  vehicle  is  leaving 
the  bridge.   On  the  other  hand,  the  period  of  the  interact- 
ing force  for  the  level  surface  bridge  is  less  than  the 
period  of  the  fundamental  frequency  of  the  bridge  and 
the  period  of  the  corresponding  acceleration  at  mid-center 
span  corresponds  to  the  lowest  three  natural  periods  of 
vibration  of  the  bridge.   Resonance  does  not  occur  in  this 
case . 

Table  4.17  and  Figure  4.19  show  the  effects  of  the 
roughness  iimpiitude  on  acceleration.   The  numerical 
acceleration  results  are  obtained  by  using  amplitudes  of 
roughness  which  range  from  0.0  to  1.0  inch.   The  number 
of  half  sine  waves  is  kept  constant  and  the  remaining 
parameters  are  the  same  as  before.   Two  different  numbers 
of  half  sine  waves  are  selected  and  classified  into 
tv;o  cases.   In  case  A,  the  two  half  sine  waves  are  selected, 
which  do  not  affect  the  acceleration  of  the  bridge,  while 
in  case  B,  four  half  sine  waves  are  used,  which  have  most 
effect. 
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2  HALF  SINE  WAVES 

^   HALF  SINE  WAVES 

O  NODE  1 

A  NODE  2 

+  NODE  3 

X  NODE  4  / 

♦  NODE  5  / 

^  NODE  6  / 

«  NODE  7  / 

)K  NODE  8  / 

X   NODE  9  ^ 
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SURFACE 


Figure  4-19   Effect  of  Amplitude  of  Roughness  on  Acceleration 
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It  can  be  seen  from  the  results  that  accelerations 
for  the  two  half  sine  waves  are  almost  constant  reyardless 
of  the  ^lmplitude  of  roughness  while  the  accelerations 
for  the  four  half  sine  waves  gradually  increase  for  larger 
values  of  the  amplitude.   This  is  in  good  agreement 
v/ith  the  results  presented  in  Article  3.5  for  simple 
span  bridge.   Comparing  the  tv/o  parameters  which  represent 
the  surface  roughness,  one  can  see  that  the  effect  of 
roughness  amplitude  on  acceleration  depends  upon  the 
number  of  half  sine  waves.   The  amplitude  affects  the 
acceleration  only  if  the  number  of  half  sine  waves  does. 

It  should  be  remembered  that  these  results  were 
obtained  for  only  this  kind  of  bridge  and  vehicle  models 
vv'ith  given  velocity.   However,  it  can  bo  seen  at  this 
stage  that  the  surface  roughness  of  the  three  span  continuous 
bridge  is  the  most  important  parameter  that  affects  the 
acceleration. 
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CHAPTER  V 
ACCELERATION  STUDIES  OF  TWO  SPAN  CONTINUOUS  HIGHWAY  BRIDGES 

5 . 1   General 

The  accelerations  of  two  span  continuous  highway 
bridges  as  affected  by  the  important  parameters  are 
reported  in  this  chapter.   The  bridge  was  idealized  as 
a  two  span  continuous  flexible  beam  subjected  to  the  moving 
three-axle  vehicle  with  constant  velocity.   The  method 
of  analysis  and  the  description  of  the  computer  program 
used  have  been  presented  in  Articles  2.3.3  and  2.3.4 
respectively.   The  numerical  results  of  accelerations 
investigated  at  the  concentrated  masses  of  the  bridge 
are  in  dimensional  form  so  that  it  will  be  convenient 
to  compare  these  results  with  the  suggested  acceleration 
criteria  for  human  response. 

As  in  the  previous  chapter,  the  parameters  affecting 
the  acceleration  can  be  classified  into  the  following 
four  groups : 

a)  Solution  parameters. 

b)  Bridge  parameters  including  the  span  length  and 
the  girder  flexibility. 

c)  Vehicle  parameters  including  the  axle  spacing 
and  the  vehicle  speed. 
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d)  Initial  conditions  of  vehicle  and  bridge  including 
the  initial  oscillation  of  the  vehicle  and  the 
surface  roughness  of  the  bridge. 

Throughout  this  chapter,  the  bridge  is  considered 
to  be  an  I-beain  type  composed  of  steel  girders  and  a 
reinforced  concrete  deck  with  non-composite  section. 
It  has  a  prismatic  cross  section  and  equal  span  lengths. 
Damping  is  considered  in  both  the  bridge  and  vehicle. 
The  bridge  and  vehicle  models  are  referred  to  in  Figures 
2.15  and  2.11  respectively.   The  computed  accelerations 
of  the  bridge  occur  at  the  specified  concentrated  masses 
and  the  maximum  acceleration  value  together  with  the 
corresponding  front  axle  position  at  that  time  are  considered. 

The  natural  frequencies  of  two  span  continuous  bridges 
can  bo  evaluated  by  substituting  the  values  of  A.  obtained 
from  Figure  4.1  into  Eq .   (4-1).   Most  of  the  vehicle 
parameters  were  obtained  from  the  information  for  typical 
vehicles  presented  in  Table  4.1. 

5.2   Solution  Parameters 

Table  5.1  shows  the  convergence  of  the  acceleration 
at  mid-right  span  due  to  the  various  values  of  the  number 
of  integration  steps.   The  acceleration  is  evaluated 
at  the  different  positions  of  the  front  axle  of  the 
vehicle  on  the  bridge.   The  bridge  model  used  has  two  80  ft. 
spans,  with  a  damping  factor  of  0.02  and  a  level  bridge 
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surface.   The  vehicle  model  is  a  three  axle  smooth  vehicle 
v/ith  US  20-44  loading  travelling  at  60  m.p.h.   By  con- 
sidering the  accelerations  obtained  by  using  from  600 
to  2200  integration  steps,  one  finds  that  the  2000  inte- 
gration steps  is  sufficient  for  the  stability  of  the 
solutions. 

The  other  solution  parameter  to  assure  the  accuracy 

of  results  is  the  number  of  mass  concentrations.   It 

(29) 
was  suggested  by  Huang      that  the  results  obtained  from 

the  bridge  m.odel  with  n  =  4  were  sufficiently  accurate 

and  the  natural  periods  of  this  model  were  close  to  those 

of  the  continuous  beam.   Therefore  the  value  of  n  =  4 

wi]l  be  used  for  the  bridge  model  throughout  this  chapter. 

5 . 3   Bridge  Parameters 

5.3.1   Effect  of  Span  Length 

Table  5.2  shows  the  characteristics  of  the  two  span 
bridge  models  with  three  different  span  lengths,  60  ft.  -  60 
ft.,  80  ft.  -  80  ft.  and  100  ft  -  100  ft.   The  cross  section 
of  th(2  bridge  model  is  shown  in  Figure  5.1(a).   The  bridge 
is  of  I-beam  type  consisting  of  six  wide-flange  steel 
girders  and  a  7.5  inch  thick  concrete  deck.   The  bridge 
has  a  roadway  width  of  44  ft.  and  is  designed  for  HS  20-44 
loading. 

Figure  5.1(b)  shows  the  single  beam  with  the  mass 
concentrations  representing  the  tv/o  span  bridge  models 
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Table  5.2    Maximum  Accelerations  of Two  Span  Continuous 
Bridges  with  Different  Span  Lengths 

6  Girder  Bridge    44  ft.  Roadway  Width    n=  4    c/c  -   .02 
Level  Surface    HS  20-44  Loading    60  m.p.h.        ^ 
Smooth  Vehicle 


Span  Length  (ft)  Girder  I*  (in"^)  Wj^l"  (k/ft)  fj^f+Ccps)  T^^  (sec) 

60-60  33WF118   55200      5.46  3.53      0.28 

80-80  36V7F170       82800               5.77  2.36               0.42 

100-100  36WF245  116400      6.22  1.73      0.58 


Maximum  Acceleration  (in/sec^) 
Node 

60-60      x/Lj^  80-80       x/Lj^  100-100       x/L)-, 


1 

-10.20 

0.71 

-15.39 

0.48 

-15.87 

0.50 

2 

24.68 

0.38 

-18.34 

1.15 

-13.94 

1.17 

3 

17.87 

0.38 

-15.47 

0.50 

-12.14 

0.47 

4 

17.86 

0.51 

16.52 

0.44 

11.82 

1.10 

5 

25.86 

1.16 

19.82 

0.43 

-17.12 

0.83 

6 

-21.37 

0.84 

15.88 

0.45 

-14.35 

0.85 

*   I  is  total  moment  of  inertia  of  cross  section. 
+   Wj_,  is  weight  per  linear  length  of  right  span. 
it  fj^  is  fundamental  natural  frequency  of  bridge. 
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described  above.   Also  shown  is  the  vehicle  model  with 
IIS  20-44  loading  and  smooth  moving  condition.   It  is  also 
assumed  that  the  bridge  has  a  level  surface  and  the  vehicle 
travels  at  a  constant  speed  of  60  m.p.h.   Bridge  properties 
for  different  span  lengths  and  the  remaining  vehicle 
parameters  are  listed  in  Tables  5.2  and  4.1  respectively. 
It  should  bo  remembered  that  the  bridge  and  vehicle 
models  shown  are  not  the  complete  models  used  in  the 
analysis.   Figures  2.15  and  2.11  show  the  complete  models. 

Maximum  accelerations  at  the  concentrated  masses 
or  nodes  are  evaluated  by  using  bridge  models  with  different 
span  lengths  and  are  listed  in  Table  5.2.   The  corresponding 
position  of  the  front  axle  as  the  maximum  acceleration 
occurred  is  also  shown.   These  results  are  p]otted  in 
F'igure  5.2.   It  can  be  seen  that  the  accelerations  increase 
for  the  shorter  spans.   This  is  in  agreement  with  the 
results  of  simple  span  and  three  span  continuous  bridges 
in  Articles  3.3.1  and  4.3.1  respectively.   It  can  be  noted 
that  the  magnitudes  of  accelerations  of  the  two  span 
bridges  are  about  1.5  times  larger  than  those  of  three 
span  bridges  for  equal  span  length.   In  order  to  compare 
the  accelerations  of  the  two  span  bridges  with  those  of  the 
simple  :5pan  bridges,  the  single  axle  load  is  used  as 
the  vehicle  model  for  the  two  span  bridges.   The  results 
further  discussed  in  Article  5.4.1  indicates  that  the 
accelerations  of  the  two  span  bridges  due  to  single  axle 
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Figure  5.2    Effect  of  Span  Length  on  Acceleration 
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load  are  1.5  times  larger  than  those  due  to  three  axle 
vehicle.   For  equal  span  length  and  the  same  single  axle 
load  as  the  vehicle  model,  it  can  be  seen  that  the  accelera- 
tions of  the  tv;o  span  bridges  are  still  less  than  those  of 
the  simple  span  bridges. 

It  is  also  noted  that  the  accelerations  of  the  two 
midspan  nodes  are  greater  than  those  of  other  nodes . 
Figure  5.3  shows  the  history  curves  of  the  accelerations 
of  the  two  midspan  nodes.   The  periods  of  both  accelerations 
correspond  to  the  lowest  three  naturiil  periods  of  vibration 
of  the  bridge.   The  maximum  accelerations  of  nodes  2  and 
5  occur  when  the  center  of  gravity  of  the  total  load  is 
either  on  the  mid-left  span  or  on  the  mid-right  span. 

5.3.2   Effect  of  Girder  Flexibility 

Figure  5.1(c)  shows  the  bridge  model  used  for  the 
study  of  accelerations  as  affected  by  using  different 
values  of.    ciirder  stiffness.   The  bridge  has  two  60  ft. 
spans  with  a  damping  factor  of  0.02  and  a  level  bridge 
surface.   The  remaining  parameters  are  considered  to  be 
the  same  as  before.   Four  different  girder  sections 
are  used  to  evaluate  the  acceleration  at  each  node. 
The  properties  of  each  section  were  obtained  from  the 
information  given  in  Reference  42  and  listed  in  Table  5.3. 
The  33WF118  is  supposed  to  be  the  proper  design  section 
of  tliis  bridge  model. 
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Table  5.3   Maximum  Accelerations  of  Two  Span 

Continuous  Bridge  Obtained  by  Using 
Different  Girder  Flexibilities 

60-60  Span  Bridge    n=4    c/c^.  =  .02    Level  Surface 
HS  20-44  Loading    60  m.p.h.    Smooth  Vehicle 


Girder    1^^   (in^)  w  (k/ft)   A  (in2)    f^^  (cps)    Tj^  (sec) 


33WF118 

5900 

.118 

34.8 

3.53 

0.28 

30WF116 

4930 

.116 

34.2 

3.34 

0.30 

30WF99 

4000 

.099 

29.1 

3.18 

0.31 

2  7WF84 

2830 

.084 

24.8 

2.94 

0.34 

Node 


Maximum  Acceleration  (in/sec^) 


33WF118   x/Lj^   30WF116  x/L^       30WF99  x/L^      27WF84   x/L^ 


-18.20       0.71  15.81       0.40         21.59       0.40         22.54       0.40 


2 

24 

63 

0 

38 

21 

23 

0 

39 

24 

51 

0 

39 

21 

70 

0 

39 

3 

17 

87 

0 

38 

15 

84 

0 

38 

20 

90 

0 

38 

-23 

06 

0 

55 

4 

17 

86 

0 

51 

-19 

39 

1 

11 

22 

87 

0 

53 

28 

78 

0 

54 

5 

25 

86 

1 

16 

-22 

96 

0 

86 

-26 

79 

0. 

88 

26 

99 

0 

53 

6 

-21 

37 

0 

84 

-20 

99 

0 

87 

-23 

99 

0 

88 

25 

75 

0, 

51 
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Maximum  accelerations  at  nodes  evaluated  by  using 
the  four  different  girder  sections  are  listed  also  in 
Table  5.3  and  are  plotted  as  a  function  of  girder  moment 
of  inertia  in  Figure  5.4.   As  might  be  expected,  the 
results  show  that  the  increased  flexibility  of  girder 
does  not  affect  the  overall  acceleration  level  of  the  bridge 
appreciably  although  the  moment  of  inertia  of  girder 
v;as  reduced  more  than  half  of  the  original  value. 

5 . 4   Vehicle  Parameters 

5.4.1   Effect  of  Number  of  Axles  and  Axle  Spacing 

Figure  5.5  shows  the  typical  vehicle  models  with 
different  number  of  axles.   The  details  describing  the 
parameters  of  typical  vehicles  were  given  in  Table  4.1. 
Ihe  bridge  model  and  the  remaining  parameters  are  the 
same  as  in  Article  5.3.2.   All  typical  vehicles  are 
assumed  to  be  smoothly  moving  with  a  constant  speed  of 
60  m.p.h.   The  maximum  accelerations  at  nodes  subjected 
to  different  typical  vehicles  are  compared  in  Table  5.4 
and  shown  in  Figure  5.6.   As  might  be  expected,  the  magni- 
tudes of  acceleration  for  two  axle  and  three  axle  typically 
are  about  the  same  but  they  are  about  two-tliirds  of  the 
magnitudes  for  single  axle  loads. 

Table  5.5  and  Figure  5.7  show  the  effect  of  vehicle 
axle  spacing  on  bridge  accelerations.  The  vehicle  model 
has  the  IIS  20-44  loading  with  the  practical  values  of  axle 
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Table  5.4    Maximum  Accelerations  of  Two  Span 

Continuous  Bridges  Subjected  to  Different 
Numbers  of  Axles  of  Vehicle 

60-60  Span  Bridge    33WF118  Girders    n=4    c/c  =  .02 
Level  Surface    Smooth  Vehicles    60  m.p.h.      ^ 


Node 


Maximum  Acceleration  (in/sec^) 


3-Axle      x/Lj3  2-Axle      x/L^  1-Axle      x/L^ 


1 

-24.41 

0.46 

-24.31 

0.36 

-28.73 

0,19 

2 

-23.15 

0.44 

-23.80 

0.34 

31.75 

0.04 

3 

-21.58 

0.43 

-21.11 

0.33 

29.57 

0.  03 

4 

-17.81 

0.35 

17.73 

0.32 

26.43 

0.  16 

5 

17.68 

1.13 

19.26 

1.03 

30.89 

0.15 

6 

-18.60 

0.79 

-17.74 

0.69 

28.88 

0.14 

Table  5.5   Maximum  Accelerations  of  Two  Span 

Continuous  Bridges  Subjected  to  Vehicles 
with  Different  Axle  Spacings 

60-60  Span  Bridge    33V7F118  Girders   n=4    c/c  =.02 
Level  Surface    72  kips  Total  Load    3-Axle     ^ 
60  m.p.h.    Smooth  Vehicle 


Node 

Maximum 

Acceleration    ( 

in/sec 

2) 

l2=  20  ft 

x/I^ 

l2  =  25  ft 

x/I^ 

I2  =  30  ft 

x/Lj, 

l2=35  ft 

y. 

■■/h, 

1 

-24.41 

0.46 

-23.18 

1.17 

-18.20 

0.71 

15.93 

n, 

,40 

2 

-23.15 

0.44 

28.54 

1.25 

24.68 

0.38 

15.52 

0, 

.41 

3 

-21.58 

0.43 

-21.81 

0.48 

17.87 

0.38 

-13.14 

1  , 

,  1  9 

4 

-17.81 

0.35 

20.16 

1.13 

17.86 

0.51 

-14.24 

1  . 

,  1  1 

5 

17.68 

1.13 

28.68 

1.14 

25.86 

1.16 

17.23 

n. 

54 

6 

-18.60 

0.79 

-22.02 

0.83 

-21.37 

0.84 

-17.22 

1. 

29 
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Figure    5.5        Vehicle   Models   with   Different  Numbers 
of  Axles 
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spacing  of  trailer  varied  from  20  ft.  to  35  ft.   The  axle 
spacincj  of  the  tractor  is  kept  constant  at  12  ft.    All 
remaining  parameters  are  considered  to  be  the  saime  as 
before.   From  the  results  shown,  the  magnitudes  of 
accelerations  are  largest  when  the  axle  spacing  is  equal 
to  25  ft.,  or  in  non-dimensional  form,  the  axle  spacing 
ratio  has  the  value  of  0.42.   The  axle  spacing  ratio  in 
this  case  is  the  ratio  of  the  axle  spacing  to  the  length 
of  eitlior  span  of  the  two  span  bridge.   The  value  of 
0.42  is  in  good  agreement  with  the  results  presented 
in  Article  4.4.2,  where  three  span  bridge  has  maximum 
accelerations  for  an  axle  spacing  ratio  between  0.37  to 
0.43. 

5.4.2   Effect  of  Vehicle  Speed 

Table  5.6  and  Figure  5.8  show  the  effect  of  vehicle 
speed  on  acceleration  of  the  two  span  continuous  bridge. 
The  bridge  m.odel  was  shovm  in  Figure  5.1(c)  and  has  the 
same  parameters  as  before.   The  vehicle  model  has  HS  20-44 
loading  and  is  smoothly  moving.   The  vehicle  speed  ranges 
from  30  m.p.h.  to  70  m.p.h.   Generally,  the  accelerations 
of  the  bridge  increase  for  increased  values  of  speed. 
However,  the  accelerations  of  nodes  1  and  6  decrease  for 
the  speed  of  50  m.p.h.   It  is  believed  at  this  stage  that 
they  v;ere  influenced  by  the  other  parameters  such  as  the 
period  of  oscillation  of  interacting  force  and  the 
period  of  the  bridge  itself. 
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5.4.3   Effect  of  Frequency  Ratio 

The  maximum  accelerations  evaluated  by  using  the 
different  values  of  frequency  ratio  are  shown  in  Table 
5.7  and  Figure  5.9.   The  frequency  ratio  is  the  ratio  of 
vehicle  frequency  to  the  natural  frequency  of  the  bridge. 
In  this  study,  the  frequency  ratio  is  considered  to  be  in 
the  range  between  0.5  and  2.0,  the  bridge  model  used  is  two 
60  ft.  spans  and  has  a  fundamental  natural  frequency  of 
3.53  cps .   The  vehicle  frequency  is  the  tire  frequency 
of  rear  axles.   It  is  assumed  that  the  vehicle  oscillates 
on  its  tires  only  and  the  coefficients  of  friction  of 
suspension  springs  have  the  value  of  infinity.   The  vehicle 
has  HS  20-44  loading,  and  is  smoothly  moving  at  a  constant 
speed  of  60  m.p.h.   Generally,  the  magnitudes  of  accelera- 
tions at  the  nodes  are  about  the  same  for  all  values  of 
frequency  ratio  although  the  accelerations  of  midspan 
nodes  for  f  /f,  =  1.0  increase  slightly.   It  may  be  con- 
cluded that  the  accelerations  of  the  bridge  are  not 
materially  affected  by  the  frequency  ratio. 

5 . 5   Initial  Conditions  of  Vehicle  and  Bridge 

5.5.1   Effect  of  Initially  Oscillating  Vehicle 

As  previously  explained  in  Article  4.5.1,  the  vehicle 
usually  has  an  initial  oscillation  due  to  approach  pavement 
irregularities  and  a  possible  discontinuity  at  the  abutment, 
The  parameters  representing  the  initial  oscillation  are 
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classified  as  the  amplitude  of  initial  force  variation  C., 
the  initial  phase  angle  G.,  and  the  initial  coefficient 
of  friction  ii  .  .   The  results  of  the  three  span  bridges 
in  Article  4.5.1  showed  that  the  bridge  acceleration 
v/as  not  affected  by  initial  phase  angles  and  initial 
coefficients  of  friction.   It  is  believed  that  the  accelera- 
tions of  the  two  span  bridges  are  also  not  affected  by 
these  two  parameters  because  the  method  of  analysis  for 
the  two  span  bridges  is  similar  to  that  for  three 
span  bridges.   Therefore,  the  amplitude  of  the  initial 
force  variation  C.,  is  the  only  parameter  studied  here. 

As  in  the  study  of  the  three  span  bridges,  the 
amplitudes  of  initial  force  variation  C,  are  assigned 
to  have  the  values  of  0.0,  0.15,  0.30  and  0.50.   The  maximum 
accelerations  at  the  nodes  obtained  by  using  different 

values  of  C.  are  shown  in  Tabel  5.8  and  Figure  5.10.   The 
1 

bridge  is  the  same  as  before  and  the  vehicle  model  has 
IIS  20-44  loading  v;ith  initial  oscillation  and  a  constant 
speed  of  60  m.p.h.   The  values  of  C.  are  considered  to 
bo  the  same  for  each  axle.   The  initial  values  of  the 
frictional  force  and  phase  angle  for  each  axle  are  all 
taken  equal  to  zero. 

It  can  be  seen  from  the  results  that  the  maximum 
accelerations  increase  for  the  larger  amplitude  of  initial 
oscillation.   For  the  values  of  C.  equal  to  0.15  and  0.50, 
the.  magnitude  of  maximum  acceleration  of  all  nodes  increases 
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respectively  1.3  and  1.5  tirpes  larger  than  those  for 

C.  =  0. 

1 

Table  5.9  and  Figure  5.11  show  the  effects  of  initial 
phase  angle  difference  of  the  50%  initially  oscillating 
vehicle  on  bridge  accelerations.   It  is  assumed  that  the 
vehicle  has  ©1=02=  0°  and  0   ranges  from  0°  to  180". 
All  remaining  parameters  are  the  same  as  before.   The 
results  show  that  the  magnitudes  of  accelerations  are 
about  the  same  for  all  values  of  phase  angle  difference 
except  the  accelerations  of  some  nodes  are  larger  for 
the  values  of  phase  angle  difference  equal  to  90"  and 
150°,  due  to  the  behavior  of  interacting  forces.   It 
can  be  noted  that  the  magnitudes  are  still  below  50  in/sec' 
which  docs  not  affect  the  human  response  according  to 
Reference  5. 

5,5.2   Effect  of  Surface  Roughness  of  Bridge 

As  previously  explained  in  Article  4.5.2,  the  number 
of  half  sine  waves  and  the  roughness  amplitudes  are  the 
two  parameters  assumed  to  represent  the  surface  roughness 
of  higliway  bridge.   In  Figure  5,12  are  shown  the 
profiles  of  bridge  models  having  one,  two  and  four 
half  sine  waves.   The  peak  amplitudes  of  half  sine 
waves  are  kept  constant  at  0.5  in.   The  bridge 
model  is  6  0  ft.  -  60  ft.  span  bridge.   The  vehicle  model 
has  HS  20-44  loading  and  is  smoothly  moving  at  a  speed 
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Table  5.9    Maximum  Accelerations  of  Two  Span 

Continuous  Bridges  Subjected  to  50% 
Initially  Oscillating  Vehicle  with 
Varied  Phase  Angle  Differences 

60-60  Span  Bridge    33WF118  Girders    n=4    c/Cj.=  .02 
Level  Surface    HS  20-44  Loading    60  m.p.h. 
Oscillating  Vehicle    0^^=0.50    p  j_  =  0    0^^  =  O2  =  0° 
A0  =  0o-Ot 


Maximum 

Acceleration 

(in/sec 

2) 

Node 

AO 

=  0° 

x/Lj^ 

AG  =  30° 

x/Lj, 

AO  =  60° 

x/Lb 

AG  = 

90° 

x/I^ 

1 

32. 

.72 

1.29 

-33.74 

0.54 

30.97 

1.29 

33. 

,24 

0.62 

2 

38. 

.79 

1.30 

-41.03 

0.52 

38.50 

1.30 

40. 

.50 

1.27 

3  ■ 

-34. 

,21 

0.52 

-39.26 

0.52 

-35.87 

0.52 

-32. 

,39 

1.19 

4 

36. 

.18 

0.51 

-34.67 

0.40 

-34.88 

0.40 

-36. 

,54 

0.40 

5 

42. 

.25 

0.50 

-40.44 

0.41 

41.51 

0.50 

49. 

,94 

0.50 

6 

41. 

.00 

0.49 

-40.34 

0.42 

-40.50 

0.42 

48. 

,00 

0.49 

Maximum  Acceleration  (in/sec^) 


Node 


A0=120°    X/L5       AO  =  150°    x/L]3    AO  =  180°    x/Lj^ 


1  29.98  1.29  -49.77  0.54       46.55  0.62 

2  38.42  1.30  49.40  0.61       45.39  0.63 

3  -37.52  0.52  38.32  0.60       33.25  0.38 

4  -34.79  0.40  -36.27  0.40    -37.40  0.40 

5  -40.71  0.41  -42.10  0.41       45.23  0.50 

6  -40.45  0.42  45.90  0.49       47.28  0.49 
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or  60  m.p.h.   All  remaining  parameters  are  the  same  as 
before. 

The  maximum  accelerations  at  the  nodes  are  evaluated 
by  using  the  number  of  half  sine  waves  varied  from  0  to 
10.   The  results  are  listed  in  Table  5.10  and  shown 
graphically  in  Figure  5.13.   The  results  are  similar  to 
those  of  the  three  span  bridge  in  such  a  way  that  the 
accelerations  are  maximum  for  the  value  of  4  half  sine 
waves.   It  can  be  explained  in  Figure  5.14-1  and  5.14-2 
v;hich  shows  the  histories  of  accelerations  at  mid  right 
span  and  of  interacting  forces  of  the  second  axle 
respectively.   In  each  figure,  the  results  obtained  by 
using  the  level  surface  condition  and  four  half  sine 
waves  are  compared. 

It  can  be  seen  that  the  interacting  force  of  the 
vehicle  is  influenced  by  the  surface  roughness .   Its 
magnitude  and  period  of  acceleration  have  been  changed. 
The  period  of  oscillation  of  the  interacting  force  for  the 
value  of  four  half  since  waves  is  coincident  with  the 
fundamental  natural  period  of  bridge.   Therefore,  resonance 
between  the  bridge  and  the  vehicle  occurs.   The  magnitudes 

of  accelerations  are  much  larger  than  the  suggested 

2 

value  of  100  in/sec   and  this  result  can  cause  the 

unpleasant  condition  to  the  bridge  users. 

Table  5.11  and  Figure  5.15  show  the  comparison  of 
the  effect  of  varied  amplitude  of  roughness  on  acceleration 
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Table  5.10    Maximum  Accelerations  of  Two  Span 

Continuous  Bridges  Due  to  Different 
Numbers  of  Half  Sine  Waves 

60-60  Span  Bridge    33WF118  Girders    n=4    c/Cj,  =  .02 
HS  20-44  Loading    60  m.p.h.    Smooth  Vehicle 
Amplitude  of  Roughness =  0 . 5  inch 


Maximum  Accelerat 

ion  (in/sec^) 

Node 

NS* 

=  0 

x/Lb 

NS 

=  1 

x/Lj^ 

NS  =  2 

x/Lb 

NS  =  4 

x/Lb 

1 

-18, 

.20 

0.71 

74, 

.93 

1.29 

-62.47 

1.12 

261.57 

1.30 

2 

24. 

.68 

0.  38 

74. 

.41 

1.30 

-70.74 

0.86 

-298.01 

0.98 

3 

17. 

.87 

0.  38 

70, 

.50 

1.31 

-65.16 

1.10 

169.64 

1.05 

4 

17. 

.86 

0.51 

73. 

.51 

1.18 

-74.12 

0.97 

-213.59 

0.86 

5 

25. 

.86 

1.16 

85. 

.16 

1.18 

-99.27 

0.98 

353.05 

1.17 

G 

-21. 

.37 

0.34 

73. 

.95 

1.16 

-98.47 

0.99 

291.74 

1.17 

Maximum  Acceleration  (in/sec2) 


Node 


NS    =    6       x/Ly.         NS  =  8       x/Ly.         NS=10       x/Ly 


1 

162, 

.59 

1. 

.26 

69. 

,01 

0. 

,37 

112. 

,10 

0. 

,69 

2 

191. 

.39 

1. 

.25 

-94. 

.69 

0, 

.49 

-59. 

.02 

0. 

,47 

3 

-122. 

.79 

0. 

.51 

-61. 

.43 

0. 

.48 

-92. 

.72 

0. 

.52 

4 

122. 

.23 

1. 

.24 

63. 

.  38 

1. 

,19 

101. 

,45 

1. 

.05 

5 

-224, 

.94 

1. 

.31 

90, 

.25 

1. 

,03 

78. 

.38 

1. 

.03 

188.14   1.32    66.15   1.02    113.47   1.02 


*  NS  =  Number  of  Half  Sine  Waves 
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for  tv;o  different  cases  of  number  of  half  sine  waves. 
In  case  A,  two  half  sine  waves  are  selected,  which 
slightly  affects  the  acceleration  of  bridge,  while  in 
case  B,  the  four  half  sine  waves  are  used,  which  has 
the  most  effect.   In  both  cases,  the  nodal  accelerations 
are  obtained  by  using  the  amplitude  of  roughness  in  the 
range  from  0.0  to  1.0  in.   It  can  be  seen  from  the 
results  that  the  maximum  accelerations  increase  for  the 
larger  value  of  roughness  amplitude.   The  amount  of 
increase  in  magnitudes  of  acceleration  for  four  half 
sine  waves  is  larger  than  that  of  two  half  sine  waves. 
It  may  be  concluded  that  the  amplitvide  of  roughness  affects 
the  acceleration  only  if  the  number  of  half  sine  waves  does 

One  can  see  that  these  results  for  two  span  bridge 
are  similar  to  those  of  simple  span  and  three  span 
bridges.   The  only  difference  among  the  three  is  the 
magnitude.   The  magnitude  of  acceleration  is  largest 
for  simple  span  bridge  and  least  for  three  span  bridge. 
It  is  also  seen  that  the  effect  of  surface  roughness  seems 
to  be  the  most  important  variable  in  terms  of  bridge 
accelerations. 
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CHAPTER  VI 
SUMMARY  AND  CONCLUSIONS 

6  . 1   Summary 

The  effects  of  the  major  parameters  on  the  accelera- 
tions of  liighway  bridges  have  been  investigated  and  compared 
to  the  acceleration  criteria  for  human  response.   Three 
different  types  of  highway  bridges  were  investigated: 
simple  span,  two  span  continuous  and  three  span  continuous 
bridges.   Major  parameters  selected  include  the  parameters 
reflated  to  the  bridge,  to  the  vehicle  and  the  initial 
conditions  of  bridge  and  vehicle.   The  numerical  data 
pr(?sented  are  derived  from  a  theory  in  which  the  bridge 
is  idealized  as  a  plate  continuous  over  flexible  beams 
for  simple  span  bridges  and  as  a  continuous  beam  with 
concentrated  point  masses  for  multi-span  bridges. 
Tlic  vehicle  has  been  represented  as  a  sprung  single  load 
consisting  of  tv;o  wheels  for  the  simple  span  bridge 
model  and  as  a  sprung  load  unit  having  one,  two  or  three 
axles  with  a  suspension  system  for  multi-span  bridges. 

Most  of  the  bridge  parameters  used  were  obtained 

(41) 
from  the  information  of  Standard  B.P.R.  bridges      while 

the  values  of  vehicle  parameters  came  from  manufacturers' 

(2)  (22)  (29) 
data  and  test  reports.  Because  of  the  very 
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large  number  of  variables  involved  and  the  considerable 
computer  time  required  for  a  solution,  it  was  impractical 
to  obtain  solutions  for  all  possible  combinations  of  the 
variables.   For  the  above  reason,  one  parameter  was  varied 
at  a  time  to  investigate  its  effect.   The  results  of  the 
investigation  are  shown  in  the  form  of  history  curves 
of  the  response,  tables  and  graphs. 

The  principal  results  of  the  study  of  simple  span 
bridge  are  summarized  as  follov;s: 

a.  Standard  B.P.R.  bridges  with  different  span 
lengths  were  investigated  and  the  maximum  accelerations  of 
each  beam  was  obtained.   The  results  show  that  the 
acceleration  of  the  bridge  decreases  for  l(3nger  span 
lengths  and  the  accelerations  of  the  exterior  beams 

tend  to  be  larger  than  those  of  interior  beams. 

b.  For  two  different  cases  of  the  transverse  position 
of  load  on  the  bridge,  over  the  edge  beam  and  over  the 
c-^nter  beam,  the  width  variations  of  the  bridge  model  do 
not  affect  the  accelerations  appreciably. 

c.  The  increased  flexibility  of  beams  does  not  have 
much  effect  on  the  accelerations.   The  results  show  that 
the  maximum  acceleration  increases  approximately  by  20 
percent  when  the  moment  of  inertia  of  beam  is  reduced 
nearly  40  percent. 
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d.  The  variations  of  the  transverse  position  of  load 
have  an  effect  on  the  acceleration  distribution  among  the 
beams.   Generally,  the  acceleration  of  an  edge  beam  is 
greater  than  that  of  an  interior  beam. 

e.  The  bridge  acceleration  obtained  by  using  a  single 
wheel  load  is  approximately  10  percent  greater  than  that 

by  using  two  wheels,  regardless  of  the  transverse  position 
of  load  on  the  bridge. 

f.  By  using  the  combination  of  two  variables,  the 
flexibility  of  beam  and  transverse  position  of  a  single 
J.oad,  it  was  confirmed  that  the  increased  flexibility  of 
girder  did  not  have  much  effect  on  accelerations  as  long 
as  the  moment  of  inertia  of  the  beam  section  was  not 
reduced  more  than  35  percent  of  the  standard  size. 

g.  The  speed  of  vehicle  affected  the  acceleration 
of  the  bridge.   The  accelerations  were  greater  for 

the  higher  speeds. 

h.   The  surface  roughness  of  the  bridge  was  idealized 
as  a  series  of  continuous  half  sine  waves  along  the  span 
length  with  equal  amplitude  for  each  wave.   The  accelerations 
of  the  bridge  as  affected  by  the  variation  of  the  number 
of  half  sine  waves  and  also  the  amplitudes  of  roughness 
v;ere  studied.   Generally,  the  acceleration  of  the  bridge 
with  a  rough  surface  was  greater  than  that  of  the  bridge 
with  a  level  surface  condition.   It  was  shown  that  the 
effect  of  the  number  of  half  sine  waves  caused  the  bridge 
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to  vibrate  severely  if  the  period  of  the  interacting  force 
as  influenced  by  the  roughness  was  equal  to  one  of  the 
natural  periods  of  bridge.   The  results  also  showed  that 
the  amplitude  of  roughness  had  an  effect  on 
acceleration  only  if  the  number  of  half  sine  waves  did. 
Finally,  the  acceleration  increased  for  larger  roughness 
amplitude .   ■ 
.     i.   By  comparing  the  amplitudes  of  acceleration 

obtained  in  this  study  to  the  suggested  value  of  100 

2 
in/sec   as  the  criteria  for  pedestrian  discomfort.   It 

was  shown  that  the  maximum  accelerations  of  the  simple 

span  bridges  were  generally  larger  than  the  suggested  value, 

For  the  study  of  three  span  and  two  span  continuous 
highway  bridges,  the  following  is  a  brief  summary  of  the 
principal  results: 

a.   Throughout  the  study,  the  spans  were  kept  in  the 
ratios  of  4:5:4  for  three  span  bridges  and  equal  spans 
for  two  span  bridges.   Both  bridge  types  had  a  damping 
factor  of  0.02  and  three  concentrated  masses  on  each  span. 
The  results  showed  that  the  accelerations  increased  for 
s^iortcr  span  lengths.   Additionally,  the  accelerations 
at  the  middle  of  any  span  were  more  important  and  had 
greater  amplitude  than  the  accelerations  at  other 
points  of  that  span. 
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b.  The  increased  girder  flexibility,  the  frequency 
ratio  and  the  span  ratio  did  not  materially  affect  the 
bridge  accelerations.     i  ■      ,  • 

c.  Generally,  the  bridge  accelerations  were  greater 
for  the  heavier  vehicles.   For  the  two  vehicles  with  the 
same  total  load,  the  one  which  had  fewer  axles  caused 
greater  accelerations.   The  axle  spacing  of  trailer  for 
three  axle  vehicle  did  not  have  much  effect  on 
accelerations  although  the  acceleration  was  greatest  for 
the  values  of  axle  spacing  ratio  ranging  between  0.37 

to  0.43. 

d.  The  vehicle  speed  did  have  an  effect  on  accelera- 
tions, with  higher  speeds,  causing  greater  bridge 
accelerations. 

G.   The  effect  of  initial  vehicle  oscillations  due 
to  the  bridge  surface  irregularities  or  the  sharp  discontin- 
uity at  the  approach  pavement  v;as  investigated.   The  initial 
oscillations  included  the  amplitude  of  initial  force 
variation,  the  phase  angle  difference  between  the  axles 
and  the  initial  value  of  the  frictional  force. 

Generally,  the  bridge  accelerations  when  subjected 
to  the  initially  oscillating  vehicle  were  about  1.5  times 
larger  than  the  accelerations  of  the  bridge  ^^7hen  subjected 
to  the  smooth  vehicle.   The  higher  the  amplitude  of 
irregularities  at  the  approach,  the  greater  was  the  accelera- 
tion.  The  initial  phase  angle  difference  and  the  initial 
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frictional  force  did  not  have  much  effect  on  the 
acceleration.  '  • 

f.  The  surface  roughness  of  bridge  was  idealized 
as  a  series  of  continuous  half  sine  waves  in  order  to 
compare  the  results  to  those  obtained  for  simple  span 
bridges.   Each  span  was  considered  to  have  an  equal  number 
of  half  sine  waves.   The  accelerations  of  the  bridge 

with  a  rough  surface  were  greater  than  those  of  the  bridge 
with  a  level  surface  condition.   The  bridge  can  severely 
vibrate  when  the  period  of  the  interacting  force  influenced 
by  the  roughness  is  equal  to  one  of  the  natural  periods 
of  vibration  of  the  bridge.   Finally,  the  amplitude  of 
roughness  had  an  effect  on  acceleration  only  if  the  number 
of  half  sine  waves  did  and  the  acceleration  increased 
for  larger  roughness  amplitude. 

g.  The  magnitudes  of  accelerations  as  affected  by 

bridge  and  vehicle  parameters  for  a  level  surface 

2 
condition  were  less  than  the  suggested  value  of  100  in/sec  . 

The  only  parameter  having  an  effect  approaching  the  criteria 

for  human  response  for  multi-span  highway  bridges  was  the 

bridge  surface  roughness. 

6.2   Conclusions 

For  simple  span  bridges,  the  amplitudes  of  accelera- 
tions which  psychologically  disturb  the  pedestrian  were 
predominately  affected  by  the  span  length  of  the  bridge. 
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the  weight  and  speed  of  the  vehicle,  and  the  surface  rough- 
ness.  The  parameters  that  provided  the  minor  effects 
were  the  girder  flexibilities  and  the  transverse  position 
of  load.   For  two  and  three  span  continuous  bridges,  the 
magnitudes  of  accelerations  were  larger  than  the  recommended 
limit  of  comfort  only  when  the  surface  roughness  of  the 
bridge  was  taken  into  account.   High  strength  steel  girders 
could  therefore  be  used  for  highway  bridges  since  the 
effect  on  the  bridge  accelerations  was  relatively 
insignificant . 

From  the  information  presented  in  this  study,  a 
basis  for  design  regulating  the  satisfactory  vibration 
characteristics  for  human  comfort  can  possibly  be 
established.   The  design  procedure  would  appear  in  the  form 
of  estimating  the  values  of  maximum  accelerations 
produced  under  the  most  unfavorable  but  likely  combinations 
of  parameters  involved.   The  major  parameters  such  as  the 
speed  of  vehicle,  the  span  length  and  the  natural  frequency 
of  vibration  of  bridge  were  combined  in  the  nondimensional 
form  of  speed  parameter  a.   This  parameter  v;as  expressed 
in  the  form: 

vTj^ 

"  '-^   -217 

where  v  is  the  speed  of  vehicle,  T.  is  the  fundamental 

b 

period  of  vibration  of  bridge  and  L  was  either  the  length 
of  the  simple  span  bridge,  or  the  length  of  center  span  of 
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the  three  span  continuous  bridge,  or  the  length  of  either 
span  of  the  two  span  continuous  bridge.   Some  examples  of 
the  prediction  of  acceleration  of  highway  bridges  as  a 
function  of  speed  parameter  and  weight  ratio  are  shown 
in  Figures  6.1  and  6.2.   The  weight  ratio  is  the  ratio 
of  the  weight  of  vehicle  to  the  weight  of  bridge. 

Further  study  is  recommended  to  investigate  the 
effect  of  surface  roughness  which  was  the  most  important 
parameter  in  producing  an  adverse  effect  on  human  response. 
The  measured  profiles  of  various  bridges  in  the  field  are 
needed  for  better  correlation  in  this  regard.   A  basis 
for  design  to  control  the  vibration  characteristics 
within  the  comfort  level  should  also  be  further  studied 
and  proposed  for  bridge  design  specifications.   Other 

types  of  highway  bridges  should  be  investigated,  such 

..      .  -,      .       •  ^      ^  ^     (14)  (18)  (20) 

as  the  cantilever-type,  since  many  test  reports 

show  that  this  type  produces  the  most  adverse  effects. 
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5   GIRDER   BRIDGE      SIDE  RflTI0=0.8     FREQUENCY   RflT10=l.Q 
TWO  WHEELS  OVER   BEAM   1      LEVEL   SURFACE 


Figure    6.1        Acceleration   of    Simple   Span   Bridge  with 
Level   Surface 
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5  GIRDER  BRIDGE     SIDE   RflTIQzQ.S     FREQUENCY   RflTia=1.0 
TWO  HHEELS   OVER  BEAM   1      ROUGH  SURFACE 


Figure    6.2        Acceleration   of    Simple   Span   Bridge  with 
Rough   Surface 
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